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INTRODUCTION 

9 

This manual describes in detail the operating and 
set-up procedures for. the multi-satellite multi-arc GEODYN 
Orbit Determination Program. Additionally, all system output 
is described. 

*• n . 

The GEODYN Program is the nucleus of the entire 
GEODYN System. ‘ It is a definitive orbit and geodetic 
parameter estimation program capable of simultaneously 
processing observations from multiple arcs of multiple 
satellites. 

GEODYN has two modes of operation: the data 

reduction mode and the orbit generation mode. 

• In the data reductiom mode the following para- 
meters can be estimated from satellite tracking data: 

• Six orbital parameters representing the 
inertial position and velocity of- the 
satellite at some given epoch, 

• Physical constants related to the atmos- 
pheric drag and solar radiation forces 
affecting the motion of single satellites, 

t Tracking station coordinates relative to 

the earth 1 s center of mass, 

• Tracking instrument errors, a zero-set 
bias, and a timing bias, and 

o Geopotential gravity model coefficients. 

• Surface density parameters. 



» ’ 


In the data reduction mode the six orbital para- 
meters are always estimated for each arc; if requested, the 
remaining parameters are estimated as options. If accurate 
values for the physical constants or geopotential coefficients 
are known a priori , these can be used in the system as fixed 

values. 

In the orbit generation mode, the program will, given 
all of the orbit parameters a priori , compute satellite positions 
and velocities at times later than and/or prior to the epoch. 

The GEODYN Program uses the terms inner and outer 
iterations. On inner iterations only individual arc 
parameters may be estimated. Individual arc parameters 
.include the six orbital elements, drag, solar radiation, 
geopotential resonance coefficients, instrument biases and 
tracking station timing biases. 

On outer iterations geodetic parameters are estimated 
using data from all arcs. Geodetic parameters include tracking 
station coordinates, geopotential gravity coefficients, and 
♦ surface densities. 


11 



SECTION 1.0 ' ' : 

• CEODYN INPUT 'CARDS 

Specific problems to be solved by the GEODYN program 
are defined by input cards, which are separated into two 
categories: 

1* COMMON SET - These are input cards containing 
information that is common to all of the arcs 
being processed in the problem. 

2. ARC SET - A collection of ARC SET cards is 
necessary to define each individual arc. 

t 

Thus the GEODYN input deck will consist of one 
Collection of cards from COMMON SET and a collection of 
cards from ARC SET for each arc. 

Each set contains mandatory and optional cards. 

Mandatory Cards - These are input cards that must 
be present for each problem. 

Optional Cards - These are input cards that are 
optionally used. 



1.1 


COMMON SET 




COMMON SET defines that information which is common 
to all arcs and is necessary for the problem description. 

Included in COMMON SET are four groups to be arranged 
in the following order: 

. . 1* Run description. (Mandatory) 

2. Option cards which may be used to exercise 
GEODYN’s capabilities to estimate tracking 
station positions, modify the earth force 
model, estimate geopotential coefficients, 
estimate surface densities, update the atmos- 
pheric density function, and modify the lunar 
and/or planetary solar gravity effects. (Optional) 

3. Station coordinate cards which may be used 
to alter GEODYN's built-in complement of 
tracking stations. (Optional) 

4. COMMON SET termination card. (Mandatory) 

<m 

The following pages describe these groups more 
thoroughly. 


SSEKF-V 
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1.1.1 COMMON SET 

GROUP 1 - RUN DESCRIPTION • • „ 

9 

TITLE CARDS 



6EOPOTENTIAL RECOVERY RUN SOLVING FOR 2 PAIRS OF COEFFICIENTS 
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COLUMNS FORMAT 


DESCRIPTION 


1-80 20A4 




Run description to be printed 
on title page; use three title 
cards at the beginning of 
each. run. 


Three title cards are Mandatory . 


1 . 1-2 




1 . 1.2 


August 13, 1973 


COMMON SET 

GROUP 2 - OPTION CARDS 


s 


Each Group 2 card is identified by the name 
beginning in column 1 and read under an A6 format. 
The name is followed by 9 fields: 


Field Columns 

Field Format 

7, 8, 9, 10 

411 

11-25 

D15.8 

26-40 

D15.8 

41-55 

D15.8 

56-70 

D15. 8 

71-80 

DIO. 8 

The Option Cards are: 

• - 

BODIES 

RSTART 

CORREL 

SEQDAT 

. . EARTH 

STAEST 

EDIT 

STEP 

ELCUT 

SURF 

FLUX 

TIDES 

RECOEF (COEF) 

TOLS 

RESID 

SURDEN 

BPART 

CENTER 

J65 

J650 

J71 

POLE 


1.1-3 



COMMON SET 

GROUP 2 - OPTION CARDS 
BODIES 


]3*>>£ 




0 0 0 0 
I J 4 S < 

mm 



2222 

r r 

3 3 3 3 3 3 


44’ 444 

: :>Y 

S 5 5 5i 5 


S' 566 6 


317777 


888983 

i 

9 5 9’ 9 5 


o.si 





p ocooooooooooc 

I li 0 M 13 IS I? U n JO 21 21 2224 2: 

tin niiiuin 

{222222222222222 

33 333333333333 

144444444444444 
355555555555555 
i 6 S 6 6 6 S 6 6 S 6 B 6 S £ 
[r 7 7 7 7 ' 7 7 7 7 7 I 7 7 7 

8 ‘ ' 8 3 8 8 8 8 8 8 8 8 2 


1 5 h 9 9 9 S 9 9 9 9 3 9 9 3 9 

r :j 13 14 14 IS II « H JO :t 21 JJ2* 


Vs 3 3 3 - 333 J ; 3 3 i c * ? 3 i i - 1 2 5 3 3 >3 3 ' = 3 3 3 2 3 53 3 3 3 3 3 J 3 3 ? 3 2 2 1 3 ; . _| 
l: . , • i ; 4 4 ; 4 ft 4 4 f 4 4 |&f| : ':t 4 4 4 4 ; g 14 4 4 4 4 4 4 4 5 4 4 4 4« « 4 4 U4 ■ ;j 
' 5 Vc 3 5 i s s s' S : 13 5: i a'5 a 5 : 5 * 3 1 $ 3 £5 5 5 i 5 5 5 i i 5 s 3 5 5 5 3 r . 1 5 5 : a 5 ^ 1 5 * • | 

» 5 5 sts im jWf « :u « : s s ts 2555 3 :im« i *? « * - fi ; = * 5 s s .: j 

K" , : >-r •■•:•'■ ■:: ii".- V-!'-' : ' ■ ■ ' T 7 » 7 7 t 1 • 'i 


; ;im; 175 1; if 1 .; 

's' flf 1 f s' 5 A 2 f : Z0i^ i '.3 i t 2>:S.^ 3 S * 

ja- ;.>V o 

;. a V- ^ »• * ;» - ; • + '*’■ \ * * *. * . . « ri -. . - - •„ -•> -5- '» ;; 


a cf-> - ac iA £ * « v 1 «• 5 ft - i f- V *; ? - :• =■ *+'■ < ♦ * - ‘ 


NAME COLUMNS 

BODIES 1-6 


FORMAT 

A6 

II 


DESCRIPTION 

Modifies non-earth gravity 
perturbations . 

Indicates body for which 
gravity perturbation is to 
be modified: 

1 # - MOON 

2 - SUN 

3 - VENUS 

4 - MARS 

■r fproducibility of T&F 

■ original page is poor 


1.1-4 



BODIES (Cont.) 

NAME COLUMNS FORMAT DESCRIPTION' 

. . 7 (Cont.) 5 - JUPITER 

. . . • ‘ " '6 - SATURN 



II Number * of perturbing 

■ ’• • .bodies to be applied. 

. (Maximum of 6) 




Default - *2 . 

This number may be set on 
any one of the BODIES cards 
. used. 


11-25 D15.8 Value of ratio of mass of 

• f t 

• body to mass of earth. If 

• • • ratio is zero perturbing 

effect will be zero. 
Default mass ratios: 


M./M = 0.01229997171 GO 65C7 

M-/M = 332945.5619264376 

M 3 /M c = 0.81S000322933493S 

M./M - ^.0.1074468525270073 

M,/M = 317 .SS09303 

M-/M„ = 9S. 14905175 

6 e 


IF CARD OMITTED: Lunar and solar gravitational perturbations v:ill 

be applied using default mass ratios unless modified by BODIES 
card in ARC SET. . 

• 

To omit a particular body for gravitational purposes a BODIES 
card representing that body must be included with the correspond 
mass ratio set to zero. 

*For example, if tlie number of perturbing bodies is 4 then the 
effects of bodies 5 (Jupiter) and 6 (Saturn) will not be included 


See also ARC Set - OPTION CARDS - BODIES 

1 . 1 - 4.1 \ ‘ 

I « 
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COMMON’ SET 


GROUP 2 - OPTION’ CARDS 
BPART 
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NAME columns format description 

BPART 1-5 AS Specif ics^ that 

a) The B-Matrices create 

* * .in this run v/ill be for” 

: * * - in segments if necessarv 

w 4 

b) That no outer iterati 
adjustment of common par 

/ ... . . ters will be performed i 

: . to of the arcs ci*eates a B- 

- *• Matrix. 


t 

i 

1 
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BPART (Cont.) 


NAME 


COLUMNS 
1-5 (cont.) 


FORMAT 






DESCRIPTION 

c). Solution for the arc 
parameters will be performed 
on the last iteration only 

■ » 

•if sufficient core is avail - 

1 able. ♦ 


7*8 


12 


Number of Segments to be utili: 
in formation of the B-Matrix. 
This is only an approximate 
number and may be modified by 
GEODYN. If values greater thar 
zero are placed in columns 11-2 
.26-40, or 41-55 then columns 7 - 
will be ignored. 


9-10 


12 


Number of measurements that v;ii 
be utilized in each operation 
on the normal equations. The 
time consumed for B-Matrix form 
tion by segments will be invers 
ly proportional to this value. 
The core required will be linen 
related to this value. Default 
20 . 


11-25 D15.8 

: 


Total amount of core allocated 
for the GEODYN execution in 
units of 1024 bytes. 




I.l-S.l 



••• • V ' BPART (Cont.) '• •• .. 1 

m P * # * 

' ♦ ‘ * * ' V 

NAME COLUMNS* FORMAT * DESCRIPTION 

26-40 D15.8 . Total amount of core allocated 

V v .. ’• ; '' . >. for the GEODYN execution less 

. * .*. • - that amount required for all 

'*■: • '' •*' I/O buffers. Units of 1024 byt 

• . V # ; . Vv.* J If this value is greater than 

/ m - 9 ' \ ’. zero, columns 7-8 and 11-25 wil 

• . *■ ... be ignored. 

•/ 41-55 D15.S Total amount of core to be 

i c * 

■/ utilized for B-Matrix formation 

/ i • • 

Units of 1024 bytes. If this 

* . • Value is greater than zero, 

columns 7-8, 11-25 and 20-40 v;i: 
be ignored. 


56-70 D15.8 Minimum E-Matrix Segment size 

in Double Words . The blocksize 
’ ' ■ . specified for unit 10 must be S 

. times this- . value. E-Matrix se**- 

REPRODUCIBILITY OF THE * mont sizes will be integral 

ORIGINAL PAGE IS POOR* ' ’ * multiples of this number. 

• g Minimum value for this number 

■ * ■ . is 226. Default value for th 

: number is 220. 


IF CARD OMITTED: 


Beatrices will not be formed in segments. 



!• 1 - 5. 2 


r 


H' 


Increases in running time due to segmented formation cf 
the B-Matrix are related to the parameters specified cn 
this card as follows. *•••-* 

o directly proportional to the number of segments 
(which is proportional to the amount of core), 


o inversely proportional to the blocksize for unit 
• : 16 , * * ; . 


o inversely proportional to the value in columns 
* 9 - 10 '. 




1 . 1 - 5.3 


COMMON SET 

GROUP 2 - OPTION CARDS 
CENTER 


August 11, 1973 
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* 



Name 


Columns 


Format Description 


CENTER 1-6 


A6 


Specifies change of central 
body in the forcing function. 


7 II 


11-25 . D15.8 


Central body indicator 
** 1 Moon 
« 2 Sun . 

** 3 Venus 
*= 4 Mars 
e 5 Jupiter 
« 6 Saturn 

GM C - Universal gravitational 
constant times mass of the 
central body. 


REPRODUCIBILITY OF Thi? 
ORIGINAL PAGE IS POOR 


* 

1.1-5. 4 


r 



Name 


Columns 


Format 


26-40 

D1S.8 

41-55 

D15. 8 

56-70 

D15.8 


Description , v 

f A r " Mean equatorial radius 
of the central body. 

f c “ Flattening of the central 
body. Used only for surface 
density positioning. 

Value greater than zero indi- 
cates trajectory output will 
be in mean of 1950 coordinates. 


IF CARD OMITTED: Central body will be Earth. 


Default values for Central Body Only. 


Potential Expansion 
Max Max 


Body 

GM r 


is. 

Deg- 

Ord . 

MOON 

4. 902778D+12 

1. 73809D+06 

. 378D-03 

4 

4 

SUN 

1 . 32712 5D+ 2 0 

6. 96000D+ 08 

. 500D-04 

4 

4 

VENUS 

3. 248602D+14 

6. 16881D+ 06 

o 

• 

o 

4 

4 

MARS 

4. 282843D+13 

3. 380422D+06 

. 1 OSD- 01 

4 

4 

JUPITER 

1. 267077D+17 

7. -13715540+07 

. 667D- 01 

4 

4 

SATURN 

3. 7926535D+16 

6. 0401128D+07 

. 105D- 00 

4 

4 

NOTES : 

This card should not be used 

if central 

body is 

Earth 


When this card is utilized the following inputs are 
referred to the central body: • • 

• Columns 7-8, 9-10, 56-70, 71-80 on the EARTH card, 

• RECOEF card information, 

• SURF and SURDEN card information, 



• Variance/covariancc matrix of elements as sped-, 
fied on YARCOV card, and 

• Epoch elements unless otherwise indicated on 
element cards. 

When the central body is not Earth, the drag and Earth 
tidal perturbations are automatically shut off. 

See. also Section 2.0, JCL for Unit 1 when CENTER card 
is used. 
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NAME 

COLUMNS 

FORMAT 

DESCRIPTION 

CORREL 

1-6 

A6 

""Specifies the__ correlation 
coefficients between the 
coordinates of a station to 
be adjusted. 

0 

7 

11 

«=1 - correlation coefficients 
input on earth-fixed X,Y,Z. 

• 



*= blank - correlation co- 
efficients input on latitude, 
longitude, and height. 


11-25 

D15.8 

.Station number. 


26-40 

D15.8 

Correlation coefficient 
between X and Y or latitude 
and longitude. 


41-55 

D15.8 

Correlation coefficient be- 
tween X and 7. or latitude 
and height. 


1 . 1-6 


M 



f 

- CORREL (Cont.) 

* 

r « 

NAME COLUMNS FORMAT DESCRIPTION 

D15.8 Correlation coefficient 

between Y and Z or longi- 
tude and height . 

The coordinates of* the station are assumed 
to be uncorrelated. 

NOTES: The maximum number of master and unconstrained 

adjusted stations* (see STAEST card in this Section) 
is 381; one CORREL card may appear for 
each such adjusted station. If no CORREL card 
appears for such an adjusted station, uncertainties 
in the a priori station position will be assumed 
uncorrelated. Constrained adjusted stations which 
are not master stations are assumed to have the 
same correlations as the master station. 

If a CORREL card appears for a station there must 
also be a STAEST card present for that station. 

If the station is a constrained station it must 
be a master station. 


* The terms constrained, unconstrained and master stations 
refer only to adjusted stations (i.e., those for which 
STAEST cards are present) . A group of constrained stations 
are stations which are adjusted such that the same station 
position corrections are to be applied to each station in 
the group. One station in each group is arbitrarily desig- 
nated as a master station and all other stations in the group 
must refer to that master station. An unconstrained station 
is an adjusted station which is not a member of any group 
of constrained stations. 


56-70 


IF CARD OMITTED: 


1,1-7 
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DPOLE 
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NAME 


DPOLE 


COLUMNS 


1 -5 


FORMAT DESCRIPTION 

AS Modifies the coefficients of 

the diurnal nutations of the 
pole* . 


7 II Specifies the disturbing body 

=1 moon 
= 2 sun 

11-25 D15.8 Multiplier of the mean anomaly 

of the moon (£,) . 


REPRODUCIBILITY of tee 

oXal page is POOR 


* 

Note: The DPOLE option card may be used to expand or reduce 

the existing tables in BLOCK DATA. 


1. 1-7.1 



DPOLE (Cont.) 


NAME COLUMNS 

FORMAT 

DESCRIPTION 

26-40 

D15.8 

Multiplier of the mean anomaly 
of the sun (2, 1 ) . 

41-55 

D15 . 8 

Multiplier of the mean angular 
distance of the moon from its 
ascending node (F) . 

56-70 

D15.8 

Multiplier of the mean elonga- 
tion of the moon from the 
sun (D) o 

71-80 

D10. 8 ■ 

Multiplier of the longitude 
of the mean ascending node of 
the moon's orbit (ft) . 

SECOND CARD 



COLUMNS 

FORMAT 

DESCRIPTION 

1-10 


Must-be left blank. 

r 

11-25 

D1 5 « 3 

Coefficient for the rotation 
axis. j 

26-40 

D15.3 

Coefficient for angular 

momentum. j 

! 

t 

Coefficient of Euler angle 
perturbations in longitude. 

r 

41-55 

D15.3 


9 


1.1-7. 2 


DPOLE (Cont.) 


The second card must be present if the DPOLE option 
is used. If the second card is left blank, then there will 
be a zero contribution for the fundamental argument multipl 
as specified on the DPOLE card. 

Any number of DPOLE cards may be used. 
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CROUP 2 - OPTION CARDS 
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NAME 

EARTH 


COLUMNS 

1“5 


. FORMAT 

AS 


DESCRIF-nON" 

Modifies the geopotential 
and/or introduces ne:; earth 
constants. 


7 -8 


12 


REPRODUCIBILITY of 1<Hi? 
ORIGINAL PAGE Is /qOR 


9-10 


12 


Maximum degree coefficient to 
be used in the geopotentiai 
model plus 1. 

Default - 23 # Minimum = 4 

Maximum order coefficient to 
be used in the geopotentiai 


i 

i 


iia 


- fm 


AUjjUac II, iy/o 


EARTH (Cont.) * 

# ' 

NAME COLUMNS FORMAT DESCRIPTION 

9-10 (Cont.) model plus 1. 

Default = 23, Minimum - 2 

a 

11-25 D1S.8 GM - The gravitational constant 

times the mass of the earth. 

‘26-40 D1S.8 .a e - Semi-major axis of the 

earth. 

41-55 D15.8 Inverse of the earth flatten- 

ing. (1/f) 

56-70 D15.8 Odd value greater than zero - 

All geopotential coefficients 
set to zero except those which 
appear on RECOEF cards which 
follow this card. 

^ Value greater than 1.0 - 

Requests that geopotential 
coefficients printed in run 

heading be in normalized form 
rather than the default de- 
normalizcd form. 

‘71-80 D18.8 Fractional uncertainty in GM of the 

central body (see CENTER card) . 

For value other than zero, GM of 
the central body will be adjusted,. 

The default values of the earth constants as used by 
GEODYN are: 

GM = 3. 986013D + 14 m 3 /sec 2 

“ 6378155. m. 
e 

l/f= 298.255 

1.1-9 


M 


EARTH (CONT.) 

t 

j 

* 

IF CARD OMITTED: The stored SAO 1969 geopotential field 

containing terms up to degree and order 22 will be used 
with the default values of the Earth constants unless f 
modified by RECOEF cards in COMMON SET or EARTH and 
RECOEF cards in 'ARC SET. 

The set-up for two' body Earth computation is 
explained in ARC SET, OPTION CARD — EARTH. 

See ARC SET, OPTION CARD -- EARTH for individual 
arc specification. 
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. REPRODUCltiiLl'iY OJ.’ the 

NAME COLUMNS FORMAT ... DESCRIPTION 0RIGINAL PAGE ® POOR 

EDIT 1-4 A4 Used to change initial outer 

iteration RMS. 

7 ' Value greater than zero indicates 

that automatic editing of arcs 
will be performed. 


26-40 D15.8 Initial outer Iteration RMS. 

Used for outer iteration con- 
vergence test. Default = 1000. 


4 55 D15.8 First iteration edit criterion 

for electronic bias passes. 
Default - use same editing 
as other data uses. 


IF CARD OMITTED: Default RMS used, and no automatic editing will be 

performed. 


1.1-11 


r 
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COMMON SET 

GROUP 2 - OPTION CARDS 

• t 

ELCUT 



NAME COLUMNS FORMAT 

ELCUT 1-5 AS 




11-25 


D15.S 


DESCRIPTION 

Specifies the elevation 
angle below which measurements 
will be deleted from the 
solution. 

Elevation cut-off angle 
in degrees. 

Default = 0° 


IF CARD OMITTED: All measurement with elevation angles less 

than 0° will be edited from the solution 
unless modified in ARC SET by ELCUT card. 

Elevation cut-off angle may be negative. 

See also ARC SET, OPTION CARD - ELCUT. 
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COMMON SET 

GROUP 2 - OPTION CARDS 
.FLUX 
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NAME 


COLUMNS 


FORMAT 


DESCRIPTION 


FLUX 


1-4 


A4 


Adds and modifies solar 

and/or magnetic flux values 
in built-in tables. 


11-25 


D15.8 


Date of flux values (in 
YYMMDD ). 


26-40 


D15 , 8 


Value of solar flux. If 
blank, value in table will 
not be altered, f. 


10.7 


1.1-13 


r* 



. I ■ * 


NAME 



IF CARD 
Note: 


I 


FLUX (Cont.) 


COLUMNS FORMAT DESCRIPTION 

41-55 D15.8 Value of magnetic flux. If 

. ’ . . •/* blank, value in table will not 

. ' . . -be altered. Input value is th 

v' : •. - '- daily sum of the eight three- 

> .• \ • . hourly values kp. 


OMITTED: Built-in tables are used. 


Flux values are daily values.. All zero values 
between the end of the built-in tables and the 
data added via FLUX cards will be filled by 
straight line interpolation between known values. 
For dates beyond the table entries and data added 
by flux cards, the most recent values arc used. 
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1973 
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NAME 


COLUMNS FORMAT DESCRIPTION 


J65 


1-3 


A3 


Changes the default atmospheric density 
model for all arcs to the Jacchia 1965a 
Static Density Model updated to 1968 
(herein designated J65) . 


IF CARD OMITTED: Default will be J65 unless elsewhere modified to 
J650 or J71. See also Epoch Card and Option Cards J650 and J71. 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 
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J650 1-4 A4 Changes default atmospheric density model to 

Old J65 (herein designated J650) . This 
model differs from J65 in several aspects [1]. 
It should be noted that these differences con 
stitute errors in the J650 model. The J650 
. ■ - model is available .in GEODYN only because thi 

model has previously been the only model used 
by GEODYN. 


IF CARD OMITTED: Default will be J65 unless elsewhere modified to J650 

or J71. See also Epoch Card and Option Cards J65 and J71. 

[lJ"Notes on the Application of 10.7 -.a Solar Plux Data in the Jacchia Sta 
Diffusion Models of the Upper Atmosphere and Related Topics," Kolf 
Research and Development Corporation, July 1973. 
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* NAME COLUMNS FORMAT 


DESCRIPTION 


J71 1-3 A3 


Changes default atmospheric density model 
for all arcs to the Jacchia 1971 Static 
Density Model (herein designated J71) . 


IF CARD OMITTED: Default will be J65 unless elsewhere modified to 
J650 or J71. See also Epoch Card and Option Cards J65 and J650. 


REPRODUCIBILITY OF TEE 
ORIGINAL PAGE IS POOR 
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NAME 


LOVE 2 


COLUMNS 


FORMAT DESCRIPTION 

A5 Requests that solid earth 

tide coefficients of the second 
and third kind be modified. 


0, this card pertains to 
both lunar and solar 
effects. 


original 


1, this card pertains to only 
lunar effects. 

2, this card pertains to only 
solar effects. 
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LOVE 2 (Cont.) 


NAME 

COLUMNS 

FORMAT 

DESCRIPTION 


11-25 

D15.8 

Love number of the second 
kind to account for radial 
displacement. h 2 


26-40 

D15 ♦ 8 

Love number of the third kind 
to account for horizontal 
shearing. Z^ 

IF CARD 

OMITTED: 




Default values are as follows: 

l 2 = 0.075 


THE 


ORIGINAL page IS v 


h, = 0.600 
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NAME COLUMNS FORMAT DESCRIPTION 


POLE 1-4 A4 


Used to modify and/or request the adjustment 
of the true pole of the Earth. 


7 II * =0 No adjustment to be made. 

>0 Adjustment requested. 

*=2 Solution constrained to adjustments on 

the great circle passing through the pole a 
defined by the longitudes specified in 
columns 41-70 of the second card in this 
group. 

11-25 D15.8 X of the pole in seconds of arc. 

26-40 D15.8 Y of the pole in seconds of arc. 


reproducibility of the 
1.1-14.4 ORIGINAL PAGE IS POOR 
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POLE (Cont.) 

t 

COLUMNS FORMAT DESCRIPTION 

« 

41-55 D15.8 X of the pole in seconds of arc per day. 

" 4 

# 

56-70 D15.8 Y of the pole in seconds of arc per Gay; 

SECOND, CARD 

1-6 A6 Blanks must be present. 


11-25 D15.8 , Start time for application of' polar coordin- 

ates in form YYMMDD . DDD 

26-40 D15.8 Stop time for application of polar coordinate 

in form YYMMDD. DDD 

41-55 .D15.8 A priori correlation coefficient between X an 

Y or (see First Card column 7) longitude of 
constraint for X and Y. 


56-70 D15.8 A priori correlation coefficient between X an- 

* 

Y or (see First Card column 7) longitude of 
constraint for X and Y. 


THIRD CARD (read only if adjustment in requested ; see First Card colum 


1-6 


A6 


Blanks must be present. All other fields on 
this card must contain values greater than ze 


11-25 D15.8 A priori uncertainty in X. Same units as X. 


26-40 D15.8 A priori uncertainty in Y. Same units as Y. 



POLE (Cont.) 


August 11, 197 J 


COLUMN'S 

FORMAT 

DESCRIPTION 





41-55 

D15.8 

A priori uncertainty in 

« 

X. 

Same 

units 

• 

as X . 

56-70 

D15.8 

A priori uncertainty in 

c 

Y. 

Same 

units 

as Y. 

IF CARD OMITTED: 

Default 

values of the pole will 

be 

used . 




Maximum number of POLE cards 337. 

Coordinates of the POLE valid for midpoint of time span. 
Overlapping time spans not allowed. 


c 
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NAME COLUMNS FORMAT 

RECOEF 1-6 A-6 

or 

COEF 

7-8 12 

9-10 12 


DESCRIPTION 

Modifies and/or requests the 
estimation of coefficients 
in the geopotential model. 

Degree of C and S coefficients, 
(n index) 

Order" cTfjC and S coefficients, 
(m index) 

REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 


1.1-15 



NAME 


COLUMNS FORMAT 

11-25 015 .8 

4 

26-40 015.8 

41-55 . 015.8 


> ; ' 
DESCRIPTION 

A priori (or starting) value 
of C coefficient. 

A priori (or starting) value 
of S coefficient. 

Standard deviation of C coefficient. 


56-70 


• Coefficient will not be adjusted 
if value is zero. 

• Value is fractional uncertainty if 
;>l«0D-5 ( </TjOD- 5 for J2) and 
standard deviation equals fractional 
uncertainty times nominal* value of 

a priori coefficient. The fractional 
uncertainty is the standard deviation 
- divided by the a priori value. 

Value is standard deviation if 
' <1 .OD-5 (/l.OD-5 for J2). 

Standard deviation of S coefficient 

Coefficient will not be adjusted if 
value is zero. 

Value is fractional uncertainty if 
>1.0D-5 and standard deviation equals 
fractional uncertainty times nominal* 
value of a priori coefficient. 

Value is standard deviation if 

<1.0D-5. 
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April 8, 1974 


RECOEP (Cont.) 

NAME COLUMNS FORMAT DESCRIPTION 

71-80 DIO. 4 =0 Values on this card will be 

interpreted as being unnormalized. 

>0 Entries in columns 11-40 of 
this card will be ignored. A priori 
values of coefficients are obtained 
from current geopotential model 
including all modification prior to 
this card. Input standard deviations 
are un-normalized. 

<0 Values on this card will be 
interpreted as being normalized. 

IF CARD OMITTED: The geopotential used will be. the GEM 1 

Model unless modified by EARTH card or modified in ARC SET 
by option cards, EARTH and REC0EF. 

See also ARC SET, OPTION CARD - RECOEF. 

The geopotential coefficient denormalization equation 
used is : 


c m = 

c m * 

D m 

n 

n 

n 

S* - 

* 

D m 

n 

n 

n 
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where 




** yf (4n 


2 ) 


(n • roJT 

(n + mj I 


for m f o 


*= V 2n + 1 for 


m * o 


and 


CjJ* S n are un-normalized coefficients 


are normalized coefficients 
* Nominal values of geopotential coefficients 


are 


1*0 - 


(10' S ) D® 


n 


A Priori Value of Adjusted Coefficients 




is desired that the starting values of geopotential 
coefficients be different than the ' a priori values for those 
coefficients then the RECOEF card should be used as follows: 

• RECOEF cards requesting the desired coefficient 
adjustments should be included in the setup 
deck. The a priori values will be indicated on 
these cards. 

• For each coefficient pair for which the starting 
value is to be different than the a priori value 
an additional RECOEF card should be in the setup 
deck. This RECOEF card will indicate the start- 
ing value and come later in the deck than the 
corresponding RECOEF card requesting adjustment. 
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ORIGINAL PAGE IS 


NAME COLUMNS FORMAT 


DESCRIPTION 


RESID 


1-S 


AS 


Controls residual printout 
for all arcs. 


II 


0 * Indicates no residuals are 
requested for any arc. 


1 - Indicates that residuals 
are requested on the first 
inner iteration of the 
first outer iteration for 
all arcs. 


1.1-19 


RESID (cont.) 


NAME COLUMNS 

7 (Cont.) 


IF CARD OMITTED: 


See also ARC SET 


FORMAT DESCRIPTION cont. 

2 - Indicates that residuals 
are requested on the last 
inner iteration of the last 
outer iteration for all 
. arcs. 

- Indicates that residuals 
are requested on the first 
inner iteration of the first 
outer iteration and the last 
inner iteration of the last 
outer iteration for all arcs. 

4 - Indicates residuals are 

requested on all iterations 
for all arcs. 


Residuals are printed out for the first 
and last iterations of each arc unless 
modified in each arc by the RESID card 
in ARC SET. 


, OPTION CARD - 


DpQTT) 

REPRODUCIBILITY of the 
ORIGINAL PAGE IS POOR 
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COMMON SET 

GROUP 2 - OPTION CARDS 
RSTART 



NAME 

RSTART 


COLUMNS FORMAT 


1-6 


A6 


DESCRIPTION 

Requests the restart of a 
previous run and/or requests 
that a restart tape be written 
in this run. 


7-8 


12 


Unit number of input restart 
tape* If zero or blank, job will 
start from the beginning. 


9-10 


12 


Unit number of output restart 
tape* If zero or blank, no 
output restart tape will be 
written. 


1 . 1-21 


r 



* RSTART (Cont.) 

NAME COLUMNS FORMAT DESCRIPTION 

11-25 D15 . 8 Number of arc in which termination 

occurred. If zero or blank,’ the 
job will start from the beginning. 

IF CARD OMITTED: The job will start from the beginning and 

no output restart tape will be written. 

NOTE: Set up changes may not be made such that the size 

or configuration of the normal matrix is altered 
. when restarting. 




G 
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COMMON SET 

GROUP 2 - OPTION CARDS 
SEQDAT # 



NAME 

COLUMNS 

FORMAT 

DESCRIPTION 

SEQDAT 

1-6 

A6 

Assumes arcs are in 
chronological order, 


IF CARD OMITTED: Makes no assumption as to order of arcs. 

If SEQDAT is included in set up, arcs absolutely MUST 
be in chronological order. 

This option is included to reduce running time by not 
rewinding the tape between arcs unless unit number changes. One 
data point will be lost between arcs when this option is used. 


REPRODUCIBILITY OF THE 
. . ORIGINAL PAGE IS POOR 
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* COMMON SET 
GROUP 2 - OPTION CARDS 
STAEST 
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1 ;; ;i *1 n u ;i n *t it 


NAME 


COLUMNS 


FORMAT 


DESCRIPTION 


STAEST 


1-6 


A6 Request estimation of a 
station position. 


II =1 - standard deviations 
input on station earth- 
fixed X,Y,2. 

*=b-lank - standard deviations 
input on station latitude, 
longitude and height. 


11-25 

26-40 


D15.8 Station number. 

D15.8 Standard deviation cf X in 

meters or latitude in seconds 

of arc. 


1.1-24 



; * March 1, 1973 

• # . *• • • * 

• • • * • 

• STAEST '(Cont.) 

NAME COLUMNS FORMAT DESCRIPTION 

41-55 D15.8 Standard deviation of Y in 

meters or longitude in seconds 
. * /% . ; * . of arc. . 

56-70 D15.8 Standard deviation of Z in # 

; • ; ..meters or height in meters. 

71-80 F10.4 Master station to which the adjust- 

• ** 

ment of this station is to be 
.. . *y • constrained. If station is 

•; * unconstrained or a master station 

.. leave blank.* , 

♦ 

IF CARD OMITTED: Station positions will be held fixed. 

NOTES: A STAEST card mu s t be present for each station to be 

adjusted. If two or more stations are constrained 

'V 

one of these stations mu s t be specified as the master 
station by leaving columns 71-80 blank. 

A maximum of 381 unconstrained and master stations may 
be adjusted. The total number of adjusted stations 
including constrained stations may hot exceed 3S1. 

.STAEST card’s may not be present for stations that 
are not in the tracking complement. 

*See CORREL card in this section for definitions of 
constrained, unconstrained , and master . stations . 

See also STAPOS card. REPRODUCIBILITY OF THE 

. ORIGINAL PAGE IS POOR 
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COMMON SfiT 

GROUP 2 - OPTION CARDS 
• STEP 


4 


poooo o«a •: Ci : ? r -.i z* o 
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REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 
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■NAME. 

STEP 


. CO LUMNS. 

1-4 


FORMAT DESCRIPTION 

A4 Modifies variable step 

integrator error bounds 
and/or maximum and minimum 
step sizes. 


26-40 


D15.8 Error tolerance for step 

size decrease. Default is 
0.25D-4 m. 


41-55 


D15.8 Error tolerance for step 

size increase. Default is 
0* 2SD-1Q m. 


1.1-26 


r 


STHP (Cont.) 


March 1,1' 


.NAME 


COLUMNS 

56-70 


. 


FORMAT ' 

D15.8 


DESCRIPTION 


Maximum step size in variable 
step mode. Default is 400 
seconds. , 


71-80 


DIO. 5 Minimum step size in variable 

Step mode. Default is 12,5 
* 

' seconds . * ‘ 


IF .CARD OMITTED: Default values giVen above arc used 

unless modified by STEP card in ARC SET. 


NOTE : 


Only those bounds for which non-zero values appear 
on the STEP card will be modified. 


See also ARC SET, OPTION CARD - STEP. 




% 






i 
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• • COMMON SET 

GROUP 2 - OPTION CARDS ’ • 


stinnF.M 



NAME COLUMNS ' FORMAT DESCRIPTION 

■ ' - * 1 — » * ■ ■ — ■ — i ■ . . — .... , * 

SURDEN 1-6 . A6 Specifics number of intervals ir 


* latitude and longitude used for 
. surface area integration of sur- 

face densities. The trivial car 

# 

* . degenerates to a Mas con ♦ 

. 7-8 12 Number of latitude intervals. 

Default =1. 

.• 12 Number of longitude intervals. 

Default « 1. 

m 

• ' REPRODUCIBILITY OP twp 

IF CARD OMITTED: Default values used. ORDINAL PAGE Is Poor 

T 


1 1 - 0 



GROUP 2 - OPTION CARDS 
SURF 

9 
















SURF (Cont.) 


NAME COLUMNS FORMAT . . DESCRIPTION 

41-55 D15.8 Value of geocentric latitude at 

center of master block* in 
degrees. 

• ; 

56-70 D15.8 Value of east longitude at center 

of master block* in degrees* 

A second card must follow the SURF card in the input stream. 

The format of this card is: 

NAME COLUMNS FORMAT DESCRIPTION 

1-6 A6 These columns must be blank. 

/ • . 

* 

7-8 12 Number of latitude subdivisions 

in master block.* 

Default = 1 

9-10 12 Number of longitude subdivisions 

in master block.* 

Default =1 


11-25 D15.8 Size of each latitude subdivision 

in degrees . 

26-40 D15.8 Size of each longitude subdivision 

in degrees. 

WARNING: Failure to provide the second card or to specify 

sizes of latitude and longitude subdivisions will 
result in program termination. 



t 

SURF (Cont.) 


* March 1 , 


NOTES; Surface densities are an expression of deviations frer. 
the earth's gravitational potential as defined by 
♦ the spherical harmonic expansion of the geopotential 
* * coefficients. ... . f 

Surface densities are applied in GEODYN as uniform 
. . payers of mass covering areas of the Earth's surface. 
These areas are defined as the area of blocks bounded 

K 

by lines of latitude and longitude. 


*A master block is an area bounded by lines of latitude 
and longitude which may be subdivided into blocks of 
equal size (in terms of boundary lengths in degrees 
latitude and longitude) each of which has the same 
a priori value of surface density. 

A master ' block may be so- constructed as to contain 

only one block. This is the program^ default if the 

number^ of subdivisions is zero. 

% 

m 

For convenience, master blocks are used in GEODYN to 
input surface densities. 

A maximum of 674 master blocks may be specified. 


p * 
p 
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COMMON SET 

GROUP 2 - OPTION CARDS 
TIDES 
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* . * reproducibility of the 

name COLUMN'S FORiMAT DHSCRIPTION ORIGINAL PAGE IS POOR 

TIDES 1-5 AS Modifies earth tide perturbations 


TIDES 1-5 
7 


«0 f This card pertains to both 
Lunar and Solar tidal effects. 


^1, This card pertains only to 
Lunar tidal effects. 

*=2, This card pertains only to 
Solar tidal effects. 


1.1-31 


r 


I A 


* 


•NAME COLUMNS 

8 



FORMAT 



. ... r 

* j . 

DESCRIPTION . 

Indicates which tidal coefficient 
.is to be modified. 

■Value ' Coefficient 

* I Ml ft ' '**' 

‘ ;. r 1 K2 

V 2 v K3 

• : *• 3 ' K2 - Phase 

A value of zero in this column 
‘will set all coefficients to zero. 


» 


11-25 



D15.3 


Value of tidal parameter being 
modified . 

Defaults : 

. Coeff i cicn t Value 


K2 0.29 

K3 0.0 


K2-Phase 


2.5 Degrees 


26-40 D15.3 Standard deviation of tidal co- 

efficient. If non-zero value is 
specified, t i da 1^ parameter will 
be adjusted. Default values arc 
all zero , 


If card omitted : Default values given above arc used. 


REPRODUCIBILITY of the 
ORIGINAL PAGE IS POOR 



1.1-32 
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'COMMON SET 
* 

GROUP 2 - OPTION CARDS 
TOLS 





NAME COLUMNS 

TOLS 1-4 


FORMAT 

A4 


DESCRIPTION 

Modifies integration tolerances. 


7 II Ind icates for which satellite 

tolerances are to be modified 
for all arcs. 

■ 0 both satellites 
*= 1 first satellite 
*= 2 'second satellite 


9-10 


12 Maximum degree and order of 


BEPBODUCIBIL1TY of the 

ORIGINAL PAGE IS POOR 


1.1-33 
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• TOLS (Cont.) 
• « 

NAME COLUMNS FORMAT 

9-10 (Cont.) 


11-25 D15.8 


Critical predictor -corrector 
tolerance. Used to determine 
number of corrector iterations 
necessary at each integration 
step* If after 2 corrections 
the. local error exceeds critical 
tolerance an- error message will 
be printed (execution will be 
terminated) . 

Default = 0.25D-10 m. 

Maximum degree (and order) of 
* geo'potential coefficients for 
which constraint equations will 
be included in the estimation 
of surface densities. Default - 1 

IF CARD OMITTED: No modifications will be made to 

. integration tolerances. 


26-40 D15.8 


TT TTV OF 

1S poor 

origr^ v 


41-55 D15.3 


DESCRIPTION 

geopotential effects to be 
included in the computation 
of the variational equations. 
Default— 4, Minimum - 3. 

Desired reasonable predictor- 
corrector tolerance to be 
achieved by changing the inte- 
grator step size while in variable 
step integration mode. Used 
only in computation of new 
step size during step size 
variation. 

Default = 0.25D-7 m. 


. " . March 1, 15’ 

1.1.3 COMMON SET 

GROUP 3 - STATION COORDINATE CARDS 

i • 

l , 

* 1 

* •’ . '• V * STAPOS 


CARD TYPE 


FORMAT COLUMNS DESCRIPTION 


Starting Card A6 


1-6 The word STAPOS 

beginning in 
. . : column one. 


II 



II 


7 >0 indicates station 

ad j ustment sigmas on 
this card apply to 
• Earth centered fixed 
••• ® Cartesian coordinates. 

$ >0 indicates that no 

station coordinates 
• are to be used from 

built- in coaplcir.cn : . 



D1S.8 


D15. 8 


D15.8 


11-2S 


26-40 


41-S5 

% 


Sigma in seconds of 
arc for station latitudes 
or-- meters for station 

ECFX. 

Sigma in seconds of 
arc for station longi- 
tudes or meters for° 
station ECFY. 

Sigma in meters for 
station height or ECPZ. 


NOTE: 


***•«• in «>«-» ii -2$, 

"fine the si E .; s input H /“I '•«. 

Which STAEST cards are present win t -i tl ° n ? f ° r 
Specified on the STAEST card*? \ bf ^ lC sl S mas 

» srS?“c„s f “;j n xx: 

REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 

• m 9 


r 


GROUP 3 - STATION COORDINATE CARDS 

(CONTO 


CARD TYPE 

FORMAT 

COLUMNS 

DESCRIPTION 

Station 

A6 

1-6 . 

Station 

Name 

Position Card 
for input 

14 

7-10 

Station 

Number 

• 

A1 

11 

Sign of 

Latitude 

'J/s.s j 

C 12 

12-13 

Degrees 

Station v 

} 12 

14-15 

Minutes 

Geodetic > 

[^FIO.8 

16-25 

Seconds 

Latitude ^ 


* 


Decimal 

may appear 


* 

- 

only in 

col* 18* 

'Jls-.s 

f!3 

26-28 

Degrees 

Station v 

7 12 

29-30 

Minutes 

East ? 


v^FlO.8 

31-40 

Seconds 

Longitude^ 




Decimal 

may appear 




only in 

col. 33. 

- 

D1S.8 

41-55 

Station 

Height ini 

- 

. 


Meters 

j 


A1 

56 

7th character of 


•• 


station 

name . 

. Station 
Position Card 

A6 * 

1-6 

Station 

Name 

for Barth-Fixed 

1 4 

7-10 

Station 

Number 

Cartesian Coor- 
dinates 

D15.8* 

11-25 

x in meters 

D1 5*8 

26-40 

y in meters 


D15.8 

41-S5 

z in meters 

* 

A1 

56 

7th character of . 


station name 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR ' 
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NOTH: 


'* 4 

GROUP 3 - STATION COORDINATE CARDS 
(CONT.) 



GROUP 3 is necessary only if the tracking complement 
built into the program is to be changed or replaced. 

* f * . ‘ . ‘ * . 

The number of station position' cards must not exceed 
381 . - • ■ • ‘ 

These cards are Optional . 

% - • 

The station height is measured in meters above the 
reference ellipsoid. See EARTH card, COMMON SET, 
GROUP 2 - OPTION CARDS. 

- ■ • 

See Appendix D for built-in. station positions. 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 


1.1-36.1 
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1.1.4 COMMON SET 

GROUP 4 - TERMINATION CARD, 

♦ 

NAME COLUMNS FORMAT DESCRIPTION 


DATA 1-4 
7-8 

9 

10 


A6 

12 


II 


Designates end of COMMON SET. 

Convergence criterion* for 
outer iterations in percent. 
Default is 2 percent. 

Minimum number of outer 
iterations. Default is 1. 


II 


Maximum number of outer 
iterations. Default is 1. 


This card terminates the COMMON SET containing cards to 
define information common to all arcs and is Mandatory . 


*When 


RMS previous outer iteration - 
RMS current outer iteration 

RMS current outer iteration 


x 100 


< 


converge: 

criterioi 


convergence has been obtained. On the first outer iteration 
the values of the RMS for the previous outeT iteration is, by 
default, 1000. This can be changed by using columns 26-40 of 
COMMON SET, OPTION CARD EDIT. 

Inner and outer iterations are explained in the Introduction 
of this manual. 


1.1-37 


r 



1.2 ARC SET • 


• • 


'ARC SET contains that information which is necessary 
to define each individual arc. 



Included in ARC SET are four groups: 



Al'c description cards which are .three title 
cards, the Reference Time card, the Epoch card, 
and the Epoch Orbital Element cards. (Mandator 



Option cards which may be used to exercise the 


t 

i 


GEODYN individual arc capabilities. (Optional) 


. 3* Termination card for Group 2 - Option Cards. 

(Mandatory) 

m “ * * * -9 • 

A* Selection/ deletion cards to exercise GEODYN*' s 
data selection and deletion capabilities. 
(Optional) 


The. following pages describe these groups more 


thoroughly . 







• ‘ REPRODUCIBILITY OF THE 
• ORIGINAL PAGE -IS' POOR 
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1*2.1 ARC SET 

GROUP 1 - ARC DESCRIPTION 


CARD NUMBER 

1A 

IB 

1C 


FIELD FORMAT DESCRIPTION 


10A8 

10A8 

10A8 


Three title cards for each arc. 
Columns 1-80 of each card are- 
printed as identification on 
all output from the arc. 


2 


3 

/ 

4A 


4B 


Year, month, day of the 
reference time (YYMMDD) . 


See Section 
1.2. 1.3 
3D24.16 


Epoch card. 

Element card one (Position 
vector] . 


3D24.16 


Element card two (velocity 
vector) . 


These cards are Mandatory for each arc and are more fully 
explained on the following pages. 


If two satellites are in the arc cards 4A and 4B should 
be repeated specifying the second satellite epoch elements. 

(See EPOCH card.) 


REPRODUCIBILITY OR THE 
ORIGINAL PAGE IS POOR 
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• ARC SET 

GROUP 1 - ARC DESCRIPTION 
TITLE CARDS 


' “'I " *' 


5 DAY GEOS- 1! ARC — ADJUSTING SOLAR RADIATION 
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CARD 

COLUMNS 

FORMAT 

DESCRIPTION 

1A 

1-80 

10A8 

Identification to be printed 

IB 

1-80 

10A8 

on all output from arc. 

1C 

1-80 

10A8 

- 


Three title cards are Mandatory for each arc. 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 


1.2-3 


. 


ARC SET 
# 

CROUP 1 - ARC DESCRIPTION 
REFERENCE TIME CARD 


COLUMN'S 



•...FORMAT 





v< 


DESCRIPTION 


1*6 • ;‘I6 . Year, month, day of refer 

’ . •*« . ence date in YYMMDD. 


Note: This time defines the inertial coordinate 
/ system in which the program 'will integrate 

the orbit. The -reference time should precede 
or be the same as the epoch of the initial 
elements. 

- - The input epoch elements are assumed to be in 
the -inertial coordinate system defined by the 
• true equator and equinox of the reference time. 
The estimated (adjusted) epoch elements are 
also in the same coordinate system. 

It is critical to proper prorrnm operation that 
each column from 1-6 contain numerals between 
0 and 9 and that no blanhs be present in these 
columns. Columns 7-80 must remain entirely 
„* blanh. 

This card is Mandatory for each arc. 

• • • / REPROL U Uitflui’ Y {*y' 

* • ORIGINAL PAGE IS POOR 


COLUMNS 

. 1-6 

7-10 

4 

11-17 

18 

19 

20 

21-22 

23 


ARC SET 

GROUP 1 - ARC DESCRIPTION 
EPOCH CARD 


FORMAT 

16 

14 

F7.4 

II 


II 


II 


2PF2.0 


II 


\ 

DESCRIPTION 


Epoch date of the initial ele- 
ments of the arc in YYMMDD. 

Hours and minutes of epoch in 
HHMM. 


Seconds of epoch in SS.SSSS. 

Maximum number* of inner iterations 
on all outer iterations except the 
first for this arc. The default 
is 1. 

Maximum number* of inner itera- 
tions for this arc on the first 
outer iteration. The default is 1. 

Minimum number of inner iterations 
for this arc. The default is 1. 

Inner iteration convergence criteria 
for this arc. The default is 2 per- 
cent. 

Tape format flag. 

1 - indicates GEOS-1 Data Center format 
0 - indicates binary DODS data tape 

2 - indicates GEOS-2 Data Center format 

3 - indicates PCI2 format 

4 - indicates SIMULATED Data format 

5 - indicates GEOS-C Data format 

6 - indicates RAP Binary Data format 


COLUMNS 

24-29 

30-33 

34-40 

•41-46 

47-50 

51-57 

58 

59-64 

65 

66-72 

77-80 


EPOCH CARD (£ont.) 
DESCRIPTION 


Cutoff date for data in YYMMDD. 

Hours, minutes of cutoff time for 
data in HHMM. '• 

Seconds of cutoff time for data 
in SS.SSSS. 

Start date. for data in YYMMDD. 

Hours, minutes of start time for 
data in HHMM. 

Seconds of start time for data 
in SS.SSSS. 

II Number of satellites. Default - 1. 

17 . DODS Data- Base element set number. 

II DODS Data Base element retrieval 

flag (See Appendix A) . 

17 Satellite ID for element retrieval. 

A4 This field is used to alter the atmos 

pheric density model used for the arc 
Models available are: 


FORMAT 

16 

14 

F7.4 

16 

14 

F7.4 



August 11, 1975 


* MNEMONIC DESCRIPTION 

J65 Jacchia 1965a modified to 

1968 (see option card of same 
name). 

J650 Old Jacchia 1965a modified to 

1968 (see option card of same 
name) . 

J71 • Jacchia 1971 model (see option 

card if same name) . 

This card is Mandatory for every arc. 


1 . 2 - 6.1 

p 



EPOCH CARD (Cont.) 


NOTES: The Epoch and start and stop times for data must 

be specified in the UTC time system. The default 
for the data start time is Epoch* No data prior 
to the start time or later than the cutoff time 
. will be accepted. Orbit integration will be initiated 
from Epoch or the data start time, whichever is 
earlier. The data start time may be used during an 
orbit generator run to specify the earliest time for 
orbit integration. A cutoff time specification is 
mandatory only if the data is being read from the 
DODS data base. ( See Appendix A.) 

t . If the data start time is earlier than Epoch, integration 
proceeds backward from Epoch to the start. Epoch is then 
reset to the start time, and integration resumes forward. 

Inner and outer iterations are explained in the Intro- 
duction of this manual. 


*If column 18 contains a 1 and' column 19 is zero or 
• blank the defaults for the maximum number of inner 
iterations will be 10 for the first outer and 1 for 
subsequent outer iterations. 



* 

• ARC SET * 

GROUP 1 - ARC DESCRIPTION 
POSITION VECTOR 
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COLUMNS FORMAT 


DESCRIPTION 


1-24 D24.16 x Component ) 

25-48 D24.16 y Component ) in meters 

49-72 D24.16 2 Component j 

of satellite position at epoch 
time, in the inertial coordinate 
system of the reference time, or 
osculating Keplerian elements. 

a (in meters), e, and 
i (in degrees) . 

REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 


1 . 2-8 



POSITION VECTOR CCont.) 


COLUMNS 

73 


This card 


FORMAT DESCRIPTION 

II *0 or blank, orbit will be elliptic. 

orbit can be either elliptic or 
hyperbolic and the input elements 
are Cartesian. 

«2, orbit can be either elliptic or 
hyperbolic and the input elements 
are Kepler ian. 

is Mandatory for each arc. 





GRi 
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COLUMNS 

FORMAT 

1-24 

D24.16 

25-48 

D24.16 

49-72 

D24.16 


DESCRIPTION 1 


jt component ) 

^ component ; in meters/sec. 

£ component j 

of satellite velocity at epoch time 
in the inertial coordinate system 
of the reference time, or osculating 
Keplerian elements Q, w and M (all 
in degrees) . 


This card is Mandatory for each arc. See Position 
Vector Card - Col. 73. 


°^Zl c So?l Z™ 
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1.2.2 ARC SET 

GROUP 2 - OPTION CARDS 


; >UgU.>l X -L , X 
I 


% 


Each GROUP 2 card is identified by the name beginning 
in column 1 and read under an A6 format* The name is followed 
by 9 fields: 


Field Columns Field Format 


7, -8, 9, 10 

411 

11-2S 

D15.8 

- 26? 40 

D15.8 

41-SS 

D15.8 

56-70 

D1 5 * 8 

71-80 

DIO. 8 


The Option Cards are: 


BIAS 

BMATRX 

BODIES 

CULL 

DRAG 

EARTH 

EBIAS 

EDIT 

MAX L AT 


ELCUT 

SAT 

ORB1 

SIGMA 

ORBIT 

SIMDAT 

OUTPUT 

SOLRAD 

PREPRO 

STEP 

RECOEF (COEF) 

TOLS 

RESID 

VARCOV 

TERMEL 

REG I NT 


original page 


1 . 2-11 


Il 


ARC SET 

GROUP 2 - OPTION CARDS 
ALTMTR 



NAME COLUMNS FORMAT DESCRIPTION 

ALTMTR 1-6 A6 Invokes altimeter data process- 

ing options. 

7 .11 =1 Requests output of an alti- 

meter experiment file on the 
last iteration. Output is 
on unit 81. 


8 II =1 


REPRODUCIBILITY of the 
ORIGINAL PAGE IS POOR 


Requests application of 
altimeter attitude error 
correction. This correction 
contains a bias term which is 
always present. (applicable 
only to GEOS-C). 


1 . 2 - 11.1 



ALTMTR (Cont.) 


NAME 

COLUMNS 

FORjMAT 

DESCRIPTION 


9 

11 

e l -Indicates that an attitude 
tape for GEOS-C will be 
input on unit 80. 


IF CARD OMITTED: No such options will .be applied. 



. I 

ARC SET 

* GROUP 2 - OPTION CARDS 
.BIAS 
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NAME COLUMNS FORMAT 


DESCRIPTION 


BIAS 1-4 


A4 Requests adjustment of either 

measurement biases or station 
timing biases. 


7-10 


14 Station number for which the 
bias is requested. 


11-25 


D15.8 


REPRODUCIBILITY of the 

original page is poor 


Optional a priori estimate of 

the bias value in "natural" 

. » 

units: Length: meters 

Angles: radians 

Time: seconds 

Mathematical symbols 

b - measurement bias 
At - timing bias 


1 . 2-12 
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' BIAS (Cont.) 


NAME COLUMNS FORMAT 

26-40 D1S . 8 

41-46 16 

47-53 F7.0 


56-61 16 


62-68 F7.0 


71-80 DIO. 8 


DESCRIPTION 

Standard deviation of the 
a priori bias in the same 
units. 

Optional year, month, and 
day of the time at which 
the bias is to begin. 

Optional hour, minute, and 
second of the time at which 
the bias is to begin. A 
decimal point must appear in 
column 55 . 

♦ Optional year, month, and day 
of the time at which the bias 
is to end . 

Optional hour, minute, and 
second of the time at which 
the bias is to end. A decimal 
point 'must appear in Column 68 

Bias type - same as those 
listed under SIGMA 
cards - for measurement biases 
.or zero for timing biases. 


IF CARD OMITTED:. No measurement or station timing biases 

will be adjusted. 


1.2-13 



'BIAS .(CONT.) 


.Note: The normal procedure for bias adjustment is for 

independent biases to apply for each satellite 
pass. If this is the desired type of adjustment, 
columns 41-70 should be left blank. When start 
and end times for biases are explicitly indicated 
they will override the pass-by-pass assignment; 
however, if there are passes after the stop time, 
new biases will be assigned for the passes. 




ORIGINAL PAGe1sK)<>b E 
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ARC SET. 

GROUP 2 - OPTION CARDS 
BMATRX 
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NAME 


BMATRX 


COLUMNS FORMAT 


1-6 


A6 


REPRODUCIBILITY OF THE 
DESCR9Py¥Sf L PAGE IS POOR 

Indicates that a B-matrix 
tape is to be written on 
unit 71 for this arc. 


11-25 D15.8 B-matrix number to be written 

on the tape. This number must 
be greater than zero and less 
than 10000. 

26-40 D15.8 >0 Requests that the B-matrix 

.also be printed on unit 6. 


IF CARD OMITTED: No B-matrix tape will be written. 


1.2-15 


w 
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BMATRX (Cont . ) ' • 

* 

NOTES: If any arcs write a B-matrix tape then the maximum 

number of outer iterations will be one and there 
will be no outer iteration adjustment* .s 

Adjusted parameters that will be included on the 
B-matrix tape are the following: 

Geopotential coefficients adjusted as common 
parameters. (See COMMON SET, OPTION CARD - 
RECOEF) 

Master tracking stations and unconstrained 
adjusted tracking stations. (See COMMON SET, 
OPTION CARD - STAEST) 

* 

Drag and solar radiation coefficients, Cp, C^, 
and C R . (See COMMON SET, OPTION CARDS - 
DRAG and SOLRAD) ► 

Satellite elements. (See ARC SET, Run Descrip- 
tion - Position and Velocity Vector cards) 
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0 

ARC SET * 

■ GROUP 2 - OPTION CARDS 
BODIES 



* NAME COLUMNS 

BODIES 1-6 

7 


FORMAT 

A6 


II 


DESCRIPTION 


Modifies non-earth gravity 
perturbations . 

Indicates body for which 
gravity perturbation is to 
be modified 


1 

2 

3 

4 


- 'MOON 

- SUN 

- VENUS 

- MARS 


EEPRODUCIBILITY of TH£ 
ORIGINAL PAGE IS POOR 
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BODIES (Cont.) 

NAME COLUMNS FORMAT 


7 XCont.) 


8 II 


, 11-25 D15.8 


M 2 /M = 332945.5619264376 
M,/M = 0.81S00032293349SS 

^PKODUGBIUTY OF THE mX - 0.1074468525270073 

St PAGE IS POOR M s /M e = 317.8809303 

*V W e = 95.14905175 

IF CARD OMITTED: Lunar and solar gravitational perturbations 

will be applied using default mass ratios 
unless BODIES card is used in COMMON SET. 

To omit a particular body for gravitational purposes a 
BODIES card representing that body must be included with the 
corresponding mass ratio set to zero. 

*For example, if the number of perturbing bodies is 4 then the 
effects of bodies 5 (Jupiter) and 6 (Saturn) will not be 
included, etc. 

See also COMMON SET, OTTION CARD- -BODIES. 


DESCRIPTION 

5 - JUPITER 

6 - SATURN 

Number* of perturbing 
bodies ‘to be applied 

(maximum of 6) « 

Default = 2. 

This number may be set on 
any one of the BODIES cards 
used in the arc. 

Value of ratio of mass of 
body to mass of earth. If 
ratio is zero perturbing 
effect will be zero. 

Default mass ratios: 

Mj/M e - 0.01229997171006507 
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• ARC SET. 

CROUP 2 - OPTION CARDS 
CULL 
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,WA AVAf/MV.W.V.VA.^ J 


NAME 


COLUMNS 


FORMAT 


DESCRIPTION 


CULL 


1-4 


A4 Used to delete individual 
measurements from the 
solution by zero weighting. 


11-25 


D15.8 The number* of the first measurement 
in a series to be deleted. 


26-40 


D15.8 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 


The number* of the last measurement 
in a series to be deleted; may 
be left blank if there is 
only one measurement in 
the series. 


1.2-19 


CULL (Cont.) 


NAME COLUMNS FORMAT 


DESCRIPTION 


41-55 D15.8 The number* of the first 

measurement in second series 
■ to be deleted. 


56-70 D15.8 The number of the last measure 

ment in second series to be 
deleted; may be left blank if 
there is only one measurement 
in series two. 


IF CARD OMITTED: No measurements are deleted. 


NOTE: A maximum of 50 cull cards may be used per arc. 

The series may fall in any order, and two series 
per card are not necessary. » 


*The numbers which correspond to each measurement are 
printed at the far right of the residual printout. 
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NAME 

DELAY 

COLUMNS 

1-5 

FORMAT 

A5 


7 

11 



Z' 

8 

11 


DESCRIPTION 

Changes transponder delay 
curves . 


Transponder type - 

1 VHF GRARR 

2 S-Band GRARR Channel A 

3 S-Band GRARR Channel B 

4 S-Band GRARR Channel C 

5 USB/DSN 

6 ATSR 

7’ C-Band Coherent 
8 C-Band Non-coherent 

t - h 

=n where n is the n relay 
for satellite to satellite 
data as follows: 


1 . 2 - 20.1 




DELAY (Cont.) 


NAME 


COLUMNS FORMAT 

8 (Cont.) 


DESCRIPTION 

n Transponder for SST 

1 uplink delay at relay 
satellite 

2 delay at tracked satellite 

3 downlink delay at relay 
satellite 

This column must be left blank 
for data other than SST data. 


9 II The number of the tracked 

satellite as specified by the 
order of the elements following 
the epoch card (default = 1). 

=1 if 1st set of elements belong 
to the tracked satellite. 

=2 if the 2nd set of elements 
belong to the tracked 
satellite. 


10 II The number of the relay satellite 

as determined by the order of the 
elements following the epoch 
' card (required only for SST). 

=1 if 1st set of elements belong 
to the relay satellite. 

=2 if 2nd set of elements belong 
, to the relay satellite. 

11-25 D15.3 Transponder delay bias (in 

meters) (D^) 


1 . 2 - 20.2 



DELAY (Cont.) 


NAME 

COLUMNS 

FORMAT 

DESCRIPTION 

• 

26-40 

D1S.3 

Transponder delay coefficient 
(in seconds) of the linear 
variation with 


41-55 

D1S.3 

Transponder delay coefficient 
(in sec /m) of the quadratic 
variation with p . (D^) 


The total transponder delay in meters is computed as 
follows: 

DELAY = D x + D 2 p + Dp 2 
and the correction for range data is 

p true ~ p observed ~ 2 DELAY 


where: 

0 * 
p = range rate in m/sec 

Observed = observed range 
p true = correcte d observed range 


1,2-20.3 


April 8, 1974 


• 

DELAY 

(Cont . ) 



Example : 

- 



' 

given: 

• 

ft 

# 




(a SST situation) 

EPOCH 

CARD 




X 1 

V 

Z 1 
Z 1 | 

Elements for 
► tracked satellite, 


*1 

*i 

*1 ) 

| GEOS-C 


X 2 

• 

Y 2 

Z 2 ] 

Elements for 
■ relay satellite. 


X 2 

*2 

*2 ! 

ATS 


Since the order of the elements determines the satellite 
number*, the satellite number associated with the relay 
satell ite=2 , and the number associated with the tracked 
satellited . 

Suppose we wished to change the transponder delay curves 
for the signal received: 

at the relay satellite from the station 
at the tracked satellite from the relay satellite 
• at the relay satellite from the tracked satellite 
(this signal goes back to the station) 


The satellite number should not be confused with the satellite 
ID as specified on the SAT option card (see 1.2-46) in Columns 
41-55. 

o 


1 , 2 - 20.4 



XX w , X s / H 


DELAY (Cont.) 


The appropriate numbers 
follows : 

in columns 8, 

9, 10 would be as 

Column 8 

Column 9 

Column 10 

Changes Delay for the 
Signal Received: 

• 1 

1 

2 

at the relay sat from 
the station 

2 

1 

2 

. at the tracked sat from 
the relay sat 

3 

1 

2 

at the relay sat from the 
tracked satellite 



1.2-20.5 


DELAY (Cont.) 


In a SST situation the transponder delay can be 
modified : 

at the RELAY SAT when signal comes from station 

and/or 

at the TRACKED SAT when signal comes from 
RELAY SAT (2) 

and/or 

at the RELAY SAT when signal comes from the 
TRACKED SAT (3) and will be sent back to the 
station (4). 


r 


i 

• ARC SET 

GROUP 2 - OPTION CARDS 
DRAG 
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NAME COLUMNS FORMAT 


DESCRIPTION 


DRAG 1-4 


A4 Modifies drag application. 


II Satellite number (1 or 2) # * If 
less than 1 or greater than 2, 
applies to both satellites. 


11-25 


D15.8 


^PUODUCIBIUTY OP THE 
b po< * 


Drag coefficient or zero; for 
zero, drag is not applied. De- 
fault is 2.3 for satellites with 
a perigee height less than 800 
kilometers. The default for all 
other satellites is zero. 

Cq drag coefficient 


1 . 2-21 



' DRAG (Cont.) 


NAME 


COLUMNS 


FORMAT 


DESCRIPTION 


26-40 


D15.8 Standard deviation of drag 

coefficient or zero. If 
a non-zero standard devia- 
tion is specified the drag 
* coefficient will be adjusted. 

The default is - 0.5 for satel- 
lites .with a perigee height 
less than 800 kilometers. 

The default for all other 
satellites is zero. 


41-55 


D15.8 Value of drag coefficient rate 
of change per day. Default is 


zero. 


drag coefficient rate 


56-70 


D15.9 Standard deviation of in 
same units as C^. If a non- 
zero standard deviation is 
specified and if is greater 
than zero, then Cp will be 
adjusted . 


IF CARD OMITTED: Drag will be applied with default C D 

only for satellites with perigee height 
less than 800 km, 

NOTE: Cjj will not be applied or adjusted unless C^>0 . 


^Satellite number refers to either the first or second satellite 
and is applicable only th satellite-satellite trackine. 
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ARC SET 

GROUP 2 - OPTION CARDS 
. EARTH 
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KAMI: 

EARTH 


COLUMNS 




9-10 


FORMAT DESCRIPTION 

A5 Modifies size of gcopotentinl 

to be used for this arc. 

12 Maximum degree coefficient 

to be used .in the gcopotentia: 
Inodcl plus 1. 

Default = 23, Minimum = 4 

• • 

% 

12 Maximum r.Icr coefficient to 

be used in the geopotential 


1*2-23 



• .. EARTH (Cont.) . 

• . * . * 

• » 

NAME COLUMN'S FORMAT DESCRIPTION 

9-10(Cont.) - . model plus 1. 

•. - Default = 23, Minimum = 2 

. * '• . • •* ' ’■ A' 

■* t i‘\ •• ; ■ * ; -• 1 • . 

IF CARD OMITTED: Effects- of the geopotential will be 

. ■ . .applied using default values for degree 

*'••• . and order unless EARTH card is used in 

V- COMMON SET. 


Set-up for two-body orbit computations is explained on the 
following page. 


Sec also COMMON SET, OPTION CARD 


EARTH. 






o 


t 
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• EARTH (Cont.) 

* 

Two-body orbit computations may be performed in 
GEODYN by using the EARTH, BODIES, TIDES, and SAT cards 
as follows: j 

, * ' 

EARTH card 

Columns 7-8 =03 

Columns 9-10 = 02 

Columns 56-70 = 1. 

BODIES cards (two) 

. Card 1 - 

Column 7 =1 

Column 8 =2 

Column 11-25 = 0. 

Card 2 - ' 

Column 7=2 

t 

Column 8 =0 

Column 11-25 = 0. 

TIDES card 

Columns 7-80 blank 

one SAT card for each satellite 

Columns 11-25 = 0. * 



f 


. ARC SET 

GROUP 2 - OPTION CARDS 
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NAME 


COLUMNS 


FORMAT 


DESCRIPTION 


EBIAS 


1-5 


A5 


Requests removal of electronic 
biases from estimation orocess 


7-10 


14 


Station number for which elec 
tronic bias will be removed. 


11-25 

26-40 

41-55 

56-70 

71-80 



Measurement type numbers* 
. for which biases will be 
removed. 


*See Measurement Types and Sigma Defaults chart under 
COMMON SET, OPTION CARD - SICMA. 


1.2-26 


^SSSfSSi 





„ UCIUDC i 

. t 

. *w 

EBIAS (Cont.) 

f * 

IF CARD OMITTED: No electronic biases removed. 

. * 

• The value of the bias is not printed but the bias 
is removed from the estimation process and the statistics 
printed in the residual summary. An ' E f will appear in 
the column following the ratio to sigma for measurements 
used in electronic bias -estimations . 


No‘ a priori information is used. 

There is an independent bias for every pass. 
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1 

ARC SET 

GROUP 2 - OPTION CARDS 
. EDIT 



NAME COLUMNS FORMAT DESCRIPTION 

• . EDIT 1-4 A4 Removes from the estimation 

process any measurements 
with residuals greater than 
a specified number of stan- 
dard deviations. 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 
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EDIT (Cont.) 




NAME 

COLUMNS 

„ c 

FORMAT 


DESCRIPTION 


11-25 

D15.8 

This is 
plier.* 

the editing multi- 
Default is =3.5 


26-40 

D15.8 

This is 
Default 

the initial RMS. 
is E r = 200. 

• 

. 41-55 

D15.8 

First iteration edit criter 
ion for electronic bias 
passes. Default - use same 
.editing as other data uses. 

IF CARD 

OMITTED: 

Default editing, multiplier and RMS used. 


•Explanation : The product of the multipler and the RMS 

is the editing criterion. After the first 
inner iteration the RMS of all weighted 
measurements for the' previous iteration is 
used to determine the editing criterion. 
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* ARC SET 

tjROUP 2 - OPTION CARDS 
ELCUT 



NAME 


COLUMNS 


FORMAT 


DESCRIPTION 


ELCUT 1-5 AS 


11-25 D15.8 


REPRODUCIBILITY OE TRP 
ORIGINAL PAGE IS POOR 


Specifies the elevation 
angle below which measure- 
ments will be deleted 
from the solution. 

Elevation cut off angle 
in degrees. 

Default = 0° or value 
used on ELCUT card in 
COMMON SET . 


IF CARD OMITTED: All measurements with elevation angles less 

than 0° will be deleted from solution unless 
ELCUT card is used in COMMON SET. 

Elevation cut-off angle may be negative. 

See also COMMON SET, OPTION CARD - ELCUT. 
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reproducibility op THF 

NAME COLUMNS FORMAT DESCRIPTION ORIGINAL PAGE IS POOR 


MAXLAT 1-6 A6 Requests print out of the maximum latitude ob- 

tained by the satellite during each revolution. 
The OUTPUT card must be used in conjunction wit 
this option (see notes below) . 

7 II =1 Requests that the Kepler Elements be 

punched in PCE format (see Appendix C) at the 
times of maximum satellite latitude. 


11-25 D15.8 


Acceptable latitude tolerance in degrees 



26-40 D15.8 Tolerance for reduction of time interval betwcc 

successive latitude comparisons (units degrees) 

f 2- 


1 . 2 - 30.1 


MAXLAT (Cont.) 


COLUMNS FORMAT DESCRIPTION 

41-55 D15.8 Tolerance for switching back to the 

nominal latitude comparison interval 
(units degrees), 

IF CARD OMITTED: This option will not be invoked. 

NOTES: In order to determine the needed input for this option 
an orbit generator run should be made, outputting a 
Cartesian ephemeris, (which contains latitude as well) 
over one revolution of the satellite at intervals of one 
second (see ORBIT and OUTPUT cards) . can then be 

chosen as the difference between the two latitudes 
immediately preceding the maximum latitude value on 
one revolution. The output rates of the orbit generator 
and MAXLAT data reduction runs should be set equal. 

When the tolerance has been met the time interval 
between successive latitude comparisons becomes 1/10 
of one second. Near maximum latitudes are then defined as 
beginning when the difference between successive latitude 
comparisons becomes less than and ending when the 
difference becomes less than 

The normal satellite ephemeris printout as requested on 
the OUTPUT card is not generated under the MAXLAT option. 
Instead, only one line of Cartesian ephemeris is printed on 
each revolution at the highest latitude to within 1/10 
of one second; this latitude is determined when successive 
latitude comparisons within the period of maximum latitudes 
first become negative. Thus, in order to insure the proper 
functioning of the MAXLAT option, the tolerances on 
and cannot be too small. Suggested initial values for 
and c can be determined by: 

1.2-30,2 


1 


MAX L AT (Cont.) 

.1 X c 2 

€ = — , 

1 10 x OUTPUT RATE 

* c 3 = ’ ^1 

I£ the OUTPUT card requests the Keplerian ephemeris print- 
out, the Keplerian ephemeris over the period of maximum 
latitudes will be generated. If, in addition , the punching 
of Keplerian elements is requested on the MAXLAT card, the 
times of these punched elements will correspond to the 
lines of Cartesian ephemeris printout. 

Use of this option will result, in a considerable increase 
in program execution time. 

When this option is used the following must be true 
C 1 > 0 » C 2 > 0 9 C 3 < 0 


1.2-30.3 







' 0RB1 (Cont.) 


NAME 


Note 


COLUMNS 

FORMAT 

¥ * 

DESCRIPTION 


11-25 

D15.8 

Year, month, day 
(YYMMDD . ) 

Start 

26-40 

D15.8 

Hour, minute, seconds 

time 

•' 

• 

(HHMMSS.) 


41-55 

D15.8 

Year, month, day 
(YYMMDD.) 

Stop 

56-70 

D15 . 8 i 

Hour, minute, seconds 

time 

71-80 

DIO. 8 

(HHMMSS . ) 
Satellite ID** 



An 0RB1 .tape may be generated during a data reduc- 
tion run or an orbit generator run. The start time 
must be equal to or greater than the start time for 
data as specified on the Epoch card. 

*See Appendix C of this manual. 

9 

**I£ a satellite ID is specified on the 0RB1 card 
that ID will appear on the ORB1 tape instead of 
the ID on the Epoch card or the data. 
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ORBIT 

1-5 

AS 

Requests orbit generation 




mode of operation. 


7-8 

12 

FORTRAN- unit for writing 
an RV tape. If no unit is 
specified the tape will not 

. 


* 

be written. 


11-25 

D15.8 

Year, month, day (YYMMDD.) 
of start time for ephemeris 
printout and RV tape generation. 


26-40 

D15.8 

Hours, minutes, second (HHMMSS.) 
of start time. 















- . 

• 

# 

r 

ORBIT (Cont.) 

* 

NAME 

COLUMNS 

* t 

FORMAT 

DESCRIPTION 

• 

41-55 

• 

D15.8 

* 

• ' 

Year, month, day (YYMMDD.) 
of stop time for ephemeris 
printout and RV tape generation. 


S6-70 

D15.8 

Hours minutes , seconds 
(HHMMSS.) of stop time. 

NOTES: 

Use the 

OUTPUT card 

to specify the time interval for 


listing the satellite ephemeris. 

If an ORBIT card is present for any arc the maximum 
number of ‘outer iterations must be specified as 1 on 
the COMMON SET, Termination Card. 

Normal ephemeris output is inertial position and 
velocity plus groundtrack on the printer. The co- 
ordinate system is true of date. This output may be 
requested in true of the reference date by use of the 
OUTPUT card. 

An RV tape is optional and contains time in the UTC 
system, and position and velocity in the same co- 
ordinate system as the printer output. 

The start time for ephemeris output may be greater 
than or equal to the start time for data as specified 
on the Epoch card. ■ 

IF CARD OMITTED: Data reduction mode of operation will be 

requested.. 
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OUTPUT 
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NAME 


COLUMNS 


FORMAT 


DESCRIPTION 


OUTPUT 


1-6 


A6 Requests additional output 

or specifies the time interval 
for an orbit generator printout. 


7 


II 


reproducibility of the 
ORIGINAL PAGE IS POOR 


> 0 indicates that satellite 
XYZ and groundtrack are 
requested for iterations 
as follows , 

* 1 ;first inner on first outer 
(mandatory for orbit 
generator run) 

c 2 last inner on last outer* 
e 3 first inner on first outer 
and last inner on last outer 


1,2-35 



Output (Cont.) 


COLUMNS 

8 


9 

10 


11-25 


FORMAT DESCRIPTION 

II A value greater than zero 

specifies that a ground- 
track tape is to be written 
on FORTRAN logical unit 18. 

The groundtrack tape is 
written on the first iteration. 


II =1 requests a binary residual 

tape output on the last 
.iteration. Default is no 
tape. 

II =1 Requests that satellite 

« Keplerian element 

_ ephemeris be printed 

at the same times as 
the XYZ ephemeris. 

i 

D15.8 1. Specifies the time interval in 

seconds for orbit generator 
output. No default. 

2. A value in this field for a data 
reduction run indicates that the 
XYZ and groundtrack are to be 
printed at the specified time 
interval in seconds and printed 
at data times. Default is 
. output at data times only. 



f 


OUTPUT (Cont.) 

NAME COLUMN'S FORMAT DESCRIPTION 

* • 

26-40 D15.8 A value greater than zero 

references the ephemeris 
output to the coordinate 
system of the true equator 
and equinox of the reference 
time, 

41-55 D15.8 Optional change of unit 

• number for Keplerian ephemeris. 

: . Default = 9 . 

Unit may be changed to 6 or 8. 

56-70 .D15.8 Value greater than zero requests 

printing of measurement partials 
"during data reduction run. 

B 

IF CARD OMITTED: 

For orbit generator - There will be no printout. 

For data reduction -- No additional output requested. 
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ARC SET 

GROUP 2 - OPTION CARDS 
PREPRO 
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NAME COLUMNS FORMAT 


OFffiS 


DESCRIPTION 


PREPRO 1-6 


A-6 Requests data preprocessing. 


7 II Digit in column 7 requests 

transit time and transponder 
delay corrections as indicated 
in Table I on page 1.2-39.1. 

8 II Digit in column 8 requests refrac- 

tion corrections as indicated in 
Table II on page 1.2-39.1. 

9 II Digit in column 9 applies 

constant timing corrections. 


1*2-38 


Apu.1 o, Jiy/H 

PREPRO (cont.) 

DESCRIPTION 

Digit in column 10 requests 
optical data corrections as 
indicated in Table III on 
page 1 . 2 -39 . 2 .' 

Station number or zero; for 
zero, measurements from all 
stations will be preprocessed. 

Measurement type or zero; for 
zero all measurement types 
will be preprocessed. 

41-55 D15.8 Surface pressure in mm of Hg for 

tropospheric refraction computation. 
Default = 760. 

56-70 D15.8 Constant timing corrections in 

seconds. 

71-80 DIO. 8 Surface .temp in °C to tropospheric 

refraction correction computa- 
tions. Default = 20. 

IF CARD OMITTED: No data preprocessing is done. 


NAME COLUMNS FORMAT 

10 II 


11-25 D15.8 


26-40 D15.8 


1,2-39 
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TABLE I ♦ TRANSIT TIME/TRANSPONDER DELAY CORRECTIONS 


Value in 


Col. 7 

TT 

TTF 

TD 

VHF 

0 

NO 

NO 

NO 

NO 

1 

YES 

NO 

YES 

NO 

2 

NO 

YES 

YES 

NO 

3 

NO 

NO 

YES 

NO 

4 

YES 

NO 

YES 

YES 

5 

NO 

YES 

YES 

YES 

6 

NO 

NO 

YES 

YES 

7 

YES 

NO 

NO 

. NO 

8 

NO 

YES 

NO 

NO 


Notation for Table I: 

TT - transit time correction applied if data tape requests 
TTF - transit time corrections always applied. 

TD - transponder delay corrections applied if data tape 
requests. 

VHF - indicates VHF transmitter frequency assumed. 


TABLE II. REFRACTION CORRECTIONS 


Value in 
Col. 8 


0 

1 

2 

3 

4 

5 

6 

7 

8 


T 

I 

NO 

NO 

YES 

NO 

NO 

NO 

NO 

YES 

NO 

NO 

YES 

YES 

NO 

NO 

NO 

YES 

YES 

NO 


TF 

IF 

NO 

NO 

NO 

NO 

YES 

NO 

NO 

NO 

NO 

YES 

NO 

NO 

YES 

YES 

YES 

NO 

NO 

YES 


Notation 

T 

I 

TF 

IF 


for Table II: 

tropospheric corrections applied if data requests, 
ionospheric corrections applied if data requests, 
tropospheric corrections always applied, 
ionospheric corrections always applied. 


c 


1.2-39.1 


April 8, 1974 


TABLE III. PREPRO OPTIONS FOR OPTICAL DATA 


Value in 
Col. io 

Precession § Nutation 
only if data requests 

Parallactic 

Refraction 

Annual 

Aberration 

Diurnal 

Aberration 

0 

NO 

NO 

NO 

NO 

1 

YES 

SAO only 

NO 

SAO only 

2 

YES 

YES 

• NO 

NO 

3 

YES 

NO 

YES 

NO 

4 • 

YES 

YES 

YES 

NO 

5 

YES 

NO 

NO 

YES 

6 

YES 

YES 

NO 

YES 

7 

YES 

NO 

YES 

YES 

8 

YES 

YES 

YES 

YES 


SAO only means stations with numbers greater than 8999. 
Cards appearing later in deck will override earlier cards. 
Max. of 50 PREPRO cards per arc. 
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ARC SET # 

GROUP 2 - OPTION CARDS 
RECOEF 
(COEF) 
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NAME 

' COLUMNS 

FORMAT 

RECOEF 

1-6 

A6 

or 



COEF 

oo 
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9-10 

12 
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reproducibility oe the 

ORIGINAL page is poor 


DESCRIPTION 

Modifies and/or requests the 
estimation of coefficients 
in the geopotential model. 


Degree of C and S coefficients 
(n index) . 

Order of C and S coefficients 
(m index) . 



October 20 



RECOEF 
« «• . 

CCont.) 

COLUMNS 

‘ FORMAT 

DESCRIPTION 

11-25 

D1S.8 

A priori for starting) value 


* 

of C coefficients. 

26-40 - 

D15.8 

A priori (or starting) value 

• 

« 

of S coefficients. 

41-55 

D15.8 

Standard deviation of C co- 
efficient. 


Coefficient will not be 
adjusted if value is zero. 

Value is fractional uncertain- 
' ty if >1.0D-5 ( ^T7 od-*S for J2) 
and standard deviation equals 
fractional uncertainty times 
nominal* value of a priori 
coefficient. 

Value is standard deviation 
-if <1.0D-S (V 1.0D-S for J2). 

56-70 D15.8 Standard deviation of S co- 

. efficient. 

Coefficient will not be 
adjusted if value is zero. 

• Value is fractional un- 

certainty if >1.0D-5 and 
standard deviation equals 
fractional uncertainty times 
nominal* value of a priori 
coefficient . 

Value is standard deviation 
if <1. OD-5 . 
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RECOUP (Cont.) 

NAME columns format description 

71-80 DIO. 4 =0 Values on this card will be 

interpreted as being un-normal- 
ized. 

>0 Entries in columns 11-40 on 
this card will be ignored. A 
• priori values of coefficients 
are obtained from current geo- 
„ potential model including all 

modification prior to this card 
Input standard deviations are 
un-normalized . 

<0 Values on this card will be 
interpreted as being normalized 

IF CARD OMITTED: The geopotential used will be GEM 1 Model 

unless modified by EARTH card or modified in 
COMMON SET by option cards EARTH and REC0EF. 


The geopotential coefficient denormalization equation 
used is: 


,m _ 

t . “ 

C n * 

m 

n 

n 

n 

,m - 

s in A 

m 

n 

n 

n 


where 


m / (n-rc) 1 

D n = y (4n + 2) (n+m) ! for m i o 


>1 


1.2-42 



; October 20, 1972 


j * 

• RECOEF (Cont.) 

Djj = V2n + 1 for m = o 

F m 

- and 

(£, Si are un -normalized coefficients 
cjj, Si are normalized coefficients 
♦Nominal values of geopotential coefficients are 

- i-cl- < 

i 

A Priori Value of Adjusted Coefficients 


If it is desired that the starting values of geopotential 
coefficients be different than the a priori values for those 
coefficients then the RECOEF card should be used as follows: 

• RECOEF cards requesting the desired coefficient 
adjustments should be included in the setup 

. • .deck. The a priori values will be indicated 

on these cards. * 

• For each coefficient pair for which the starting 
value is to be different than the a priori value, 
an additional RECOEF card should be in the setup 
deck. This RECOEF card will indicate the start- 
ing value and come later in the deck than the 
corresponding RECOEF card requesting adjustment. 


* See also COMMON SET, OPTION CARD - RECOEF (COEF) # 


T 
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NAME COLUMNS 

FORMAT 

DESCRIPTION 

REGINT 1-6 

A6 

Indicates that numerical orbit integration 
is-to be performed in regularized mode. 


IF CARD OI^IITTED ; Integration will not be performed in regularized mode. 
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NAME COLUMNS 
RESID 1 - 5 


FORMAT DESCRIPTION 


AS 


REPRODUCIBILITY op TWP 
ORIGINAL PAGE IS POOR 


Controls residual printout 
for this arc. 


II 


0 - Indicates no residuals are 
requested for this arc. 


1 - Indicates that residuals 
are requested on the first 
inner iteration of the first 
outer iteration for this arc. 
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RES ID (cont.) 


NAME • COLUMNS FORMAT DESCRIPTION cont. 

2 - Indicates that residuals are 

requested on the last inner 

iteration of the last outer 

• - • iteration for this arc. 

? m 

3 - Indicates that residuals 

are requested on the first 
inner iteration of the first 
outer iteration and the last 
inner iteration of the last 
outer iteration for this arc. 

4 - Indicates residuals are 

requested on all iterations 
for this arc. 

IF CARD OMITTED; Value is 3 or that value specified on the 

RESID card in COMMON SET. 

See also COMMON SET, OPTION CARD -- RESID. 


Inner and outer iterations are explained in the Introduction 
of this manual. 
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NAME 


COLUMNS 


FORMAT 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 

DESCRIPTION 


-■ SAT 


1-3 


r> 

O 


A3 Introduces satellite area 
. and mass for drag or solar 
radiation pressure calcu- 
lations and satellite ID.* 


♦ » 

*. II • Satellite number (1 or 2)?* 

.. If less than 1 or greater than 

2, applies to both satellites. 


11-2S 


D15.8 Cross-sectional area of satelli 

in square meters. Default is 

* 

zero. 


f 


NAME 


* ’ 

SAT (Cont. } 

COLUMNS FORMAT DESCRIPTION 

26-40 D15.8 Mass of satellite in kilo- 

*• grams. Default is zero. 

41-55 D15.8 Satellite ID. 


IF CARD OMITTED : No drag or solar radiation pressure per- 

turbations will be applied even if requested 
* by OPTION CARDS -- DRAG and SOLRAD. 


*The satellite ID is necessary in the following 
situations . 

• If there is more than one satellite in the arc. 

m 

• If data from more than one satellite is on the 
data tape for the arc. 

• If the DODS Data Base is to be used. 


^‘Satellite number refers to either the first or second 
satellite in the arc and is only applicable in satellite - 
satellite tracking. 
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NAME 


COLUMNS 


FORMAT 


DESCRIPTION 


SIGMA 


. rwroducibility of tee 

AS Modifies TTlcasurem |^Q]^[yy J PAGE IS POOR 
standard deviations. 


11-25 


D 15.8 Station number or zero; for 
zero all measurements of the 
type specified in columns 26 - 
40 will be assigned the standard 
deviation on this card 


26-40 


D 15.8 Measurement type or zero; for 
zero all types will be as- 
signed the standard deviation 
on this card. 


r 


f 


SIGMA (cont.)* * 

♦ » 

NAME * COLUMNS FORMAT ‘ DESCRIPTION 

41-55 D15.8 Standard deviation of measure 

ment or zero; for zero, resi- 
duals will appear for the 

. • measurement but the measure- 

ment will not be used in the 
estimation process, o 

IF CARD OMITTED: All measurements are assigned Sigma 

Defaults. 

NOTE: If requests conflict in Sigma change decks the last 

/ of the conflicting changes takes precedence. 


See following page for Measurement Types and Sigma Defaults. 
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SIGMA (Cont.) 
MEASUREMENT TYPES AND 

t 

SIGMA DEFAULTS 



MEASUREMENT 

TYPE 

NO. 

GEOS DODS 

DEFAULT DEFAULT 

UNITS FOR 
STANDARD DEVIATION 

a 

Right Ascension 

1 

Tape 

2. 

Seconds of arc 

( 

Declination 

a 

Tape 

2. 

Seconds of arc 

P 

Range 

2 

Tape 

25. 

Meters 

• m 

p 

Range Rate 

3 

Tape 

10. 

Centimeters per 
second 

H alt 

Altimeter Height 

4 

Tape 

N.A. 

Meters 

fl alt 

Altimeter Rate 

11 . 

Tape 

N. A. 

Centimeters per 
second 

1 

Direction 
Cosine E 

5 

.3 

.3 

Mils 

m 

Direction 
Cosine N 

12 

.3 ' 

.3 

Mils 

V 

X angle 

6 

Tape 

50. 

Seconds of arc 

Y. 

Y angle 

13 

Tape 

50. 

Second of arc 

a 

A, 

Azimuth 

7 

Tape 

50. 

Seconds of arc 

E 1 

Elevation 

14 

• Tape 

50. 

Seconds of arc 




PCE 

DATA 






STANDARD 

DEVIATION 


MEASUREMENT 

TYPE NO. DEFAULT VALUE 

UNITS 


X 

15 


100. 

Meters 


y 

16 


100. 

Meters 


z 

. 17 


100. 

Meters 


• 

X 

18 


10. 

Centime ter s /Second 


• 

y 

19 


10. 

Centimeters /Second 


1 

20 


10. 

. * 

Centimeters /Second 


a 

21 


150. 

Meters 


e 

22 


IS. 

ppm 


i 

23 


10. 

arc seconds 


ft 

24 


10. 

arc seconds 


U) 

25 


20. 

arc seconds 


e 

M 

26 


SO. 

arc seconds 


SIGMA (Cont.) 


i 


-TWO STATION MEASUREMENTS 


TYPE DEFAULT 


MEASUREMENT 

NO. 

VALUE 

UNITS 

T 

‘VLBI Time Delay 

27 

1000 

Nanoseconds 

V 

VLBI Fringe Rate 

28 

100 

Microhertz 

• 

2-Way Average Range 
Rate 

f 

29 

0.1 

Centimeters/ 

Second 

P3 

t 

3-Way Average Range 
Rate 

30 

0.1 

Centimeters/ 

Seconds 
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March 


SIGMA (cont.) 

* 

NOTES: GKOS Default refers to the default if the data tape 
is in GEOS format. 

DODS Default refers to the default if the data is 

'• ■ I I !■ ) I— H » ■ . ■ ■ I 

in DODS format. , . • 

■ • . * • V ' / 

r, ‘ ■ 

Tape implies that the sigmas were obtained from the 
data tape. * . - 

The sigma for right ascension is divided by\the 

• . cosine of the declination. * 

. The sigma for azimuth is divided by the cosine of 
the elevation. * 

* * 

The sigma for direction cosines arc divided by the 
square root of 1.0 minus the measurement squared. 

The maximum number of SIGMA cards is SO per arc. 

- PCE sigma defaults are used only wlffcnr sigmas on 
tape are zero. 

Summed range and range rate measurement used for 
^ satellite- to-satellite tracking are handled 
internally as types 2 and 3 respectively. 
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ARC SET ' 

GROUP 2 - OPTION CARDS 
SIMDAT 





NAME COLUMNS FORMAT 

SIMDAT 1-6 A6 


OF THE 

ORIGINAL PAGE IS POOR 


DESCRIPTION 


Requests that data be simulated 
during the run and written on 
unit 17. 


If the run is in data reduction 
mode the data will be simulated 
on the first inner iteration of 
the first outer iteration.* This 
data will be of the same types, 
from the same stations and at 
the same times as the real data 
used in the data reduction run. 




1, 2- S3 



# 

SIMDAT (Cont . ) 


f 


s 


NAME COLUMNS FORMAT . DESCRIPTION 


•1-6 (Cont.) If the run is in orbit 

generator- mode PCE data will 
be simulated at satellite 
ephemeris output times. 

7 II |>CE -.type indicator if in 

orbit generator mode. 

*0 X,Y,Z,X,Y,Z data 
*1 a,e,i,n,w,M data 


IF CARD OMITTED: Data will not be simulated or written. 


NOTE: One file will be written on unit 17 for each arc 

in which a SIMDAT card is present. 


*See Introduction for definitions of inner and outer 
iterations . 


REPRODUCIBILITY OP THE 
ORIGINAL PAGE IS POOR 
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GROUP 2 - OPTION CARDS 
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NAME 


COLUMNS 


FORMAT 


DESCRIPTION 


SOLRAD 1-6 A6 Modifies solar radiation 

pressure application. 

7 D15.8 Satellite number (1 or 2).* If 

less than 1 or greater than 2, 
applies to both satellites. 


11-25 D1S.8 


PFPRODUCIBILITY of tiif 

ORIGINAL PAGE IS POOR 


Satellite reflectivity or zero; 
for zero, solar radiation pres- 
sure is not applied. Default is 
C R *1.5. 
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SOLRAD (Cont.) 


NAME 

COLUMNS 

FORMAT 


. 26-40 

D15.8 


DESCRIPTION 

Standard deviation of satellite 
reflectivity or zero; for non- 
zero standard, deviation the 
satellite reflectivity will be 
adjusted. Default is zero. 


IF CARD OMITTED: Solar radiation pressure is applied with 

C R = 1.5 only if SAT card is also present. 


NOTE: SAT card in ARC SET must be present for solar radia- 

tion pressure to be applied. 


♦Satellite number refers to either the first or second 
satellite in the run and is applicable only in satellite - 
satellite tracking. 
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GROUP 2 - OPTION CARDS 
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REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 


NAME 

STEP 


c 


COLUMNS FORMAT DESCRIPTION 

1-4 A4 Changes integrator step size 

and mode . 


7 . II 


«1 Satellite 
* Motion 
*2 Satellite 
Motion 
«3 Satellite 
Equations 
*=4 Satellite 
Equations 


1 Equations of 

2 Equations of 

1 Variational 

2 Variational 
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. STEP 

* 

(Cont.) * 

t 

m 

NAME 

COLUMNS 

FORMAT 

DESCRIPTION 

• 

.7 (Cont.) 

• 

■ * * 


■0 or > 4 applies to Equations 
of Motion and Variational 
Equations for both satellites. 


8 

11 

®0 Fixed step integration 
“1 Variable step integration 
®2 Variable step integration 
with stepsize 'halving and 
doubling only. 

i 


. 9-10 

12 

Integration order. Default - 11 


11-25 

D15.8 

Step size in seconds; the initial 
step size for variable step mode. 


26-40 

D15.8 

Error tolerance for stepsize 
-decrease. Default is 0.25D-4 m. 

- 

41-55 

D15.8 

Error tolerance for stepsize 
increase. Default is 0.25D-10 m. 

- ■ 

56-70 

D15.8 

Maximum stepsize in variable 
step mode. Default is 400 
seconds. 


71-80 

DIO. 5 

Minimum stepsize in variable 
step mode. Default is 12.5 
seconds. 

IF CARD 

OMITTED: 

Default values will be used. 

© 
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• 

r 



. ( 


STEP (Cont.) 

« * 

' f » 

NOTE: The step size and integrator mode defaults are the 

following: : 

If perigee height is less than 800 kilometers the 
Step. size is 75 seconds. 

If perigee height is greater than 800 kilimeters 
and less than 1500 kilometers the step size is 
100 seconds . 

The step size is linearly interpolated from 100 
seconds at 1500 kilometers perigee to 400 seconds 
at 36,200 kilometers perigee for satellites with an 
eccentricity less than 0.01. 

If the perigee height is greater than 36,200 kilometers 
-the step size is 4.00 seconds for circular orbits. 

If the eccentricity of the orbit is greater than 
0.2 the integration is done in variable step mode 
and the initial step size is 50 seconds. 

If columns 26*40, 41*55, 56-70. or 71-80 are blank 
or zero, default values will be used for the corresponding 
parameters. The default values may be modified for all 
arcs by use of a STEP card in COMMON SET. 
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ARC SET 

• GROUP 2 - OPT i ON CARDS 
TERMEL 
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NAME COLUMNS 
TERMEL 1-6 


FORMAT 

A6 


DESCRIPTION 

Requests -printing of tcr:~: 
elements after end cf clr.tr: 
a data reduction run. 


11-2S 


D15. 8 


Date of terminus elements 
YYMMDDi - 


26-40 ^ D15.8 . * Time of terminus elements 

HUMMSS. 

» * ' ' • 

IF CARD OMITTED: Terminus elements u 11 not be printed. 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR ,. t 
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ARC SET 

GROUP 2 - OPTION CARDS 
TOLS 


*S 
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NAME 

TOLS 


COLUMNS 

1-.4 

7 


r> 

0 


FORMAT 

A4 


II 


REPRODUCIBILITY OR THU 
ORIGINAL PAGE IS POOR 


9-10 


12 


DESCRIPTION? 

.Modifies integration tolerances 

Indicates for which satellite \ 
tolerances are to be modified 

for this arc« 

• » 

* 0 both* satellites 

** 1 first satellite ' 

® 2 second satellite 

(applies to satellite - 
satellite tracking) 

Maximum degree and order of 


I 


TOLS (Cont.) 


NAME 


COLUMN'S 
9-20 (Cont.) 


FORMAT 


*/. •• 


v • ••• 


DESCRIPTION 

gcopotcntial effects to be 
included in the computation 
of the variational equations, 
Default = 4, Minimum = 3. 


11-25 


D15. 8 ‘Desired reasonable predictor* 

; corrector tolerance to be 
* achieved by changing, the intc- 
- . grator step size while in 
• .* variable step integration mode, 

Oscc\ only in computation of 
new step size during step size 
variation. 

• . - Default = 0 . 2 5 D - 7 hi . 


26-40 


D15.8 


Critical predictor- corrector 
tolerance. Used to determine 
number of corrector iterations 
necessary at each integration 
step. If after two corrections 
the local error exceeds critical 
•tolerance an error message will 
be printed (execution will 
be terminated). 

Default = 0 . 25D-10ia, 


IF CARD OMITTED: Default values will be used. 


'fflnjjy 0F TH 
page is poor. 
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ARC SET 

GROUP 2 - OPTION CARDS 
VARCOV 
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NAME 

VARCOV 


COLUMNS 

1-6 


FORMAT 

A6 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 


DESCRIPTION 

Specifies that the a priori 

variance -covariance matrix (^^3 
associated with the a priori 
orbital parameters is to be 
changed. 


II Satellite number (1 or 2).* 

Indicates for which satellite 
matrix will apply. If zero 
or greater than 2, matrices 
must be input for both satellites 


8 


II -0 V ariance - covariance on epoch 

Cartesian elements. 

»1 Variance -covariance on epoch 
Keplcrian elements. 
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'VARCOV (Cont.) 

NAME COLUMNS FORMAT • DESCRIPTION 

9 II =0 Matrix diagonal only will 

be input on one card in 6E12.6 
format for each satellite for 
which the matrix is changed. 

. * This card(s) follows the VARCOV 
card. 

* 

Full 6x6 matrix will be input 
on 6 cards in 6E12.6 format 
for each satellite for which 
the matrix is changed. These 
cards follow the VARCOV card. 


IF CARD OMITTED: No change to matrix V^ will be made. 


^Satellite number refers to either the first or second 
satellite in the arc* 


r 



. t 


The a priori variance 

* 

VARCOV (cont.)' 
-covariance matrix 

is : 

s 

l.E+14 0. 

0. 

0. 

0 . 

0 . 

0. l.E+14 
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0. 

0. 

0 . 
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l.E+14 
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0 . 0 . 

0. * 
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1.2.3 ARC SET 

GROUP 3 - TERMINATION CARD FOR GROUP 2 - OPTION CARDS 

DATA 

There are four kinds of termination cards used to signify the 

end of Group 2. One of these is mandatory for every arc. 

- • * 

1. Termination Card for orbit generator arcs. 

NAME COLUMNS FORMAT DESCRIPTION 

DATA 1-4 A4 Signifies the end of Group 2. 


2 . 


Termination Card for data’ reduction arcs with data in GEOS 
format. 


NAME 

COLUMNS 

FORMAT 

DESCRIPTION 

DATA 

1-4 * 

A4 

Signifies the end of Group 2. 


a 

11 

Specification of a 1 or 3 indi- 
cates that a GEOS Data Selection 
is to' be made and that the corres- 
ponding GROUP 4 SELECT/DELETE 
cards must be present. 




Specification of a 2 or 3 indicates 

that the data used in the arc he 
* 

written onto unit 3 in GEOS II forma 
# 




TERMINATION CARD FOR GROUP 2 (Cont.) 


NAME COLUMNS FORMAT 
9-10 12 

11-25 D15.8 
26-40 D15.8 
41-56 D15.8 


DESCRIPTION 

Optional change of input unit 
for GEOS data tapes. Default 
is 20. 

Additional input units for GEOS 
data tapes. If more than one 
input unit is specified the data 
will be merged in time order as 
it is read. 


Termination Card for data reduction arcs with data in 
GEOS-C Card Image Format, RAP Binary Format, or DODS 
Binary Format. 


NAME COLUMNS FORMAT DESCRIPTION 

DATA 1-4 A4 Signifies the end of GROUP 2. 


II. A value of 1 or 3** specified 

here signifies that data 

Selection/Deletion cards follow 

for data in any of these formats. 

GEOS-C Card Image Format*** 
RAP* Binary Format 
DODS Binary Format 

A value of 2 or 3 indicates that 

the data used in this arc be 

written onto unit 3 in the same 

format that it was read. 




TERMINATION CARD FOR GROUP 2 (Cont.) 


NAME 

COLUMNS 

FORMAT 

DESCRIPTION 


9-10 

12 

Input unit for data tape. 
Default is 20 for card image 
formats and 21 for binary 
formats. 


11-25 

D15.8 

Additional data tape input 


26-40 

D15.8 

^ units. If more than one unit 



| 

is specified the data will be 


41-55 

D15.8 

j 

l merged as selected. 

4. 

Termination Card 

for data 

reduction arcs with data in PCE 


or simulated data 

formats 

♦ 

NAME 

COLUMNS 

FORMAT 

DESCRIPTION 

DATA 

1-4 

A4 

Signifies the end of Group 2. 


9-10 

12 

The data tape input unit is 


specified here. Defaults are 

20 for PCE data and 

21 for simulated data. 
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1.2,4 Arc Set 


GROUP 4 - SELECTION/DELETION CARDS 


1. 2.4.1 Selection of Multiple Data Forms 

This Selection/Deletion option is applicable to the 
following input data forms: 

GEOS-C Card Image 
RAP Binary 
DODS Binary 

If the data Selection/Deletion option has been specified 
on the card then at least one SELECT card and as many as 100 
■SELECT and DELETE cards will follow the DATA card. The end of this 
card group is indicated by an ENDALL card. 

The form of the SELECT card is: 

Name Columns Format Description 

SELECT 1-6 A6 Selects from data tape. 

7 II Measurement type category indicator. 

1 - Right Ascension (or Hour 

Angle) and Declination 

2 - Range or Range Difference 

3 - Range Rate or Range Rate 

Difference 

4 - Altimeter 

5 - Direction Cosines 

6 - X-Y Angles 

7 - Azimuth and Elevation Angles 
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Name 


Columns 


Format 


Description 


8 A1 Measurement Sub-Type Indicator 

Sub -types vary with measurement 
type cateogry and will not in general 
-need to be specified. A complete 
list of the categories and sub-types 
is specified in Appendix C; columns 
8-9 of the GEOS-C Card Image Format. 

9-10 12 Measurement Modulo Number 

An integer n specified here indicates 
that only every nth data point from 
• the stations and types specified on 
this card will be selected. 


11-15 

15 

Receiving station number if data 
either the GEOS-C format or the 
format. 

is in 
RAP 

11-16 

A6 

Receiving station name if data is 
DODS Format. 

in 

18-24 

17 

Satellite ID 


25-30 

16 

Date in form YYMMDD for start of 
data selection. 


31-36 

16 

Time in form HHMMSS for start of 
data selection. 


37-42 

16 

Date in form YYMMDD for end of 
data selection. 


43-48 

16 

Time in form HHMMSS for end of 



data selection. 
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1.2.4. 2 SELECTION AND DELETION CARDS -- GEOS FORMAT- 

• • 

If Data Selection and/or Deletion is desired 
at least one and as many as 100 SELECT and/or DELETE cards 
must follow the DATA card. An ENDALL card must terminate 
this group of cards.* . 

* SELECT 

. . * 

The SELECT card is a Data Selection card for data tapes 


in GEOS format . 

NAME COLUMNS . FORMAT 

SELECT .1-6 A6 

7-12 16 ' 

14-15 12 

16-19 14 

20-23 14 

25-30 16 

32-33 12 

34-37 14 

38-41 14 


DESCRIPTION 

Selects data from tape. 

YYMMDD* of start time for 
data selection. 

Hours of start time for 
data selection. 

Minutes and integral seconds 
of start time for data selection. 

Decimal seconds of start 
time for data selection. 

YYMMDD of stop time. 

Hours of stop time. 

Minutes and integral seconds 
of stop time. 

Decimal seconds of stop .I me . 
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SELECTION/DELETION CARDS (Cont.) 

# 

NAME • COLUMNS FORMAT DESCRIPTION 

43-47 IS Station number; all stations 

. if not specified, 

* 

49 -II Measurement type; all types 

if not. specified. 

50-51 12 An integer n specified here 

indicates that only every 
. . n** 1 data point from the sta- 

tions and types specified on this 
this card will be selected. 

/ / These columns have no meaning 

for the DELETE card. 

/ 

This card is optional. 


*If start and stop times are not specified, these times are 
assumed to be the start and stop times of the data on the 
tape. 
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DELETE 


-The DELETE card is a Data -Deletion Card for data tapes 
~-in GEOS format. 

NAME COLUMNS * FORMAT DESCRIPTION 

DELETE 1-6 A6 Same format as SELECT except 

that all times, types, and 
stations refer to data del- 

’ * etion rather than data selection 

» 

This card is Optional. 


ENDALL 


NAME 

COLUMNS 

FORMAT 

• DESCRIPTION 

ENDALL 

1-6 

A6 

Specifies end of data selection 

• 

1 % 

• 

and/or deletion cards. 


This card is Mandatory when SELECTION OR DELETION cards for 
data tapes in GEOS format are used. 


1 


FORMAT 

Key n 
b 

FORMAT 

In 

Fn.m 


Dn.m 


En.m 

An 


SPECIFICATION . 

% 

« 

m * integer numbers 

* . a blank space 

• ■ 
m 

CODE 

Integral • numbers right-adjusted in a 
field of n columns. 

Example: 25 in an 13 format: b25 

Digits with a decimal point anywhere 
in a field of n columns; digits punched 
without a decimal point assume the 
point m places from the right. 

Example: 30.1 in an F5.2 format: 

b30.1 or 30.1b or b301b 

Digits with a decimal point anywhere in 
a field of n columns; digits right 
justified in a field of n columns with 
an exponent of the form D+_XX, where XX 
is the power of 10. If there is no 
decimal point, assume the point m placed 
to the left of the D. 

Example: 30.1 in a D8.5 format: 

. b30.1bbb or 30.1D+0C or 301.D-01 

* 

Same as Dn.m except that the exponent 
when used, takes the form Ej^XX. 

n alphanumeric characters. 
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2.0 GEODYN JOB CONTROL LANGUAGE AND HARDWARE AND SOFT- 
WARE RESTRICTIONS 

J l 

* < 

2.1 Job Control language (JCL) 

Certain JCL cards are necessary for all GEODYN runs. 

These cards are: . • 

• The EXEC card specifying the name of the pro- 
gram to be executed, 

• The STEPLIB card specifying the program library 
in which the GEODYN load module is stored, and 

i 

• The DD cards specifying the characteristics of 
the FORTRAN Logical Units used by GEODYN. 



. DETERMINATION OF GEODYN CORE REQUIREMENTS 


"Since GEODYN ’s core storage requirements are dynamic, 

"the following equations -should be -used to compute the REGION 
to be specified on the EXEC statement. 

REGION = + 330 

1024 

• * 

Bytes = 226 ♦ 46N, + 146N 2 + 224Nj (N 2 + Nj) + 96N 3 N 13 

4 ISONj + 214N 4 + 4Nj + 136N 6 + 16N ? +8N 6 N g + 16N g 
4 8N g + 56N n + 40N 12 + 18N 13 + 28N 1S + 4N* S + N 19 + N 


where 


N. 

k! 

i 

N 

n 

*1 

N. 

N 


8 



N 

N 

N 


11 

12 

13 


N 


14 


maximum number of biases in any one arc. 
number of tracking stations used, 
maximum number of satellites in any one arc. 
number of tracking stations estimated, 
number of surface densities. 

- number of surface densities estimated, 
number of constraint equations for surface 
densities . 

Ng - N d (N d + 2) where N Q is the maximum degree 
of geopotential for constraint equations as 
Specified on the TOLS card in COMMON SET, 

GROUP 2 - OPTIONS CARDS. 

number of solid earth tides coefficients 
estimated (not a current capability), 
number of geopotential coefficients estimated 
as common parameters. 

N io + N i6 ; 

N 7 4 Ng - Ng ' 

maximum number of equations xchich must be 
•integrated in any one aTC. ‘(=N 3 + Maximum 

of t N eql * N eq2» 

number of master tracking stations estimated. 



f 


where 




N 


17 


N 


18 


19 

20 


N 

N, 


21 


maximum number of parameters estimated in any 
one arc plus (number of common parameters 

-estimated) . 

maximum number of geopotential coefficients 
estimated as arc parameters in any one arc. 

N 10 +N 14 +N 12 


greatest number of two station (i.e. f 3-way 
average range rate or VLBI measurement) base- 
lines associated with a single station. 

(N 2 + N 3 ) (N 3 ) (N 18 ) (23) 

maximum number of electronic biases in a 
single arc. 1 

(2) (4N 1S + 13) (N ig ) 

maximum for any arc of the following 

N S “ N INT1 + N INT2 

N INTi = 12 + 12 N eqi + 6 ^°Ei ' ^ *Vi 

+ 6 (0 yi - 1) I yi N eqi - (I yi - 1)(1 + N eqi ) 




integration order for the equations of motion 
of the i satellite. 

integrtion order for the variation equations 
of the i t ^ 1 satellite. 

1 if fixed step integration is used for 
i th satellite. 

2 if variable step integration is used for 
i** 1 satellite. 

number of force model equations which must be 
estimated for the i ^ satellite. 

• orbital elements are six equations. 

• drag, drag rate, and solar radiation are 
one equation each. 

• geopotential coefficients are one equation each. 

• one equation for each 

0 for orbit generator run. 
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CEODYN JOG CONTROL LANGUAGE 
FOR THE 36C/S5 


August' 11 f 1975 


✓✓•• « JOS# • • 

✓✓GEODYN EXEC PGM = ZCT VMGDN • REG I CN— 3B0 K 

✓/STr.PLIO 00 CSN=C-i VG IUMP#CECSTR,DI SP = SHR 

✓✓FT'MFO'U DO DSN = G WZCTVH# Ef HE.V»DI 5P = SMC,OCP=DUFNO = * 

✓✓FTC2F001 OC DSN=PC£J$,CCE=<RECFr/=F, 81*51 Zn = -520#GUFNC-i )*UNIT=23J4* 
✓✓ V0L=SERsC00Sy2 . «PAC£= (TR*. 33»>4 > .0 ISPsSHH 

✓✓FTG?F*;'l OD UN 1 T = 2* t 9 . CCE* <R£CFM= VOS *LRECL=1 r*4 #CLKSI ZE=1044) « 

✓/ LABEL - ( #ELP}»VOL = SER = GLir.ATA 

✓✓ FT 14 FT "V 1 OC 0SN = PtLS.0C5 = <FECFM=F.a LK!:1 2E=ZS20 *CUFNC=t >.UNIT=23S4. 
✓✓ VGL=S E3=COD302» SP ACE= (IRK# 221 ) #01 SP=SHP 
*✓✓ FTtEFC^I DC DONAMU-CAT AS 

✓✓FTOG FMJ DO SYSQUT=A , DC S - ( FCCFM— VO A #LRECL=1 37 , GLKSl ZE = 7?.65 ) • 

✓✓ SPACE=<CYL, 20 > 

V/FTO7F0^1 00 SYSC]UT=B,DC8=(RcCFf'=FB*l.R£CL=QC#ELKSIZe = 3DC )• 

✓✓ SPACE- < SCO. 1 JC-) 

V/FTC'SFCO 1 DO SY50LT=A #DCE = <RECFM=VBA «LPECL=127 « EL KS I ZE=726 5 > « 


// SPACE MC yl.1? ) 

„//FT^9F0‘M OD S YSOU7 = A . OCe = ( PECFK=VBA #LPECL=137# ELKS I ZE = 7 2 65) • 

✓✓ SPAC€ = (CYL. 10 I 

✓✓FT10FCM 0 0 UN I T = DI SK,DCE= (RECFM=F5*LRECL=46 *QLKS I ZE=463 )* 

✓✓ SPACE=(TRK, 3 ) 

✓/FT1IFC01 00 UN IT=0ISK«DCE=<RECFM=veS*LFECL=68*EtKSIZE=7254), 

✓✓ SPACE-CCYL.3 ) 

//FT! 2FC*0 1 OD UN 17= 23 14 ,DCB= (GLKSI ZE=?232 ,0UFNC=1 ) , SPACE “( CYL « 2G ) 
✓/FT13FC'U OD UN I T=D1 SK,OCE=(fiECFM=VS #PLKSl 2L = 1 632 # BUFNC— 1 )• 

✓✓ SPACF-( C YL ♦ 10 ) 

✓✓FTiAFf-r 1 DO UN iT-CISK.DCE-(RECFM=V3 .BLKSI ZC = 1 CC4 ♦ C*J. c iiC*l }« 

✓✓ SPACE=( C YL, 5 > 

✓✓FT15FP?1 0 0 UNlT=Z3l4,DCE=(DLKSIZE=ien? ,QUFKO = 1 ) * SPACE = ( C YU . 2 > 

✓✓FT1 £FCO 1 0 0 UN £T=23l4, DCE- (RECFM=VGS*LFECL=1 804,eLKSIZE=l£7B# 

✓✓ OLFNO =1 ) .SPACE=(C YL, 5 } 

V/FT 17F01 1 DO UN IT=OI SK , DC E s <R E CFM= VB £ «LRECL=IOO ,GL KS I ZE = 3£27 ) • 

✓✓ SPACE= ( CYL « 2 )#D13P=( *PAS£> ,0£N=GS I FCATA 

/✓FTiEFCCl DC UNlT = 0ISK,0Ce-<RECFM=F8 *LRFCL=eC* SLKS I ZE = 1 30 0 ,3UFN0=1 ) « 

// S°ACE-( C YL. 1 1 ,D I3P=< , PAS 5) , D £N= &G RNCT RK 

//FT19FCM DO UN IT- 01 SK,OC6= (RECFM=V0S,LRECL=112*BLKS IZE-1 SC-4 , eUFNO=I ), 
✓✓ SPACE = (*CYL,1 ) #D ISP=< , PASS) ,OSN=CGEORGE r 

✓/ FT 2 0 F 0 0 1 DO UN I T = 24O0— 4 » C C E- ( R£CFM= F C #LRECL= 30 • EL KS I ZE = OOO'i , CEN = 3 ) • 


✓✓ LADEL= 
✓/FT21F0O 1 
✓✓ LABEL = 
✓✓FT2 2F70 1 
✓✓ SPACES 
✓✓FT2CF0C 1 
✓✓ SPACE 
✓✓FT7lFC.r 1 
✓✓ DEN=3 ) 
✓✓SYSLDUMP 
✓✓OAT AS OC 


( , ELP ) # VUL=SER = GEC 
DO UNlT=24y A -9,CC 
( , GLP ) . VCL=SER=3CC 
OD UN t T= Cl SK . DCE = 
(2eC4,SG) ,D ISP=( « F 
CD UN IT =2314 ,CCE 
= < 2 5176,1 C) 

DO UN I T s 2 4 C O— 4 * C C 
#LAOEL=( t ,OLP) , VCL 
DC SYSOUT=A, SPACE 


S D AT A 

E= (KcCFMsVCS .LPECL = 134 *eLK3 I ZE = 1044 > • 
£D A T A 
(PECFMsV: 

ASS ) ,DSN: 

=( RECFM=i 


S ,eLKSI ZE=2E04#OUFNO=I >, 
I=CCFP1 

: UT ,ELK-5lZE = 25l52 } • 


E” < RECFM 
=SER=BMA 
= ( C YL# 1 I 


= ves *LRECL=0S1G#0LKS I2E=SD2P# 

TFX T (S? 

nT .\3C®^l is POOR 


CEOOYN SETUP GOES HERE 



I 


THE FOLLOWING JCL C 05 4P£ K*ZZ~>- 0 FC* 
THE 360/51 IN CLAC£ C c THE 
STEPLIB C FTCIFiCl CAPDS 


/✓STfcPLJl CO DSNsM2.ZCTVM.LC4 OL Tti .SI 3 Ps3hQ 

■✓/FT 01 F C 0 l 00 ZCTVM.EPI-E M.Dl S^sSHR.O CBsEUFNO=i 






REPRODUCIBILITY 
ORIGINAL PAGE IS 


OF THE 
POOR 
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Name Columns Format Descript ion 

49-50 F2.2 Relative Humidity*at Receiving 

Station (subsatelli te point for 
. altimeter) 

51-53 13 Surface pressure in mm of Hg at 

Receiving Station (subsatellite 
point for altimeter) . 

54-56 13 Surface temperature in °C at 

Receiving Station (subsatellite 
point for altimeter) - 

57-61 15 Second station number if data is in 

either the GEOS-C format or the 
RAP format. 

Transmitter for average range rate 
‘ Reference for VLBI 
Relay for SST 


57-62 A6 Second station name if data is 

in DODS format. 

Transmitter for average range rate 
Reference for VLBI 
Relay for SST 


64-70 


17 Satellite ID for Relay Satellite in SST 


• - 71-72 F2 . 2 Relative Humidity* associated with 

second station. 

73-75 13 Surface pressure in mm of Hg at 

second station. 


76-78 13 Surface temperature in °C at ♦ 

second station. 

REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 

tfOTE: ^Relative Humidity - If 40$ input 40 

'■* - If 100$ input 1. 

1 . 2 - 67.2 

~r 



✓CL FOR GEOOYN UNIT 1 WHEN CENTER 
CARO IS USED 


✓✓FTC IF 1 OO UN IT *2 ACC*— 4 • CC E= <RcCFM=VDStLRECL=9T? « BL KS IZE-7294* 
✓✓ DEN- 2 } * LA8EL=( «NL)« VCL=£ER=XXXX 



GEODYN FORTRAN LOGICAL UNITS 

. ' V 

UNIT SPACE PURPOSE ' DESCRIPTION 


Ephemeris 


GEODYN uses the JPL ephemeris 


2 DODS Data • 

Base 

Observation File 


3 Data Tape 

Output 


4 DODS Data 

Base 

Element File * 


5 Card Input 


6 (CYL, (10,1)) Printer’ 


REPRODUCIBILITY of the 
ORIGINAL PAGE IS POOR 


7 (CYL, 1) Punched 

Output 


8 (CYL, 10) Printer 


for the moon, sun, Venus, 

Maxs , Jupiter, and Saturn . 
Necessary in all runs 

On 360/9S only. See Appendix 
A. 

See Arc Set, Group 3 - Termination 
Card of Group 2. 

On 360/95 only. See Appendix 

A. 

System card input unit; used 
for GEODYN Input Cards. 

System printer output unit. 

10 cylinders of primary and 
15 cylinders of secondary 
space are allocated. 1 
cylinder will hold 1060 
full lines of print. 

GEODYN punches updated para- 
meter values. The one cylinder 
of primary space which is 
allocated will hold 1800 
punched cards. 

Used for satellite Cartesian 
ephemeris printout during data 
reduction runs. 
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♦ 

* .UNIT 

space . 

PURPOSE 

• ti 

• 

• o 

n 

* ' m 

. , • _ DESCRIPTION' 

9 

• 

(CYL, 10) 

Printer 

-Used for satellite Keplerian 



# 

. m • 

„.ephcmeris printout. 

10 

(TRK,3) 

Scratch 

• 

Necessary for all data reduction 

* * * ♦ 



runs. ; Allocated primary space 




of 3 tracks will hold 4 95 

# 

. 

*■ ‘ * • * m 

. • * . 

• records. GEODYN writes a max- 


♦ •; * - ' 

• ♦ 

© r . 

imum of 000 records. 

• • , 

11 

(CYL, 3) 

Scratch 

« 

Necessary for all data reduction 



« 

runs. Allocated primary space 


. . 

* 

of 3, cylinders will hold ?060 

• 


* 

-observation data records. Suf- 



• 

ficient space for the maximum 

. 

* 

• 

number of data records in a 

• 

• * 


single arc must be allocated. 

12 

(CYL, 20) 

Scratch 

Necessary for all data reduction 

= ... 

REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR . 

runs. 20 cylinders of orir.arv 

• < 

. space arc allocated. 1 cylinder 
of space will hold 17 20 . 

. 

« 


servation data records. 

13 

(CYL, 10) 

Scratch 

Necessary for all runs. 10 


. • 

• » 

cylinders of primary space 

4T 


% 

are allocated. 1 cylinder 

. 


** 

of space will hold SO records. 

* 

* 

* 

Each arc writes five records. 

14 

(CYL, 5) 

Scratch 

* 

Necessary for all runs. S 


• 


cylinders of primary space 


‘are allocated. 1 cylinder 
of -space Will hold SO records. 


« 

• 

The mmbor of records per arc 



is variable. 



UNIT 

SPACE 

PURPOSE 

15 

(CYL,2) 

Scratch 

* • 

9 

16 

(CYL.S) 

Scratch 



Simulated 

Data 

Output 


18 


Groundtrack 


19 Binary 

Residual 

Tape 

• reprqducibility of the 

ORIGINAL PAGE IS POOR 


20 


Data Tape 
Input 


21 


Data Tape 
Input 


DESCRIPTION 

Necessary for all runs. 2 
•cylinders of primary space 
->re allocated. 1 cylinder 
of space will hold 80 records. 
The number of records per 
arc is variable. 

Necessary for all runs. 5 
cylinders of primary space 
are allocated. 1 cylinder 
of space will hold 80 records. 
The number of records per 
arc is variable. 

GEODYN can simulate tracking 
data or PCE data on option. 

If a groundtrack tape is 
-requested, JCL will be 
required for unit 18. 

A DD card is necessary for 
each file to be written. 
GEODYN writes one file for 
each Binary Residual tape 
requested. 

GEODYN default unit for input 
data tapes in GEOS format 
or PCE format. , 

GEODYN default unit for input 
data tapes in DODS format or 
Simulated Data format. 
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UNIT 

22 


30 


71 


23-29 

31-70 

72-99 


August 11, 1973 


SPACE 


(29176, 

10 ) 


PURPOSE 
ORB1 Tape 


Partitioned 
B -Matrix 
Scratch 
File 


- DESCRIPTION 

.A DD card is necessary for each 
file to be written. GEODYN 
writes one file for each 0RE1 
tape requested by an ORB1 card. 

Necessary only for those runs 
where the BPART option is used. 
Space allocated is sufficient 
for 10 records to be written. 
B-Matrix size in bytes must be 
less than the number of allocated 
recrods times 29176. 


B-Matrix 

tape 


A DD card is necessary for each 
file to be written. GEODYN writes 
one file for each arc exercising 
the BMATRX -option card. 


Option I/O These units are not regularly 
units used by GEODYN and may be used 

as additional data tape input 
files or as RV tape output files. 


REPRODUCIBILITY op TOr 
omginal page b°S e 
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-THE -6EODYN OVERLAY STRUCTURE IS 


August 11,. 1973 


OVERLAY ESTIM 

INSER T EST IM1 »E5T I M, PCSVEL 
OVERLAY ESTIM 

INSERT BFORMl .OFURM, BSCLVE .DREAD. F.fcRITE 
OVERLAY JACCH l A<PEG ION ) " 

INSERT DENSTY 
OVERLAY J7SJ65 
INSERT DEE 
OVERLAY J71J65 
INSERT D 6 50 
"OVERLAY U71J6S 
INSERT D 7 1 

OVERLAY PR IMEtREG ION > 

INSERT A CFlUX , JA.NTHO. AL 1ST .RFTMCD .FLUXM, FLUXStCMOCEL* GMOOEL 
INSERT INTFLX .FLUXAP 
OVERLAY PRIME 

INSERT 0 I AS * I NO LPT .NEfcAKC. RKNTPR « CCEFL * S I G8LK 
INSERT SRFBLK , POEM. E2 IAS. NCSTA«,FARTQ 
OVERLAY CAT AR E AO 

-INSERT DCOSRD .DATBSE . D COELF .OCOO AT 
OVERLAY DAT AR E AC 
INSERT GFQSRC 
OVERLAY SATCLCS 

INSERT SATCLC * 

OVERLAY SATCLCS 

INSERT S ATCL2 *SATC?1 . SATC22 
OVERLAY CAT AREAS 
INSERT PCERD 
OVERLAY OATAREAD 

INSERT SIMkD ' . 

OVERLAY PRIME 

INSERT NCNAME .CHARLY .CBRCWN.SQOANT .SQANT?. , CSTh ET . XYZ ♦ MCC3NGR . DRGOLK 
INSERT PR0CS1 .PROCES.CCVAO* ,CO‘*ELL .0A7APD .CATROl • O'* AG, CRAG: 

INSERT TRUPOL . SUMMRV . C L \ F.CT • GNOTRK’. X Y23U T , C ST AT , GFH«~ A 1 * G^M'iAN 
INSERT STAIFl , STA InF .SMSTaT .TRUER. TRUEST , V£ VAL . VCV AL 1 .VMAT 

INSERT ST A 1 F ?. • F *F1 .OBSDCT. C E SDTt , O = « I T .OF 3 17* 

INSERT PRECCT .PRCDC * » RE £ PA R . RES PR 1 .CGRREL . COM AO J 
INSERT SCRDEN .SLRON . .PCEN. FDEN2 

INSERT TfcOSTA .T*Q ST • . R A S CMP .UPDO WN • eSCOMF , 3S CFP1 
INSERT REFION.kFJCNI 
OVERLAY SECCNO 

INSERT C CMPAR . STA INP. ARCFAR. STAPOS 
' INSERT AREAS « CELTAZ, I NCE NT 
OVERLAY SECOND. - 

INSERT ECRAV.0R8TAP.SUNGRV.7IDAL.PEFCCR 

INSERT CCEF .COM . INT.'iP. INTERF.CSLI F 


OVERLAY TH IPO 
INSERT START 
OVERLAY THIPD 
INSERT CO*CGF 
OVERLAY THIRD 

INSER T ANTENA .DELAY, OILRN, OCFPLft .ICNREF.FCLE 
OVERLAY THIRD 
INSERT REARG 
OVERLAY THIRD 


^BBTOBUeiBILITY OF THE 
ORIGINAL PAGE IS POOR 
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. . » 

* * * 

August 11, 1973 

* 


INSERT TTPGRB .CEOOYN 
OVERLAY THIRD 

INSERT GEGIPH. ,'AVGPOT 
OVERLAY SECOND 
INSERT UPCATE 
OVERLAY EXTRA<RCGIONJ 
INSEPT ERROR 
OVERLAY EXTRA 
INSERT OnOl 
OVEFLAY EXTRA 
INSERT HEADER 
-OVERLAY EXTRA 
INSERT BMTfcRT 


REPRODUCI&ILITY OF THE 
ORIGINAL PAGE IS POOR 
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2.2 Hardware and Software Restrictions 

Normal operation of GEODYN requires a large scale 
TBM 360 Computer with 500K bytes of user accessable core, 
one 2314 direct access disk unit, four 9 track tape drives, 
one 7 track tape drive, one high speed ccrd reader, and one 
high speed printer. For some applications, GEODYN can 
operate with only one 9 track tape drive and no 7 track 
tape drive. 

The current GEODYN program is executable under 
versions 14, 16, and, 18 of the IBM 360 operating system. 

The compilation of GEODYN requires an IBM FORTRAN IV 
-'Level G compiler and for more efficient operation also 
an IBM FORTRAN IV Level H compiler 'with level 2 optimization. 

GEODYN contains four software restrictions. These 
restrictions and the subroutines affected are: 

1. Time difference tables must be kept up to date 
in subroutine TDIF. 

2 . Polar motion data tables should be kept up to 
date in subroutine POLE. 

♦ 

3 . Solar and Geomagnetic flux tables should be 
kept up to date In subroutines FLUXM and FLUXS. 

4 . Tables of the right ascension of Greenwich in 
blockdata subroutine CONSTS must be kept up to 
date . 

There are no GEODYN hardware restrictions. 

REPRODUCIBILITY OF THE 
■ ORIGINAL PAGE IS POOR 


2.0-11 



. f 

0 

3.0 GEODYN JOB SUBMITTAL ' ' 

i 

-*To prepare a GEODYN deck for job submittal the following 

..procedure is recommended: 

1. Define* the problem to be solved. 

2 « Prepare a deck of the necessary GEODYN Input 
Cards . 

3., Prepare the JCL necessary for the desired GEODYN 
execution. 

4 * Assemble the deck and submit to the computer. 

t 

3^ 1 Preparing the necessary GEODYN Input Cards . 

1. Determine those parameters of the problem 
definition which must be defined by COMMON 
SET cards (i . e ., gravity model modification , 
lunar and/or solar gravity modification, geo- 
potential coefficient recovery, station estimation 
parameters ,* modifications to the atmospheric 
density model, and modification to earth constants.] 

2. Determine which COMMON SET cards are necessary 
and prepare these cards. 

3. Determine how many arcs will be required and what 
parameters must be defined by each arc. 

4 . Determine which ARC SET cards are necessary 
for each arc and prepare these cards. 

5. Order the cards prepared in steps 2 and 4 so that 
they conform to the four groups within each 

set and place these sets in the proper order 
as follows: 
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COMMON SET - Run description (Group 1) 

COMMON SET - Option caprds (Group 2) 

COMMON SET - Station position cards (Group 3) 

COMMON SET - Termination card (Group 4) 

ARC SET - Arc description (Group 1} for ARC 1 

ARC SET L Option cards (Group 2) for ARC 1 

ARC SET - Termination card for Group 2 (Group 3) 
for ARC 1 

ARC SET - Selection/deletion cards (Group 4} for 
ARC 1 

« 

« 

ARC SET - Arc description (Group 1) for last arc 

ARC SET - Option cards (Group 2) for last arc 

. ARC SET - Termination card for Group 2 (Group 3) 
for last arc 

ARC SET - Selection/deletion cards (Group 4} for 
last arc 

NOTES: Always required for orbit generator mode is the ORBIT 
card. 

Usually required for orbit generator runs also is an 
OUTPUT card. 

STAEST, CORREL, BIAS, RESID, CULL, EDIT, PREPRO, SIGMA, 
VARCOV, and ARC SET, GROUP 4 - SELECTION/DELETION CARDS 
should not be used in orbit generator runs. 

An ORBIT card may not be used in a data reduction run 
' under any circumstances. 

A problem may consist of any combination of arcs in 
orbit generator mode and/or in data reduction mode. 

If any arcs are defined a :>eing in orbit generator 
mode, the maximum number of outer iterations must be 
one, and no common irameter estimation may be requested. A 
common parameter e: .mation problem must define all arcs as 

being in data reduction mode* 
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Preparing the Necessary JCL Cards. 

♦ 

1* Determine which JCL cards are always necessary 
for a GEODYN run. 

2. Determine which -JCL cards are required for the 
•specific mode of operation for the GEODYN run. 

3* Determine what additional JCL cards are necessary 
for special output requests. 

4 . Combine all necessary JCL cards into their proper 
order as specified in section 2.0. 

♦ 

GEODYN Deck Assembly and Job Submittal. 

1. Place the GEODYN Input Card deck into the JCL 
deck after the DATA 5 DD card and follow with a 
delimeter card. 

2. Submit the job to the computer with the proper job 
identification slip. 



I 


4.0 GEODYN EXAMPLE JOBS 

• k 

On the following pages is a comprehensive set of 
~~GE0DYN example job setups. GEODYN Example Jobs 1, 2, and 
3 include the setup .examples for the GEORGE, DELTA, 0RB1 
Conversion, and GROUNDTRACK programs. By combining the 
execution of the GEODYN Analyses and Graphics Support 
Programs with GEODYN execution, the most common usage 
of the support programs is demonstrated. 
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4.1 Example One 

. * 

This sample run is a 4 arc station estimation run 
using station constraints. Included are two 2 day .GEOS-I 
arcs and two 2 day GEOS- II arcs. GRARR biases and timing 
errors are estimated. Following the GEODYN execution is 
a GROUNDTRACK execution for arc 4. 

This job was run on . the 360/91 and requires core 
•and time as follows: 

CORE CPU I/O 

GEODYN step 394K 6.23 0..96 

GROUNDTRACK step 376K 0.03ra 0.02m 
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CEOOTN EXAMPLE JXH3 ONE * '* 

•• - . v - :/• •• ‘ * ■ 

- _ * g , * • ♦ • . . . \ ' ■ f - 

./✓••• * JOG . « • • • 

//GEODYN EX Z C PGWs7.CTV*'f,DM»P* Gl ON=300>C 

✓/STE»Lin f-o osn=m 2 . zct vM.LnAtL i n. r i sp=smi? 

//FTOjrooi on DSN = M2 *Z CT VM .E3H; M ,D f SP= SHR tDCD=UUFNO= } 

//FT05FC01 tO I>0 NAME=DaTA5 

//FTOSF JOl DO SY5CUT=A , CC G=( REC FM=VI? A.L^ECL®! 37 , PLKS 1 ZE = 7265 ) « 

// S 3 ACE=< CYL. 20 ) 

//FT07F00 1 DO FYGPUT=n,DCS=( ^ECF«=FB.LRECL= 00 * 0 LKSlZe-e 00 ) * 

// s=Aca=< ;aoo, ioo) 

V/FT 0?F 00 1 0 0 SY5CUT=a ,CCD-( 5 EC FM-Vt? A *LRECL = 1 37 « HLKS I ZE-726S >, 

// S 3 ACC=(C YL *10) 

✓/FT 09 F JO 1 DO SYSt>UT=A » CC l?— ( RL'C F P= Vtf A , LRECL = 1 37 , DLKS 1 ZE=7265 ), 

✓✓ S 3 ACE = ( C YL *10) 

//FT10F001 DD UV IT-DI SK,OCO=( HECFM=Fp, LRECL»A6 *OLKSI ZCsiACO) » 

✓/ S^ACC- ( TRK, 3 ) 

✓✓FT11F001 00 UN IT=Dl SF , DCl)=< R= CFM= VPS , LRECti = € 6 » PLKS I ZE=7294 ) , 

✓✓ S»ACF.= ( CYL, 3 ) 

/✓ FT 1 2 F 0 0 1 0 D UM IT-231 A ,DCB=( ?.LK£I7I:=7232 ,QU-^0 = 1 ) .SPACE “(CYL, 20) 

//FT 13r001 DO UN IT “1>1 SK,DCn = ( RE C F V “ VC , ML K SI ZE ~ 1 0-32 , GUi"MO = l J , 

✓/ S‘-»ACF.= ( C YL, 1 0 ) 

//FT11F001 DO UN I T “01 f.K , CCG- ( PE CF V “VS , OLKSI 7.E -l 5 O A , C*UFNO“ 1 ) * 

// S* ACC-C CYL ,5) 

//FT 15 r 00 1 OD UN I T = ? 3 l /« , DCG = ( 3L K SI Z C = 1 GO 0 ,OUF MR = 1 ) ,5-- PACE a ( CYL, 2 ) 

// r i i 5 r 0 0 i DD U h I 7 = 2 5 l A , uC H- { r«£C? r-VcS » LAE C. L = i dO-i , i:L L$ i 7 . •; - i 6 un , 

// RJFMU=2 ) .SPACt =(C YL. £ > 

//FT 13 F 00 1 DO UN IT = 01 SK, DC.5=( t’ E C “ '*=FH , L9 ECL = 3 0 • OLKSI 7E = 1 POO # OUFt.’O^I ) , 
|/^ S' 1 * ACE = ( C YL , 1 ) , P 1 5°=( .PASS) , OSMaGGRNDTPK . ^ 

//G0.FT20FC0 1 OD UN J T “ 2 4 0 0 - 5 , L A fl f L = .( ,9 ) * VOL = 5E R =OAT A1 , 

// DC(<= ( F E CfX = ‘ n , LPF CL = £0 ♦ •PL*' SI 7E = 320 } O » DCN=3 J 

//GO* 17 ! 2 5‘- 00 1 DD U M I T = 2 4 0 0— 5 , L A f* EL “ ( • -TUP ) * VOL =5 V. P = DAT A 2 • 

// DCH= ( r *J CFv,-po , LAF.CLa 60 , ^LN SI ZL = 32 00 0 i DC M = 3 ) * 

//go.- T 2cr-co i ro un i r= zaoc-?, l*ofl = ( ,fl° ) , vol“CC«-oata 3 , 

// f>CH= (PCCKMfFJJS «LRECL=90.PLK5I ZE = >?000tDCN=3) 

//SYSJDUMP DD SYSOUTSA ,SPACE=(C YL. I ) 

//DAT AS DD * 


STATION ESTIMATION USING 
Z CEOS -I 2- DAY ARCS WITH 
SOLVING FD9 GPTARR UlASCS 
5TALST 10 37. 

STAEST 1126. 

STAFST 10A2. 

STAEST 10 21. 

STAPOS 

|nnntNI02I 332549.7672 
JFTMY5 1022 2f.32S3.lA00 
|QOM(*;i 102A-3 I 232 D. i!7t*6 
1SAT AG 102^- 3 J ;i‘j:1.7G36 ' 
JM03AV1D3D 3 51 AA7. t;937 
1 JDLUI '<! 10 3 1-2 553 I . a 3 v 7 
|NF VP. I U 32 A7AA2O.2660 
S C OL 1: 510 3) G A 52 1 7 . *:r> 5 1 
I Gfr. IK 103A 
IKNFTL I 0 35 
IftOStW 10 37 


AC 121 . t> 3 O u 
3960 
3 512 7.2/*>fi 
I0H0KL 1 0 35- 3537 32 . <,:U / 


2 GEOS— I 1 ?-OaY ARCS AMD 
STATION Cr.sSTF A I NTS AND 
AND TIMING ERRORS 

5 . , 

5 « 

5 . 

5 * 

-53,693 
-42.000 
130. 403 
713.290 
C 76. 254 
1540. «77 
45.000 
1 06. 367 
203. 162 
50.41 
P4 9 . 5 3 3 
949.073 


5. 

0 . 

5 . 

202544*. 61 30 
27 * F 4 , t 6 0 0 
1 36521 5.1400 
2651553.6617 
243 500.9161 
2 74226.2077 
3071 6 A 6. 1 A 0 0 
212 ^26 .7602 
2625919. 5072 
2051* 7.934 
277 741.1600 
141*571 4, 1.501 


50 , 

so; 

SO* 

so* 


1037 • 
1037 . 


REPRODUCIBILITY OF THE 

ORIGINAL PAGE IS POOR 
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1 


t2f* 32 -=.^'-7-: -* v 

3771671 5- Z*' .“6 5 

-643059).2S» J7363 

-637V. €7271 1 15C M C 

3769. <.*€1973 35*56 

619.-+ )J> 1 251* 273S3 

SAT 1 1.22 

i*»2. 5 


PFEPROll 

$ . 


SJCMA 

1* 2«1 



«• 3.0 

* 

SICVA 

2« 10** 


SI CM A 

3* 12« C 


61 AS 1126 0. 

1C. 

2. 

eiAs ii26 c. 

C.10 

9. 

CATA 120 

- 

t) 

SELECT 66C3Z6 

C6C32S 


^NCAU_ 


- 

APC 2 



cns*i OPTICAL 04 TA 

" 


NOV«»«pd 2i.l«t3 TWO 

CAY A°C 


651124 



6*1 124 Ci COC C .COC-7 5 

1 


—3 637 1 07 • *2<, 6 = 4^523 

-€454 7*3 3.23672 135 

277251.53456 *6173 

31C 1. *121655 31 ! 63:? 

-2291 .2 7511 65 19 7 32 

-6472.291 921 J71379 

SAT | 1.22 

172. 5 


S 1 G v A 

1. 2.0 


S I O'"’ A 

6 . 2.0 


DATA 120 


' 

SCLECT651 12* 

6*1126 

1 

ENCALL 

■ — ■ 


ARC 3 

* 


CFOS— I I GRA=R DATA 



FEBRUARY 2 5 * 1 56 3 TWO 

CAY ARC 


jfcCC22S 



66C225 5 

2 


-G. ?7 2 56?7£ 2 9^64 9 6 AO 

C7 0.2 266476633^567 9 70 

C7 -C.71i6JfCj JU7723CD 07 

-0* 1C 5 <.2e5 4 2A3 755 5T0 

C4 C. 6 4C 52624992 / 'C74 60 

C 4 O.^il 2> i7o J6C76775D 04 

SAT 1 1.23 

211. 631 


SIGMA 

2* 10. 


SIGMA 

3. 12.6 


BIAS 1126 0. 

10. 

2. 

eiAS 1126 0. 

C.10 

3* 

G IAS 11 23 0. 

C.OC £ 


EIAS 1122 0. 

0.00 5 


CATA 125 



S£L E CT € ? C 2 2 5 

6CC22 7 

2 S 

SELECT66C225 

66C22 7 

3 6 

ENCALL 


P^RQDUCIBIUTY OF THE 

APC A 

. 

ORIGINAL PAGE IS POOR 

CEOS-I I OPTICAL CATA 



MAY 22.196c TwJ CAY 

ARC 

. 

CSC 5 22 



ESC 5 22 5 

2 


.421721 7. 71C51303A 

20471 63.062645225 

—6063 4 J 3 .947 2720?e 

• 3226.91 2 3C7 7A676C 

-S1B7 .0 7 3? 57 4 2 50 24 

-370*. 26 5 94. IOS757 

SAT 1 1.23 

e 

ill. 631 

♦ 

• 


4 . 0-4 



• % • 


f 


SXCma t« 2*0 

S1CVA Am 3m*. 

OUTPUT l 

oq^P^ai 1 J«0 , 

CATA 126 * . *. 

STL"CT6C:5^2 6«fS2A I 

fNOAU. i 

// EXZC L IN*. ;J3 • «' ”G ICN • OC— ?CCK 

, //L1N< • ^ Y5L l N ^0 « 

iNCLUC'i LOAOLi t( ZC«LCCA< I 
INCLUt C L C a OL l ■?» ( 2 C*- J Z *P L ) 

INC LUC' LC ACUl i< ZC- G* TVP) 

■TNTPy *\AIn 

✓✓GO » e Tl 1F<H 1 CC CSUsCGKNGTqir *0 I SP= C CL3 •CiL £T E I 
✓✓GO .DATA 5 r D - 

JCrS^NlfZT 3=12 7. 27? e 277 74 1*16wC 'JAS.PSJ 

r*o 

PLOTS 

TIMT frCC522- 6PCS24. . 

LNDPLT 
CAT* 
l AST 






t 

4*2 Example Two 

r 

« 

-This sample run is a five and -one half day GEOS-II 
data reduction run which estimates 2 pairs of resonant 
geopotential coefficients and a solar radiation parameter 
Following the GEODYN execution is a GEORGE execution 
which computes ROSRAN range zero set biases and timing 
errors. , 


This job was run on the 360/91, and requires core 
and time as follows: 

CORE CPU I/O 

GEODYN step 338K 4.16m 0.83m 

GEORGE step 326K 0.11m 0.06m 




* BEPRODUCIBILITY 
original page is 


OF THE 
BOOR 
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Mdrch 1, 1975 


GEODYtJ EXAMPLE JCIfV TWO 


* V • ' • 

joh* * • 

//GCODYN EXEC PG«= ZCTVMGDN.‘q = GI ON=350K ' . 

//STEPL1H tO rSN=V2 . ZCTVM .IOAPL CO,D 1 SP=SHR 
//FT 01 FO 0 1 OD f:SM= *-’2 ♦ z CT VM . n PHEM ,D1 5P-SMP. *nCO=OUFNO= 1 
//FT05F00 1 DO t‘DN A Ml' =■ D A T A S 

DD SYSOUTnA * CCC = ( * tCFM= VPA« LPECL ”1 37 < RLKS IZF = 7 26G ) % 

// S»ACE=< CYL. 2 0 ) 

//FT07F00J 30 G YiiO UT= 0 « DC C = < R EC F M=F 0 *LP.ECL- / ?0 » BLKS I 2E=C 0 0 ) • 

// S^ACL* = t ?00, 1 00 ) - 

//FTO3F00! DO SYSOUTsA f OCO=< R£C FM=VOA,LRf;CL = l 37 , OLKS IZ C =:726S ) • 

// S*ACE = ( CYL. 10 ) " 

//FT C?F00 1 0 0 SYSOUT = A *DCE=( RECF.v.= VF?A.LRF.CL = l 37 » ELKS I ZC = 7 2 6 5 ) , 

// S^ACCMCYL. 10) 

//FT J0FO0 1 00 UNITsDI 5K* OCH-( PECFKsF B « L7 ECL=4 6, OLKSI ZE=460) « 

// S= , ACt = ( r PK, 3 J 

//FTJjroOl DO UN IT =01 SK ,DC3 = < R£CFM= VOS «LRECL-6C »OLKS ] 2- =7 294 ) , 

// S=»AC«:=< CYL, 3) 

//FTJ2F00 1 00 UN IT = 231 4 ,DCB=C 9LKSI ZC-723?,9«J-.N0=1 ) ,SP*CE=( CYL, 20) 
//FT13f 00| DO UNITED! SX,OCO=< Rrcrwrvs.rUKSi ZE-~I 63’ .[JU-NO-l ) , 

//. S«»ACC»CCYL* 10) 

//FT UFOOJ DO UNlT = DISi<,CtO=(R£CF^V5,»JLKSJ Zz.~l A 04 .QU^.'CJ^l 
// S D ACf = { CYL . 3) 

//FT I6F001 DO UN l T = 231 4 *DC9=( - LKSl 7Es 1 300 .OU^WOnl ) *S n #'CE = CCYL»2J 
//FT i5 e 00! CO UM !! -2J * 4 ,OC‘i^{ 2 = CFY-VPS i L' 3rr ' s <«^,»> **c t 7 e- , c nn 
// CUFNOs !}, SPACE =( CYL . 0 ) ‘ ~ ’■'* 

//FT I ? 0 0 i CD UNITED! SK , OC^ ( RfCFv = VES * LREC L = M 2 .9LKSI 7 z^l -04 . 0UFN3= I) . 
// SPACE-' t CYL . 1 J. D I SP = < , PASS) , DSN=f.C.-_ljiiGc 

//FT20F00 I DO UN IT =2400- 4 , DC3- ( RECF tt^FQ . LRECL-t 0., 3LKSI ZE r f;000, 0EN=3 ) * 

// LA DEL = { ( l.'L^ ) « V0L-SLP-G“0S0A 7 A 
//SYSUDUMP 00 SYSOUT-A ,SPACE=(: YL, 1 ) 

//daYag DO * - 

A S 1^2 DAY GCOS-II DATA REDUCTION RUN w I TM THE ESTIMATION 

OP 2 PAIRS OF PE SONANT COEFFICIENTS AND A SOLAR P AO 1 AT 1 ON \> A R AMF. T £ R 

FOLLOWED -3Y A GEOFGE run COMPUTING ROSMVN CRAr.p TIMING ERK0»S 


nncoEFj^i^ 
flECOC=14 13 
OATA 2 

DATA P.EOUCTION 

CF.OS-J 1 GRARP DATA <? 

FERPUAf-.Y 21* 19 63 5 J /2 DAY ARC 

650221 

6*0221 5 2 
10751 4 3. 1 4 J 9 04 07 0 
-2040* 13524 4 212024 
SIGMA 
Si gma 

sol pad 1 .6 

OUTPUT - I 
tSAT I J « 23 
PNC PRO 1 1 
DATA 1 • 

SELECT 61*0221 0000 
CM DAL- 

// f.ALC L 1 N< G‘J« REG I Cl‘1 • CO= 4 00K 
//L I N< • SY SL I U 00 * 


0*5 

0 . 5 


0..5 

0.5 


H. 

F 1 • 


-754 6 4 F. 3 *-03012 4 559 
2$ 9* 1 95 07*3 002 2 22 j 


2. 

3* 

• 1 

211.931 


6C0226 1200 


10 . 

12. 


-001 650 « 4 77? 05 1*5 51 
-692 3*057 12 701 0592 


sasss? wtbb 

* ^ 


“FP 


20 



INCLUDE l_r I Jt * 

Include ljaclim £CRC»TYa'i 

ENTP Y MAjh 

//CO •f rT : SrOO 1 00 L'Sr4=CG5;o=C^ • Cl «Ps< CLC OiL' 1£ ) 
//CO .DA TA7 DO • 

RANGE ROSrtAK 

CLOY 

pATA 

LAST 


reproducibility of the 
ORIGINAL PAGE IS POOR 


4 . 0-8 



4*3 


Example Three 


a) Arc 1 is in orbit generator mode and uses converged 
double precision elements to generate an 0RB1 tape, 

b) Arc 2 is in data reduction mode and uses the same 
starting elements as arc 1 and the same data period 
that was used to obtain the arc 1 elements. Arc 2 
uses a truncated gravity model to reduce the elements 
An 0RB1 tape is also generated by arc 2. 

*■ 

c) A DELTA execution follows the GEODYN execution and 
computes and plots orbit differences from the two 
ORB1 tapes. 


d) An execution of the 0RB1 Conversion program follows 

the DELTA execution. The 0RB1 tape written in 
arc 1 is converted to a 7 track, IBM 7094, single 
precision, binary 0RB1 tape. 


This job was 
time as follows: 


run on the 360/91 and requires core and 

CORE CPU I/O 

2.36m 0 . 46ir. 


0 . 04m 
0.29s 


GEODYN step 348K 
DELTA step 200K 
ORB1 CONVERSION step 154K 


0,03m 
0 .09m 



ceooyn EXAMPLE JOD three 


March 2 , 


. . . Reproducibility 

■ ORIGINAL PAGE B poOr . 


//• » « JO? • • * . i ' 

✓/GEOOYM EXEC o^-ZCTV'iCON^- SION=.1f 0< 

//STE^LIO L'D OSNs 01 .G I HMP , GEO r »T « ( r J l SP=S*1R 

/✓FTOiROl DO OSM=GI . Z CT VM .{?=>ME M ,DI Sr»=!HP *DCB=U*UFNO= 1 
✓/FT05F00 1 OD f)ONA«E=OATA5 

✓/FT06F001 DO 5YSUUT=A .CCH = ( '<CCFM=VPA,L°.ECL = l 37. * ELKS t?E -7?£E*)» 

// ACC=( C YL . 20 ) 

✓/FT07F001 DO SYSO UT = 9 * DC B = ( 5ECF M=FU «LRECL= <30 » OLKSI Z E = SO 0 > • 

// S*ACr.= < 3 00, 100 J 

// FT 03 F 00 l OD SYSOUT- A « CCO=( RFC F M^V:TA »LRECL=l 37 « PLKSl'7.E=72G5 ) * 

✓/ S 3 ACC = ( C YL » l 0 ) 

✓/FT07P00 ! OD $YSUUT=A ,CCC--( R EC F«=V* A.LRECL=l 37« ELKS 1ZE = 726S). 

✓/ S D ACE = < C YL ♦ 10) 

✓/FTJO^UDl OD UN IT = DI SKt DCG'={ PECFM=FC)» L3 ECL=4 0 . F1LK 51 7£ =4 60 > , 

✓/ S^ACCst TP/ . 3) # 

✓/FT UFO 01 D D UN IT = 01 SK,0CI3=( PECFM=VES . LPEC L*=G H • ELKS I ?E =7 25 A >t 
// S°ACC-M CYL, 3 ) 

✓/FT 12 FO 01 00 UN lT = ?.:il A ,DCO- t PLKS1 7F*7Z3?VHJF»:0= I > .SPACE = (CYL. 20) 
✓/FY13F001 00 UN IT=-D! SK,DC3 = ( ''ECF«=VS,*1LKSI 2E-1632 .iiUFNO=l ) . 

✓/ S 3 ACF=( CYL, 10 ) 

✓/FTliFOOl DO >JN I T=D1 SK , CCM = ( RECP**-V$,P.LKS1 ZC = i 004 .PUFNO= 1 ) . 

✓/ S« ACE=l CYL, 5 ) 1 

tfr t | a i- o C 1 ijr*i » t — o 3 < * , CC^ s ( *.'1, ^ S ! Z l'" ! r-30* 1 ? u~ wris * J jSP ft CF = # C v L?2) 

✓/FTlSf^Cl DO UN 1 T ^23 1 A » DC {;=( PZCrr4= VUS » LR:.CL* 1 00 A « PLKS! 2E^ 1 £• OS . 

✓/ IV>ruO= 1 ) ,5PACK=(C YL • S ) 

✓/FT20F001 DO UN I T =• 2 A 0 0- A , DC 3 = C R = C - MsKI, LPECL“FO . OLK SI Z r *5C0 0 . 0FN = 3) » 

✓/ LA3EL=< . P.LP) . VOL=SER -GEDSUAT A __ 

✓ /FT 22F00 l 0 0 UN IT = 01 SK .DC <1 = 1 R E C F ** = VS . EL K S! 7E =2 00 A ,DUf NO-l ) • 

✓/ S»ACC=( 2404.50 ).0I SP=( . 3 AS? ) ,OSU = f.O?«U A 

✓/FT 22F002 CD UN IT = D J SK . PC B= ( E C F V- VS .PLK 51 ZE = 23 0 A « DUFNO= 1 ) , 

✓/ s-»acu=( 2aoa,so ) iC i sp=c .pass) iDSN=so^nm 

✓/SYSUDMi’P DO 3YS0UTsA.SPACL ? =CYLiI) 

✓/DAT AO OD * 

A TVO ARC FOLLOWED RY DELTA AMD ONE! CONVERSION 

ARC 1 IS AM O'OGEN USING CGNVFPGED ELEMENTS 

ARC 2 IS A DATA DEDUCT ION OVER THE SAME PER l OD WITH A TRUNCATED GRAVITY * 

STAPUS 

DATA I 

ARC I OP 3 I T GE NEPAT3R V/ITM 0RH1 TAPE OUTPUT 
APRIL 2? » 15 66 ONE DAY. ARC 

GEOS-l SATILL1TE STARTING ELEMENTS ARE PROM A PREVIOUS TATA REDUCTION 
6604 2? 

060423 1 

0.370035 l 5 CSS 7 1 A026 07 
0 .223066201 26 Cl 32 1 0 04 
DRAG 2,3 

SAT 1.23 

ODD I 300 660429. 

ORIUT 660420. 

OUTPUT 1 SCO. 

DATA 

ARC 2 DATA DEDUCTION Mill O 3 " I T A PE OUTPUT 
A M PII. ?'T»lvot> ONE DAY Al'C 1RUNCATED GRAVITY Pi^EL 
6E0S-1 SATILLITE SAME STARTING ELI RENTS AS PfEC.r Cl MG APC 


- • 5 » 3564 *>745 75502? COG 
6.71 544 1 5332035761 f.04 


- * 6* 5 5 M 1 4 1 2 5 3 7 0 4 1 
0.573161 2632255755 


172.5 


GAO 4 25. 
660425. 


4 . 0-10 


v* 



e* c* 2 s 

66C A ? 3 
C.37C 
0.223 
P4RTH 1616 
DRAG 
SAT 

pbfpp n l 

SIGMA 
SI GW a 
SIGMA 

orbi 

CATA 


;; 1 5cf3E7l AC2GC7 
r ( i3'.U , 6*'}l Zcl CCA 

— • 
0 4 

1 

2.3 

1.23 

172 

l.C 

1. 


2. 

C4 24 . 



c* ■»* r.4 .‘.7 j'j ?f i * < - tO'i 
7 J c 44 1 *Z i^> 337 cl CCA 


2 * t 
2.0 
o.c 

660 «>?• 


111 #l2tJ70*ltC7 
,>.*/ 31 ii* J322S> 7*? fc&62 


6 € 0 * 2 S 


.3^0 6GCA2B. 

1 

SetrCT663A2 3 
ENDALL 

// EXEC L IMKGC *Ps GtCN.GC=2ECK 
//l_IN< .SYSL IN PD * 

INCLUOF L0ADLli3( ZCTWDCLJ 
tNCLUOc Lr 4DLl?< ZCRG.IY® ) 

✓/goIft^'m do ccB=<K£CFF=Fe, t-Feci.*& .ilks xiE=*csz.oi «=i «. 

// LAP.»t. = I .ei.ol.UMIT = 24f C-7. VDL = SE e «SC-’TCH 
//G O^^TriFCZl CD O SM— CCFD l A #D IS F“ l CL>- * ASS) 

V/GC •• f T« 2 FC 0 1 CD DSNaEOMSl^.OlSPaCCCC. OILET e» 

✓/GO .DUM DC DSNs*,|_IMK.SYSLl'OC*.OISP*COLJ.O* l-TE > 

✓✓G0.PATA5 DO * 

T 61 

COMPARING A CONVERGED l DAY A PC 
VI TH TH‘J SAM£ ARC CONVERGED 
W * TM A Ti{«j.NC * Tf.O TAOTH V.-JCFL 
// <“ MCC tnt T?A-j 
/✓SOURCE • SYS IN DD * 

THE OK5J CONVERSION p C"TP AN DECK C« = S hERE 


// EXEC LINkGO _ 

//GO.PTI 3FC Cl DO CbNsfiCROl* lSF = (OLDf!)-.UT E) <-., L _, l0C 

//G0.PT11FCM CD UNlT S ?UC-7.CCE=(-:c=v^3.tUCL-21.J-<.Ut-lO-' 
✓/ DEN- I >,LA3EL=< »BLP > , VCL=S EG =SCFTCH 
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4.4 Example Four 

. * 

This sample run is a 12 day OGO-3 data reduction run. 
which estimates a drag parameter. Both data and elements 
for this run were obtained from the DODS Data Base during 
the GEODYN execution. GEODYN chose to use a variable step 

size for the orbit integration. The short amount of time 

required for this run should be noted. 

This job was run on the 360/95 and requires core ' 
and time as follows: 

- * 

Core = 364K 

CPU * 1 . 5 min 

I/O a 0.2 min 

This is a twelve day arc for which 4 iterations were 
performed for a total of 48 iteration days. The time required 
per iteration- day was 

2 seconds/iteration-day 

which is approximately five times faster than would 
have resulted from using a fixed 100 second step size for 

the integrator instead of variable step integration as GEODYN 
chose. 


original . 
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GCODYN EXAMPLE. 305 FOUR 


//. . • 4*)^» • • * , * 

//GCODYM EXEC PGM=ZCTV*CON.REC! 0N=3F0K • . 

//STF’LIB DO Df»M=Gl .GIPHP.GrOSTR t OISP = SHR * 

//FTOirOGj DO DSN=Gt • ZCTVM.EPMt M,DI S°sS^P ,DCH=5UFN0= 1 
//ftopfooi ivj dsn=po?js « um j t=23i a * voL*sERroooso 2 *c i c .p=shr • 

// S 5 ACE=< rw* (3304 ) > »CCf?=( a FC FM=F • L’L^ 312; = 3520 « MU” N'0= 1 ) 

//FT 03** 0.0] DO UN IT = 2*0 0-9 .CC5=< R EC- Ms YDS .LPECL= 1 O A .5L*S I ZE» 1 O'* A ) • 
// LA MEL - < « HLP ) • VOL — S!;R = CUT0 A T A 

//rr 041-0 01 DO DS\'=P5LS. UNI T = 2 .*J 1 a* VOL = GF.Ri:OOOS 0 2 .DIS^=CHP* 

// 5’ ACE = i T PK « C2Z1 ) ) , UC0=(REC-M=F, PLKSI ZE= 3520 , C*UFNO= 1 I 

/✓FT 05 F 001 DO DOUAML-OATA5 

//FT 06 F 00 l DO SY5QUT= A ,OCG = ( 5 IZ F H~ VH A , LOCCL = 1 37 » RLKS I ZE = 7265 } * 

✓/ S’ACE = < ZYL* 20} 

//FT07F00I U'J Srr. 0 UT = 3,CCD = { RECF P=FQ ,LRLCL= 00 , DL KS l Z 5 =Q 0 0 ) • 

// S’ ACC=< 300* 100 ) 

//FT03FO0 1 DO 5YSDUT“A ,CCE=( PEC FHr VBA , LPECL=1 27 . PLKS JZZ =7265)* 

✓/ S’ACF.= < CYL. 1 0) * 

//FT 05 F 00 1 DD S Y50U1 =A * OCE = < RECFM=VUA*LRECL = l 37 « 3LKS f ZE = 7265 ) • 

// s’/.ce={ :yl. io) * 

//FT to TOO 1 !>D UN IT = D! SK*DCI? = C nECF.v^FB. LRECL^A 6 • BLKSI 7.E = 4 60> * 

// S 3 ACCa(TflK,3) 

//FT 1 1 FOO 1 DO UNIT-DISK* OC B=( PE CFM-'VBS . LFt CL =60 * PLUS 1 Zc =729* )« 

/A S’ACF.= < CYL, 3) 

. > i t 4c* ■/«/ < u i 4 i - ^3 1 A iDwO-i DI. Iwl (.1. — • 4.D7 |D»J • *3 — I}«CFACC — *CTL.2C# 

//rTl j?00l Dt) UN IT -U) S K . D C b - ( I ' r. C r S ^ L r> i ZK. = i Oj2 tHUrNo- 1 ) * 

// S’ ACE=( CYL* 1 0 ) 

//FT 14 TOO 1 DO UN IT = D ISK,DC9=( P £ C r Ha V S . OLKSI ZE=l SO* *BU-NO = l ) « 

// S’ACL=( CYL. 5) 

//FT 1 5 F 0 0 1 DO UNIT.-231 * ,DC'i=( 3L K S I ZK = 1 ”0 0 «MU T - Kfl= l ) *^.P AC F = I C Y L. 2) 
//rr 15FOOI DO UNIT 231 A.DCD={ R;CFMaVUS*L RE C L = l C 0 4 « OLKSI Z(T = 1 £ O’T * 

// OJFNO= 1 ) ,SPACi:= (CYL . 5 ) 

//SYSUD'JVJ 00 SYSOUT-A , SPACE = ( C YL. 1 > 


ION fDR THE OGC 3 SATELLITE 
EVENTS FROM THE ODDS DATA -=$ASE 
'<£ DATA USED ONTO A TAPE ON UNIT 3 

! 3 T A I f 1 1 NG DATA FROM DOOS DATA BASS 
' A DPAC* COEFFICIENT 

AULT FOR INTEGRATOR STEP SIZE AND MODE 


REPRODUCIBIUtV of the 
ORIGINAL PAGE IS poor 


'70051 0 

'* 1000 . 

.5 

212 . 5 


r 


4.5 Example Five 


This sample run is a one arc geopotential recovery 
^run with B-matrix output. Six coefficients are being 


estimated as follows : 


Coefficient 


a priori 
value 


starting 

value 


standard 

deviation 


C(22,14) 
S(22 , 14) 


-9.25362D-27 -2.51579D-26 II of a priori value 

1.12578D-25 8.25801D-26 II of a priori value 


C(21 ,14) 
S(21 , 14) 
C(20,14) 
S(20, 14) 


SAO 1969 
SAO 1969 
SAO 1969 
SAO 1969 


SAO 

1969 

11 

SAO 

1969 

1% 

SAO 

1969 

1% 

SAO 

1969 

11 


of a priori value 
of a priori value 
of a priori value 
of a priori value 


This job was run on the 360/91 and requires core and 
time as follows: 


Core 

356K 

CPU 

3.1 m 

I/O = 

0.97m 
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CEODVN EXAMPLE 303 FIVE 


//. * • JO»). • 

//GEODYM EXEC PGMsf ZCTVMGON.Pr 91 0N=3S0< 

//STE-*I. I 3 D?> PSN- M? .ZCTVM. LfK. DL 1 Pi D I Sn=SHP 

//FT 01 F00 1 OD DSM= M2 » 2 CT VH . E 3 Hr M ( O I 5P a SSP »0C B-OUF NO— 1 

/✓FT 03 1’OO 1 00 DUKAME&DATA5 

//FT OL»rOOl DD SY5UUT=A tCCG-( ^rCFM=VC J .A«UntCL=137 , RLKS I7E=7265 )• 

✓/ S 3 AC E = ( C YU *20) 

//FT 07 F 001 DO GYf.OUT=». pCl 3 --( * E.CF «~F 0 ,LPECL= 30 « CL KSI 7 E = BOO I • 

// S»ACE=< ft 00, 1 OO ) 

//FTOSFOOX DO SYSOUT=A ,OCU=( PECF P.»VOA *LOECL-l 37 * HLKS I ZF-7265 ) • 

// S 3 ACE-I CYL. 10 ) * * • 

//FT O? FOOT 00 SYSCUT=A *CCU = < 3 EC F M= V 0 A t LRCCL= 1 37 . FLKS I ZF. = 7?d5 ) • 

// S*»ACE=< CYL. ! 0) 

// FT 1 0 F 0 0 l DO UN I T = D I SK* OCO^C cc C FHsK C , UP ECL =4 6 * R LKS1 Zc = 460) » 

✓/ S»ACE=( Tf'K. 3> 

//FT 11F00 1 00 UN 1T=D1 SK,OCR=( P = CFK= VOS «LPCCL,=68 , CLKS IZE=7294 >• 

// S 1 * ACE- ( CYL * 3 ) 

//FT 12 F 0 0 1 DD U» I T = 23 l A »DCK=C PLKSI ZF»7 232, *3UF?;n=: 1 ) . SPACE » C CYL* 2 0 ) 

✓/FT I3F00I OD U»|T±OlSK«0CO=(?rCr P»VS«'JUKSI 7E=10 32 «nU."MO=l J . 

// S»ACF=< CYL, 10) 

//FT 14100 1 OD UN 1 T=DI SK, CCD~( PECF VS • ELKSI ZE = 1 204 ,nu e NO=l ) • 

// $>ACE~( CYL, 5 > 

//FT 1 3FG0 1 OD UN IT - 23 1 iDCO-f 3LK.S1 7F« I 30 0. *JUFKO = 1 1 , 5*»*CF * ( CY L, 2 ) 

/* I t # / • VV I -/ fc' U 1 I * - ^ f 1 I L -l “ , v • — ■ —,V « u«'«- V L " 4 ^ vl 1 4 .U A.U 4 *-L “ 4 •- V - I 

// KRIf NO at 1 ) , 5S J AC i; = < C f'L . 5 1 

//FT2DF 301 OD UN I T = 2 4 0 0- 4 . OCR = < PCC“M=-*} « LPS CL =6*0 » 5LKSI ?. E = 5C 0 0 . D£ N - 3 ) . 

// LAn*.L = ( .ML* 1 ) . VOLsStRsO-CSSAT A 

//FT71 *-001 OD UN IT*240C- A ,DC-»-< P EC s «■ VUG . L^ECL-OOOn , ELKS I 7E=50 1.2, 

// D=U~3) , «AMEL= ( 1 . UL P ) ♦ VOL*?? R=P.fr*ATR X 3- * 

r — — EEPRODUCIBILITY OF THE 
*** 4 *** 9 ** ORIGINAL PAGE IS POOR 

DATA 5 roue T IO.N WITH ^MATRIX OUTPUT 


PE CO El- 2 2 1 A 
P^COC* 7 22 1 A 
RECOEF 21 1 A 
RcCOC- 20 I A 
STAPOS 
13PO IN 10 21 
|tr T MYT 10 22 
IQOUiZIM 024- 
I SAT AO 10 2:3- 
SMOJAV 1030 
IJD'UJT 1031 - 
INEWFL 1032 

lCour.ni 033 
10.Fpr.< 103/, 

I ifNKF. 10 TO 

fOOGPNt 037 
IJPCP*. 1 0 31?- 
|POSW > 1042 
1 T AM AN 10 '* 1 - 
DATA 
ARC 1 


-9 .25 3620* 
-2.31 5790* 


27 

20 


1 . 1 2S7F0* 

e. 2 sco 10 * 


25 

26 


01 


01 


332549. 
263253. 
312 J2ft. 
3 3 « 5 3 . 

3 51 v4 7 . 

2 053 1 . 

4 744 2°. 

64 0217 . 
4 15 121. 

512646. 

3 012 7. 
103732. 
3512 7 . 
19 0 32. 


767 2 
1400 
l: 760 
7 9 36 
»9 37 
4 397 
2060 
*;05l 
5300 
3~6 0 
2 793 
ofil 7 
2 94 3 

5125 


26254 4 fi, 

276 6 a; 
130521 0. 
2P9 l 953. 
243 552. 

274226. 
3 07.1 64 6. 
212 936. 
2 625 c 1 9 . 
30912 7. 

277 74 1 « 
1 4 ■>’>>' 1 4 , 
277 74 0. 

471759 . 


61 30 
1 6 0 C 
I 4 00 
66 1 7 
9161 
2077 
14 00 
7 652 
0072 
9 34 
1600 
150 1 
56 0 0 
2C 79 


• 01 
• 01 

\ 

-53. 693 
-42.000 
130.4 03 
713.690 
S76. 254 
1 540.977 
4/i. 000 
156.367 
203. 162 
90.4 1 
C49. 933 
94 9. 57 3 
350.063 
1356.777 


• 0 1 
• 01 


1 « 
1 . 


f 
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GF.O.-l GPAt-n AND OPTICAL 
— WAPCM TVsO C.A 

efCJ2*3 

*d V 3 20‘^c?c JO 3 t 

121-502 7*5 >7 > 

J77. 17371 l Iiw*l C 
jfMATP< J • • 

SAT 1 1*22 

PHtPR] 11 

SIG^A 

SIGMA 

SIGMA 

SIGMA 

CATA 120 
Sfc'LKCT 6 6032 6 


i 

« * 
P 

CATA 

Y AK Z 


3771672 

mC 5Z 2 u-7 16 6 

372 5. 

£1*7 J35 9>56 

1* 

172* !5 
2* 


1« 

3.0 

C. 

2*1 

2 • 

10. 

3*' 

12. C 

■ 

* 

6<032 £ 



CNDALL 


« 


-±*335*3. J7 3«i3 
felS.OJiiilci 7 * 5 3 



REPRODUCIBILITY of THV 
ORIGINAL PAGE IS POOR 
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4.6 , Example Six' * 

This sample run is a 2-day orbit generator outputing 


a PCE data tape 

with Kepler elements. 

A 2-body orbit was 

used and the 

core and time required on 

the 360/91 are as 

follows : 


‘ 


Core 

•— 

308K 


CPU 

= 

4.65s 


I/O 

= 

0.36m 


t 

• 


I 


REPROD UCIBILITY OF THE 
ORIGINAL, PAGE IS POOR 
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GEODYN EXAMPLE JOB SIX 


* March 1, 1973 


* . t .... 

//« • • .JOB... • * * • ‘ 

✓✓GEODYN EXE C PCM=ZCTVMGON,R=GI ON?3flO< . ».*. 

✓/5TCPLIB CD D5M = M2 .ZCTVM.LD4DL 10 .O 1 S*’ = S»m 

// FT 01 FO01 DD L'Srj='«2-ZCTVM *DI SO = SHPOC«3=BUFNO= I 

✓✓FT02F00I 00 DUMMY 

✓✓FT 05FOO 1 0f> DDNAVE-OAT A 5 

/✓FT0SF00 I 0 0 SYSOUT-A «CC2=( * EC“M=VOAiLnECL -I 37 , BLKS I 7*=! = 72 CO I * 

✓/ S*ACE-(CYU.20> 

//FTO^FOOl 30 SYSOUT = U,DCO = ( R£C“M-FD,LPECL = 30 .flLKSI ZE=300 ) , 

✓✓ S 3 ACE=I -<00» 1 00 ) 

✓✓FTOSFOO I 00 3YS r JUT = 4 ,DCU = ( »ECFM= VM A. LRECLM 37 , PLUS 17E = 7 265 ) » 

✓✓ S 3 ACF=( C YL . 1 0 ) 

✓✓FT07FOO I DO SYSOUT = A ,DCfJ = C 7 ECF M= VO A , LPECL = i 37 . DLKS 1 7z 5 ) • . 

// £ 3 ACE = ( C YL *10) 

✓✓FTJOFOOl OD UN IT^OISK, DCH=( PECF«=F8* LRfiCL =4 6 . OLKSI 7E=460>« 

✓✓ . S 3 ACE-( TRK, 3 > , 

ft F T 1 JF001 00 UN IT = DI SK, DCP=< PZCFM=V9S .LPECL=CR ,PLKS JZE = 7 2?4 } , 

✓✓ S 3 ACE“( CYL. 3 > 

✓/FT 12 ^0 0 I 00 UN IT = 23 1 .OCff=C r»t_*< SI ZC=7 2.12 «9U"N0 = 1 > ,S*>ACr = (CYL, 2 0 > 
✓✓FT 13F00 1 00 UN IT -DI 5K, CC*J=< -iz. C F.v.s VS « OLKSI ZE = 1 632 *liUFNO=] ) . 

✓✓ S»ACE~(CYI. . 10) 

✓✓ FT I A FO 0 l 0 0 UN J T = DI SK . CCLl^C c = C FI-*- VS.DLKf.l ZE =1 HO-* ,'MJFnCi« I > « 

✓✓ S 3 ACC = ( C YL • f» ) 

f fr * « j i w w a *r ks u.i i i »tj« ■*h<Cl'~! i L K 31 ZC — T C 0 C t P> U' NO ! 1 *3 •’ACL — I C . Y L « 2 5 
ft r T 1 j i 00 1 Ol - ij'** I i = cj’i DC fi — ( ^ S i L L ~ i 5 0 » i'c. k; I Zl£ s= I £* 06 t 

✓✓ OJFNO= 1 ) .SPACE = ( CYL. S ) 

✓✓FT l 3 i*oo l 3 0 UV I T = D I SK . CC?=( ?•: C F M = VHS »LRtCL = l 00,OLK5IZE=;3520) • 

✓/ S 3 ACE - I C YL « .1 ) . D 1 SPs ( , PASS > . D SN- C S I MO A TA 
✓✓SYSUI>U« 3 OO SYSOUT=A , SPACE = ( C YL, 1 J 
✓✓DAT AS DO * 

TWO *ODY 01**11 T GENERATOR RUN 
PRODUCING A CART US I AM AMO A < E 3 LCR I AN 
WITH SELLER I AN PCL : DATA SIMULATION ONTO UNIT 17 
TIDES . * 

EARTH 2 2 

onDiesi 0. 

CODIES 2 0* 

DATA 

TWO CODY OPT IT 
NOTE THAT DRAG 
BY LEAVING OUT 
6 GO 7 22 

6607 22 | 

0* 1A 1)67 07 3 92 77 7 1 1580 
0 « A \ A 152 0 7 •“ 3 1 3 7 3 6 3 SO 


OMnr? mCm UTY 

ORIGINAL PAai l 


u y OF TFfP 

fJ ® B mo 1 ? 


♦ 1 


GEN ERA TOP RUN - TWO 
A NO 5 DL AR P AO l A T I ON 
T t IE SAT CARD 

# 


DAY APC 
PRESSURE 


ARE ZEROED 


ORB I T 
OUTPUT 1 
SI MOAT l 
OATA 


660722 . 
I 360.0 . 


07 -0,502&6S070 r i0074353 07 
CA 0*4 35 74 £436337 173 10 04 

660724, 


0 • 0765 372f> 3900 05520 07 
0.307 2991445* 307200 04 



S.O GEODYN ERROR. MESSAGES . - 

The GEODYN program prints descriptive diagnostic 
error messages for the most commonly occuring program 
input errors. Most of these error messages are completely 
self explanatory. 


5.1 DODS Data 

Base Errors * 

Observation Data 

a) 

f ■ . 

UNABLE TO READ DODS DATA BASE DATA 
DIRECTORY 

GEODYN EXECUTION TERMINATED 

b) 

NO DATA AVAILABLE FOR SATELLITE REQUESTED 

SATELLITE ID IS 

GEODYN EXECUTION TERMINATED 

c) 

NO DATA AVAILABLE FOR THE TIME PERIOD 

REQUESTED FOR SATELLITE 

GEODYN EXECUTION TERMINATED 

d) 

• 

DODS DATA BASE READ ERROR OCCURED FOR 

SATELLITE 

DATA RECORD NUMBER IS 


e ) GEODYN EXECUTION TERMINATED DUE TO 

EXCESSIVE DODS DATA EASE READ ERRORS 



f 


' * • 

. % 

• . 

Messages a, b, c, and e are all followed by a 
‘ . - -hexidecimal dump of the Data Base directory and of the 

Data Base data record in core. For interpretation of 
this information see 

f, DODS Data Base Element File Description, "GSFC 
• DODS Documentation 

2. Element Data 

a) UNABLE TO READ DODS DATA BASE ELEMENT 
DIRECTORY 

GEODYN EXECUTION TERMINATED 

b) NO ELEMENT SET IN DODS DATA BASE FOR 

SATELLITE 

GEODYN EXECUTION TERMINATED 

• 

C) DODS DATA BASE READ ERROR OCCURED IN 
ELEMENT FILE FOR 

SATELLITE 

DATA RECORD NUMBER IS 

GEODYN EXECUTION TERMINATED 
LABEL = -- 

d) DODS ELEMENT RECORD ERROR ENCOUNTERED 

FOR SATELLITE - 

ELEMENT TYPE INCORRECTLY SPECIFIED IN 

RECORD 

TYPE SPECIFIED IS -- 
GEODYN EXECUTION TERMINATED 
LABEL = 


\ 
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' e) UNABLE TO FIND ELEMENT SET 

: IN DODS DATA BASE ELEMENT FILE 

PROGRAM SEARCHING FOR ELEMENT SET 
CORRESPONDING TO EPOCH 

* 4 « 

1 r 

Messages a, b, c, and d are all followed by a 
hexidecimal dump of the Data Base directory and of the 
Data Base element record in core. For interpretation 
of this information see 

M DODS Data Base Element File Description, V GSFC 

DODS Documentation 

for read errors check JCL for Data Base files. If 
JCL is correct, see DODS Data Base maintenance personnel. 


5.2 Ephemeris Tape Errors 

a) PROGRAM TERMINATED 

INSUFFICIENT EPHEMERIS DATA 

•LAST DATA POINT YYMMDD HHMM SS.SSSSS 

. b) PROGRAM TERMINATED 

INSUFFICIENT EPHEMERIS DATA 

FIRST DATA POINT YYMMDD HHMM SS.SSSSS 

Both messages a and b indicate that Ephemeris data 
was requested out of range of the data available on the 
ephemeris tape. 


a) 

b) 


.Time requested later than available data. 
Time requested earlier than available data 



f 
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CL. 3 


—User Response: 

1) Check for raispunched, misread, or out-of- 
order GEODYN Input Cards 

. * M 

2) Check for ephemeris tape read errors. 

3) . Check for correct ephemeris tape JCL (including 

• correct tape number) . 


5.3 General GEODYN Diagnostic Error Messages 

a) STATION AT TIME YYMMDD HHMM SS.SSSS 

NOT IN TIME ORDER 

b) STATION NOT FOUND 

C) STATION NOT FOUND IN FILE 

d) OVER 200 DOPPLER PASSES - TABLE OVERFLOW 

e) PROGRAM TERMINATION DUE TO ILLEGAL INPUT 

EXPLANATION: INPUT CARTESIAN EPOCH ELEMENTS 

EQUIVALENT TO KEPLERIAN ELEMENTS WITH AN 
ECCENTRICITY GREATER THAN 1. PLEASE CHECK 
INPUT ELEMENTS AND EXAMINE INPUT DECKS FOR 
MISSING OR OUT-OF-ORDER CARDS. 

f) PROGRAM TERMINATION DUE TO INAPPROPRIATE INPUT 

EXPLANATION: ADJUSTED CARTESIAN EPOCH ELEMENTS 

EQUIVALENT TO KEPLERIAN ELEMENTS WITH ECCENTRI- 
CITY GREATER THAN 1. 


S. 0-4 
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PROBABLE CAUSES: . 

} • 1) POOR FIRST ITERATION ORBIT DUE TO POOR 

STARTING ELEMENTS OR A STEP SIZE TOO LARGE 

FOR THE NUMERICAL ORBIT INTEGRATION IN THE 

‘FIXED STEP MODE. 

* « 

or 

1) POOR FIRST ITERATION ORBIT DUE TO POOR 
STARTING ELEMENTS OR INAPPROPRIATE ERROR 
BOUNDS. ’ 

2) THE EXISTENCE OF WILD DATA POINTS WHICH 
WERE NOT EDITED FROM THE SOLUTION. 

g) PROGRAM TERMINATION DUE TO ILLEGAL INPUT 

EXPLANATION: THE CARD IN THE OPTION CARD 

GROUP IS ILLEGAL. PLEASE CHECK INPUT FOR 
KEYPUNCH ERRORS OR MISSING OR MISPLACED CARDS. 


h) PROGRAM TERMINATION DUE TO INSUFFICIENT 
OBSERVATIONS 

EXPLANATION: THE SETUP AND OPERATION OF THIS 

RUN WERE SUCH THAT NO OBSERVATIONS WERE 
AVAILABLE IN THE TIME PERIOD SPECIFIED FOR THE 
SATELLITE REQUESTED. PLEASE CHECK YOUR SETUP 
DECK OR TAPE ASSIGNMENTS. 

i) PROGRAM TERMINATION DUE TO INSUFFICIENT WEIGHTED 
OBSERVATIONS 

EXPLANATION: INPUT ELEMENTS AND EDITING CRITERIA 

WERE SUCH THAT FEWER THAN 8 OBSERVATIONS WERE 
LEFT IN THE SOLUTION. 
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. j) PROGRAM TERMINATION DUE TO MISSING DATA 
CARD 

EXPLANATION: PROGRAM ENCOUNTERED END OF FILE ON 

• DATA5 BEFORE READING DATA CARD TERMINATING 
LAST ARC. PLEASE CHECK SETUP DECK FOR MISSING 
OR. MISPLACED DATA CARD OR MISPLACED END OF FILE. 

.k) ILLEGAL STATION POSITION INPUT 
■ EXPLANATION: MORE THAN ONE STATION POSITION CARD 

WAS INPUT FOR STATION 

PROGRAM ACTION: THE DUPLICATE OF THIS STATION 

HAS BEEN REMOVED FROM THE TRACKING COMPLEMENT. 
THE FIRST STATION ENCOUNTERED WITH THIS NUMBER 
WAS USED. ’ 

PROGRAM EXECUTION WILL BE CONTINUED. 

1) PROGRAM TERMINATION DUE TO INSUFFICIENT A PRIORI 
INFORMATION- 

EXPLANATION: THE ADJUSTMENT OF STATION HAS 

BEEN REQUESTED BUT NO A- PRIORI POSITION WAS 
AVAILABLE FOR THIS STATION. 


Errors a, b, c, d, and k are non-fatal. Errors 
e, f, g, h. i, j, and 1 are fatal errors. In all cases 
checks should be made to be sure that the setup deck has 
been prepared in such a manner that the problem description 
is correct. 

All of these errors are self-explanatory and in all 
cases except possibly a, b, c, and d should be corrected. 
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m) EXECUTION TERMINATING DUE TO READ ERROR ON 

RANDOM ACCESS FILE --- 

■ •> « 

ERROR BUFFER IS 

n) EXECUTION TERMINATING DUE TO INPUT ERROR 

EXPLANATION: THE SETUP OF THIS RUN IS SUCH THAT 

THE" ESTIMATION OF GEOPOTENTIAL COEFFICIENTS OF 
DEGREE (ORDER) GREATER THAN USED IN THE GEO- 
• POTENTIAL EXPANSION HAS BEEN REQUESTED. 


Both errors ra and n are fatal. and should be corrected. 


Message m will be caused by hardware errors or improper 
random access file JCL. 


For error n, check setup for errors. If a coefficient 
is to be adjusted it must be of degree and order less than 
or equal to the maximum degree and order used in the geo- 
potential expansion. - 

o) PROGRAM TERMINATION DUE TO ILLEGAL INPUT 

EXPLANATION: THE CARD IN THE OPTION CARD 

GROUP IS MISSING A CONTINUATION. PLEASE CHECK 
INPUT FOR MISSING OR MISPLACED CARDS. 

p) PROGRAM TERMINATION DUE TO ILLEGAL INPUT 

EXPLANATION: AN ILLEGAL SURFACE DENSITY INCREMENT ' 

SIZE WAS SPECIFIED. PLEASE . CHECK INPUT FOR 
KEYPUNCH ERRORS. 
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q) EXECUTION TERMINATING DUE TO INAPPROPRIATE INPUT 
'EXPLANATION: THE SETUP OF -THIS DECK IS SUCH THAT 

o _'_THE NUMBER OF SURFACE DENSITY CONSTRAINT EQUATIONS 
IS GREATER THAN OR EQUAL’ TO THE NUMBER OF ADJUSTED 
SURFACE DENSITIES.. 

T j ** ********$ NEGATIVE ARGUMENT TO DSQRT FOR 

ELEMENT $********** 

s j **********| NEGATIVE ARGUMENT TO DSQRT FOR 
£***** 

( DRAG , SOLRAD or DRGDOT may appear' in this error 
statement.) 

****** a a * *5 NEGATIVE ARGUMENT TO DSQRT FOR 
BIAS $****** 

(Errors r, s and t are not fatal errors.) 

. u) EXECUTION TERMINATING DUE TO INSUFFICIENT MAIN 
CORE STORAGE. 

THE USER SHOULD SPECIFY REGION = K PLUS 

ADDITIONAL STORAGE FOR ALL EXCESSIVELY LARGE 
I/O BUFFERS USED. 

V) WARNING: 

INPUT ON THE POSITION VECTOR CARD IN ARC 

INDICATES THAT KEPLER ELEMENTS MAY HAVE AN 
ECCENTRICITY GREATER THAN 1. 

w) DATE ON FLUX CARD OUT OF RANGE OF TABLES. 

CARD IGNORED. EXECUTION CONTINUING. 
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x j **********$ NEGATIVE argument to dsqrt for 
COMPONENT OF MASTER STATION 

■i *********** 


y) PERMISSIBLE VALUES OF IORDER ARE 5 THROUGH 15, 

VALUE PASSED WAS ***** run TERMINATED ***** 


z) sat' 

INITIAL STEP TOO LARGE 
RESTARTING WITH SEC STEP 


aa) $**********$ ACCEPTED TRANSIT TIME ERROR 

AFTER SIX ITERATIONS = SECONDS. GREATER 

THAN DAYS $**********$ 

bb) **********$ ELEVATION NEGATIVE. MAXIMUM 

CORRECTION FOR PARALLACTIC REFRACTION USED. 

£ ft* ft & ft 


CC) INTEGRATION STARTING SUMS NOT CONVERGED AFTER 
ITERATIONS. 


* EPS 
DIFF1 = 
DIFF2 = 


EXECUTION CONTINUING 

dd) $**UPD0V, f N**$ ACCEPTED TRANSIT TIME ERROR 
AFTER SIX ITERATIONS = 

SECONDS. GREATER THAN 

DAYS $**UPDOWN**$ 
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•APPENDIX A * 

DODS DATA BASE 

A 

On the 360/95 Multi-Arc GEODYN has the capability 
o£ retrieving orbital elements and DODS format data from 
the DODS Data Base. 

Additional job control language is required when- 
ever the DODS Data Base is used. Unit 02 is necessary 
whenever data is retrieved and units 4 and 02 are 
necessary whenever elements are retrieved. 

- i 

The DODS02 disk pack must be mounted whenever the 
DODS Data Base is used. 


» 


A.l DODS DATA BASE ELEMENT RETRIEVAL 

O 

If the element retrieval flag on the Epoch card 
for an arc is greater than zero GEODYN will retrieve the 
Orbital elements from the data base. The satellite ID 
of the desired satellite must always be specified. 

There are three ways of retrieving elements: 

1. By element set number. 

If the element set number is specified on the 
Epoch card’ GEODYN will search for the cor- 
responding element set. If that element 
set is not available GEODYN will use the third 
way of retrieving elements. 

2. By data. 

If the element set number is less than zero 
GEODYN will read the first observation from 
the data base observation file falling after 
epoch for the satellite specified. GEODYN 
will obtain the element set number specified 
on the observation record and then use the 
first way of retrieving elements. 

3. By epoch. 

If the element set number is zero GEODYN will 
search for an element set on or before epoch. 
If no such element set can be found GEODYN 
will search for an element set after epoch. 
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If no element set can be found for the specified 
satellite the GEODYN execution will be terminated with 
a hexidecimal dump of the element set directory and the 
element set physical record buffer. 

• If an element set not corresponding to epoch is 
found the epoch will be reset. 

If elements are retrieved from the POD S Data Base^ 

for an arc the position and velocity vector cards should. - 

not be in the deck for that arc . . 

A. 2 DODS DATA BASE OBSERVATION RETRIEVAL 

/ 

If unit 02 is specified as the data tape input 
unit and if column 23 on the Epoch card is left blank 
GEODYN will retrieve observation data from the DODS data 
base. The data selected will be all data, except those 
data types that may be deleted on the Data Type Deletion 
card, between and including the start and cutoff times 
for data as specified on the Epoch card. 

The satellite ID of the satellite for which obser- 
vations are desired must always be specified on the SAT 

card. 

If no data can be found for the satellite specified 
the GEODYN execution will be terminated with a hexidecimal 
dump of the observation data directory and the observation 
physical record buffer. 


A- 3 


I 


APPENDIX B 
DATA TAPE OUTPUT 

'1 


If during a data reduction run the following 

conditions are true: for an arc then a data tape* will 
be generated on FORTRAN logical unit 3. 

• Column 23 on the Epoch card is less 
than 3 . 

• Column 7 on the DATA card for the arc is 

/ a 

blank or zero, 

, Column 8 on the DATA card is 2 or 3, 

f 

Notes : Any odd number specified in column 8 of the 

DATA card for an individual arc will tell 
GEODYN that it must read one Data Type Deletion 
card or Data Selection cards following the DATA card. 

If a data tape is to be written the proper JCL 
must be specified for unit 3. Separate files on 
unit 3 are required for each arc generating a 
data tape. 

* Data tape - The format of the output data 

tape will be the same as the 
input data tape. 
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* APPENDIX C 

t 

'INPUT/OUTPUT FILE FORMATS 

Within this appendix are contained all input/ 
output file formats for the GEODYN System which are not 
elsewhere described in this volume. 


f 


C.l B-MATRIX TAPE FORMAT 

* ■ » 

GEODYN has the capability to. write the normal 
equations out on tape. The word format of the binary 
B-MATRIX tape is as. follows : 

Record Record 


No. 

Size (bytes) 

Description 

1 

60 


Header record. 

2 

4N+8 


Matrix labels. 

3 

8N+12 


Matrix data by rows starting 

• 


' 

• 

N+2 

8N+2 


with right hand side 

N+3 

60 


Parameter set identification. 

N+4 

4N+4 


Parameter labels 

N+5 

8N+4 


Parameter values 

N+6 

60 


Trailer record. 


Where N 

is the number of parameters. 

-• 

• 

Header Record 

Word 

Type 

Value 

Description 

1 

1 

10001 

Record Type 

2 

I 

1 to 99998 

B- Matrix number 

3 

I 

N 1000 

Number of parameters 

4 

I 

N + 1 

Number of parameters + 1 

5 

D.P. 


Total variance 

6 

D.P. 


Weighted variance 

7 

d;p. 

0.0D0 

Arc variance 

8 

i 

\ 


Number of weighted observations 
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Word 

Type 

Value 

Description 

9 


9 

I 

8 

Matrix type 

- 

10-12 

R 

blanks 

Matrix name 


* 

Length * 

60 bytes 






Matrix Labels 

Word 

lH>e 

Value 

i ■ 

Description 

1 

I 

10002 

Record Type 

2 

t 

I 

0 

Dummy 

t ■ 

3 to 
N+2 

I 


Parameter labels 


Length = 

4N+8 bytes 



* 

Matrix 

Data Records 

Word 

Type 

Value 

Description 

1 

I 

10003 

Record Type 

2 to • 
N+2 

D.P. 

• 

Elements of J th matrix row, 
starting with right hand side 


Length = 

8N+4 bytes 
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Parameter Set Identification Record 


Word 

Type 

Value 

Description 

1 

I 

10011 

Record Type 

2 

I 

1 to 99998 

B-matrix number 

3 

I 

* 

Dummy 

4 

X 

N 1000 

Number of parameters 

5 

I . 

• 

Dummy 

6 

X 

8 

Code for parameters 

7-9 

R 

Blanks 

■ ✓ r 

Matrix name 


Length 

= 60 bytes 

• 



Parameter Labels 

Word 

T XP e 

Value 

Description 

1 

I 

10012 

Record Type 

2 to 
N+l . 

I 


Parameter labels 


Length 

= 4N+4 bytes 





Parameter Values 

Word 

Type 

Value 

Description 

1 

I 

10013 

Record Type 

2 to 

D.P. 


Parameter values 

N+l 





Length * 

8N+4 bytes 
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Trailer Record 


Word Type Value 

1 I -19991 

2-14 I • f> 


Length = 60 bytes 



C. 2 BINARY RESIDUAL TAPE FORMAT * 

# 

.The Binary Residual tape is written by . GEODYN 
Specifically for input to the GEORGE graphics and 
analysis program. The word format for this tape is 
the following: ; 


WoTd 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
*17 
18 
19 


1221 

1 

I 

R 

DP 

DP 

DP 

R 

R 

R 

DP 

R 

R 

R 

R 

I 

-DP 

- DP . 
DP 
DP 


Description 

Date in YYMMDD . 

Hours and minutes in HHMM. 

Seconds . 

a 

Station name. 

Alphanumeric Measurement type. 

Value of first observation. 

First observation residual (0-C] . 

First observation residual ratio to sigma. 
Time derivative of first observation. 

Value of second observation. 

Second observation residual (0-C) . 

Second observation residual ratio to sigma 
Time derivative of second observation. 
Elevation angle in degrees. 

Station network indicator. 

Along track resid ual component 
Cross track residual component 
Argument of Longitude in radians 
Right Ascension of Greenwich in 
radians 


Optical 

Data 

Only 


Total = 108 Bytes. 


The last record contains blanks. 


C. 3 DODS 

DATA TAPE 

FORMAT 

9 


• 

HEADER 

RECORD 

Word 

Number 

• 

Bytes 

Word 

T yP e 

Description 

1 . 

8 

D.P . 

Header record flag 
(value is -1.0) 

2-6 

40 

D.P. 

Spares 

7 

4 

1 

Satellite ID 

8-19 

48 

F.P. 

Spares 


Total 100 


Logical Record Length 



» 


• 


DATA 

RECORDS 

V 

Word 

Number 

Bytes 

Word 

Type 

Description 

1 

8 

D.P. 

. Time o£ observation 

2 

8 

D.P . 

Station ID 

3 

8 

D.P . 

Observation 

4 

8 

D.P. 

Not Used 

5 

8 

, D.P. 

Ambiguity corrections 


(bits 57-60) and 
transponder channel** 
(bits 63-64) 


6 

8 

D.P. 

Observation correction 

7 

4 

I 

Satellite ID 

8 

4 . 

I 

Not Used 

9 

4 

F.P . 

Measurement standard 
deviation 

10 

4 

F.P . ^ 

| 

11 

4 

j 

F.P. 1 

S Not used 

12 

4 

F.P. j 


13 

4 

F.P. J 

1 
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f 

DATA RECORDS (Cont.) 


Word 

Number 

Bytes 

Word 

Type 

Description 

14 

4 

• 

F.P. 

. Timing error in D.U.T 

15 

4 

F.P. 

Time correction value 




in D.U.T. 

16 

2 

I * 

Not Used 

17 

2 

I 


18 

2 

I 

i 

DODS OBSERVATION TYPE 

. 



NUMBER 

19 

2 

I 

Not Used 

20 

2 

I 

Time tag indicator; 


bits 1*8, 2=A1, any 
other value = UTC; 


bits 9-16, 5=satellite 
time, any other value = 
ground time. 

21 21 Observation correction 

indicator 


22- 2 I Not Used 

23 2 1 Spare 


Total 100 = Logical Record Length 

No Trailing Record 


C-9 



NOTES: 


Ambiguity corrections; 

1 ** +18737031.3 meters 

2 » +4684257.8 meters 

3 e +936851.6 meters 

* * 


Transponder channel; 

1 ** A Channel 

2 * B Channel 

3 ■ C Channel 



t 


C.4 EPHEMERIS TAPE FORMAT , 

The ephemeris tape is a binary unformated tape con- 
taining geocentric lunar positions at half day intervals, 
heliocentric positions of the earth-moon barycenter and the 
planets, Venus, Mars, Jupiter, and Saturn at 4 day intervals, 
and the nutation in obliquity at half day intervals. These. 

* * ephemerides were obtained by precessing and nutating to true of 

date coordinates the values found on the JPL planetary ephemeris. 
All quantities on the tape are accompanied by second and 
fourth differences for use in a fifth order Everett inter- 
polation scheme. Each block of data corresponding to one 
record on the JPL ephemeris contains 8 days of data with 
the first set of quantities being the same as the last set 
-of quantities in the previous block of data. The time of 
the first set of quantities in each block of data is 
specified in the A.l time system. The tape word format 
is as follows: 

* Record Word Type Description 

1 1 I Date of first set of quantities in 

. YYMMDD . 

1 2 I Hours and minutes in HHMM. 

1 3 R Seconds. 

1 4-54 R Nutation quantities. 

•1 55-81 DP Earth-moon barycenter quantities. 

2 * 1-51 DP. Lunar quantities. 

3 1-51 DP Lunar quantities. 

4 1-51 DP Lunar quantities. 

5 1-51 DP Venus quantities. 

5 28-54 DP Mars quantities. 

* 

6 1-27 DP. Jupiter quantities. 

6 27-54 DP Saturn quantities. 


C-ll 



f 


Records 1, *5, and 6 each have lengths of 108 single 

words. 

* 

Records 2, 3, and 4 each have lengths of 102 single 

words . 


The nutation quantities consist of 17 sets of values 
for the nutation in obliquity and its second and fourth 
differences. The first set corresponds to the time on the 
block of data and each successive set is one half day later. 

The earth-moon barycenter quantities consist of the 
coordinates of the sun relative to the earth-moon barycenter 
in meters. Three sets of coordinates with their second 
and fourth differences are stored as follows with the first 
set corresponding to the time on the block of data and 
each successive set 4 days later. 

X ‘ D x> D i tY > D y> D y’ Z ’ D l’ V t 


The lunar quantities consist of the geocentric co- 
ordinates of the moon in meters. 17 sets, each one half day 
later than the previous appear as follows: 




C-12 



t 


C.S GEOS DATA TAPE FORMAT 


Columns Format Description 

1-6 .16 Satellite ID 

' f 

7 II Measurement Type - same as described 

under option card SIGMA except 4 - 
doppler data. Doppler data is 
changed after preprocessing to type 
3 for program handling. 


8 II Observation identifier. 

0 » Active 
0 f Passive 


12-13 Time Identifier. 

12 II . <5 Station Time 

• >5 Satellite Time 

13 II 1 * UT-1 

2 *= UT-2 

3 « UT-C 

4 * A.l 
Default A.l 


14-18 15 

14-17 


Station Number. 

Station number for GEOS-II satellite 


15-18 


Station number for GEOS- I satellite. 


19-34 

19-24 16 
25-28 14 
23-34 16 


Greenwich Mean Time of Observation. 
YYMMDD - Year, Month, Day 

HHMM - Hours, Minutes 

SS.SSSS- Seconds 
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COLUMNS 35-80 FOR MEASUREMENT 
TYPES 1, 6, 7, 8, 13 5 14 


Columns Format Description 

3S-53 Observations 

35-37 13 * . Hours (arc) - right ascension, 

Degrees - azimuth or X angle. 

12 Minutes (time)' - right ascension. 

Minutes (arc) - azimuth 
Decimal degrees - X angle. 

IS Seconds (time) - right ascension, 

Seconds (arc) - azimuth 

Al Sign of declination or Y angle. 

46-47 12 .Degrees - declination, elevation, 

or Y angle. 


38-39 

40-44 

45 


48-49 12 

50-53 14 

62-63 

62 II 

63 II 

72-74 13 

75-77 13 


Minutes (arc) - declination or elevation 
Decimal degrees - Y angle. 

Seconds (arc) - declination or elevation 

Equator designation 

1 * True; otherwise. Mean. 

2 = Equator of reference time. 

3 « Equator of date. 

All other values; Equator of 1950. 

Standard deviation . 

Seconds - right ascension or azim ,h. 
Degrees - X angle. 

Standard deviation. 

Seconds - declination or elevation. 

Degrees - Y angle. 



t 


Columns 

35-53 

35-50 

51-53 

64 

65-67 

68-70 

71-76 

75-77 

80 


COLUMNS 35-80 FOR MEASUREMENT 
TYPES 2, 3 r 64. 

* 

Format Description 

• Observation. 


116* Range in meters; range rate in meters 

per second; doppler in cycles per second. 


13 


II 


13 

13 

16 

13 


II 


Millimeters of range; millimeters per 
second of range rate and .001 cycles per 
second of doppler. 

Tropospheric refraction correction 
indicator. 

«=2 or 4 - correction not applied. 

=1 or 3 - correction applied. 

=9 - correction not applied. 

Measurement type - French 
^ Laser Range. 

Range standard deviation in -meters. 

Millimeters or range standard deviation. 

Refraction coefficient for French Laser Data. 

Standard deviation of range rate in milli- 
meters per second or doppler in .001 cycles 
per second. 

Transponder channel for GRARR data. 

*=0 - A channel 
=1 - B channel 
=2 “ C channel 
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* 

COLUMNS 35-80 FOR 
MINITRACK DATA 

•i 


Columns 

Format 

Description 

38-53 

’ 

Observations. 

38 

A1 

Sign of Cosine Alpha. 

39-45 

17 

Cosine Alpha (Implied decimal in 39). 

46 

A1 

Sign of Cosine Beta. 

47-53 

17 

Cosine Beta (Implied decimal in 47) . m 
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A Mod of NGSP (GEOS) Format for VLBI and Average Range Rate , 

Input to GEODYN 

% 

* 

Field Columns Format Description 

1. 1-6 16 ' Satellite ID* 

* 

1-2 Year of Launch 

. 64 = 1964 
65 = 1965 
etc. 

Order of Launch 

Component Identifier 

1 * a 

2 » b 
etc . 

2. 7-8 12 Type of measurement 

81 = time delay VLBI t 
" 82 = fringe rate VLBI 

83 = two way average range 

rate 

84 = three way average 

range rate 


3-5 

6 


3. 


9-11 


Blanks 


4, 


12-13 Time Identifier 

12 II 5 » Satellite time 

0 * Station received time 


*as per COSPAR numbering system 
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c oo 



i 


Field 


Columns Format Description ' 

9 •* 

13 II 1 * UT1 

2 * UT2 

3 * UTC 

4 * A.l 

All other values UTC 

14-17 14 Station Number 

For measurement type 81 and 82 
Reference Station Designated. 

For measurement types 83 and 84 
Receiving Station Designated. 

Network numbering systems are 


as follows: 



0C01 

0099 

USB 


0100 

- 0499 

Not used 

0500 

- 0599 

Deep 

Space (JPL) 

0600 

- 0999 

Not used 

1000 

- 1999 

NASA 

S TAD AN 


1000 - 

10S9 

Minitrack 


1100- - 

1199 

S-Band Radar 


1500 - 

1599 

Minitrack 


1600 - 

1699 

VHF Radar 

2000 

- 2999 

• U.S 

. Navy Doppler 

3000 

- 3999 

U.S 

. Air Force Opt: 

4000 

- 4999 

Interagency C-Band 
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Field 


6 . 

7. 


8. 


Columns Format Description 

5000 - 5999 U.S. Army Secor 

6000 - 6999 Worldwide BC-4 Network 

♦ 

o 

7000 - 7999 NASA SPEOPT 

_ 8000 - 8999 International Observers 

9000 - 9999 SAO 
18 Blank 

19-34 G MT of Observation 

19-20 12 Year 

64 *= 1964 

65 - 1965 
etc. 

21-22 12 Month 

23-24 12 Day 

25-26 12 Hour 

27-28 12 Minute 

29-30 12 Second 

31-34 .0001 Seconds 

35-36 Observation time correction 

35-36 12 Time correction in microseconds to 

be added to GMT. 


C-19 


» March. 1, 1973 


^Dcscript icn 

Refraction Indicator 
[ m *0 Obs. corrected for troposphere 
• and ionosphere 

.t 1 Obs. corrected for ionosphere 

but not troposphere 

‘* •2 Obs. corrected for troposphere 
* * 

-but not ionosphere 
*= 3 Obs. not corrected for tropo- 
sphere or ionosphere 

Transmitter Frequency Indicator 
« 0 Unified S-Band 
s 1 GRARR S-Band' Channel A 

* 2 GRARR S-Band Channel B 
*= 3 GRARR S-Band Channel C 

* 4 GRARR VHP 

Station Number 

For measurement types 81 and 82 
the second* receiving station 

For measurement types S3 and S4 . 
tlie transmitting station 

Same numbering convention as 
columns 14-17 , 

cIB Sca l e Factor Indicator for 
Obs erv at ions 

Measurements should be multiplied 
by 10< IE ^ to achieve nominal un 
For example : 

IE = 4 indicates R observations 
arc in nominal units. 

IE ~ 7 indicates IT observations 
arc in !Cn/scc. 

IE = 1 indicates k observations 



Field 

13. 


14. 


15. 


Columns Format 


4*4-58 


F15.10 


* >V. 


69-65 


F7, 6 


SO for If only 


66 


II 


67-69 


53.1 


70-72 


F3, 1 


Dcscrint ion 

Unsealed Observatio n 
(Decimal implied between 
columns .48 and 49 ) 

Nominal Units: 

(see column 43 description) 

T g - Km 

Vp^M/sec 

^ - M/scc ■ 

* 

Observation Siena 
(Decimal implied between 
columns 59 and 60) 

Units : 


Yf “ 
IT - 


Km 

M/scc 

M/scc 


Receiver Mount* 

I XY (Hast -West) 

^ (North -Sou th) 
Azimuth - Elevation 
Hour angle - Declinat 


2 

3 

4 


ion 


Ante nna Axis lUcmi 

• — - — <J iSMlar pPo^t £ 0 

Receiver in rZ * 

• b e t,.„ (Decimal inpii« d 

^lIHLAxiLjU fQr 

l^HIiller J.n McT^T T^TCT, 

i-Plicd between columns 71^ 72 ) 


l 


Field 


Columns 


Format 


Dcscj 5 ntion 



* 73 


II 



. . 74*80 17 


Transmitter Mount Type 
- » 1 XY (East-West) ' 

• » 2* XY (North-South) 

**3 Azimuth - Elevation 

• * 4 Hour angle - Declination 

» 1 • * 

• • 

Counting Interval for R in 

♦ -Milliseconds 




1 


A Mod of NGSP (GEOS) Format for Altimeter and Satellite to 
Satellite Data Input to GEODYN 


Field 

■ * 
Columns 

Format 

Description 

1. 

1-6 

16 

Satellite ID* 

* 

1-2 

• 

Year of Launch 

64 = 1964 

65 = 1965 
• etc. 


3-S 


Order of Launch 


6 


Component Identifier 

1 c a 

2 = b 
etc. 

2. 

7-8 

12 

Type of measurement 


92 = Summed Range ) 

93 = Summed Range / Sat 

Rate ) -sat 

94 = Altimeter Height 

and Height Rate 


3. 9-11 


Blanks 


4* 12-13 1-2 Time Identifier 


12 S « Satellite time 

0 * Station received time 
6 = Satellite time corrected 
for satellite clock 
errors 


*as per COSPAR numbering system 
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9 


K 


Field Columns Format 


Description 


13 


1 * UT1 

2 * UT2 

3 * UTC 

4 = A.l 

All other values UTC 


5, 14-17 14 Station Number 


For measurement type 94 
station is always -4. 

For measurement types 92 
and 93 Tracking Station 
Designated . 

Network -numbering systems 
are as follows: 


0001 - 0099 USB 

0100 - 0499 Not used 

0500 - 0599 Deep Space (JPL) 

0600 - 0999 Not used 


1000 - 1999 NASA STADAN 


1000 - 

1099 

Mini track 

1100 - 

1199 

S-3and Radar 

1500 - 

1599 

Mini track 

*1600 - 

1699 

VHF Radar 


2000 

3000 

4000 

5000 

6000 

7000 

8000 

9000 


2999 U.S. Navy Doppler 
5999 U.S. Air Force 
Optical 

4999 Interagency C-3and 
5999 U.S. Army Secor 
6999 Worldwide 304 
Network 

7999 NASA SPE0PT 
8999 International 
Observers 
9999 SAO 


6 


18 


Blank 



Field 

7. 


8 . 




Columns Format 


19-34 


Description 

GMT of Observation 


19-20 12 


Year 

64 = 1964 

65 * 1965 




• 

etc. 

21-22 

12 


Month 

23-24 

12 


Day 

25-26 

12 


Hour 

27-28 

12 


Minute 

29-30 

12 


Second 

31-34 

14 


10~^ Seconds 

35 

11 


Refraction Ind: 


. 0 “ Observation not cor- 
rected for troposphere 
or ionosphere. 

1 = Observation corrected 

troposphere but not fo 
ionosphere . 

2 = Observation corrected fo 

ionosphere but not for 
troposphere 

3 = Observation corrected fo 

both troposphere and 
ionosphere . 


36-54 Observation of Range, Range 

’ 4 Rate , or Height^ 




hi 


Field 

Columns 

Format 

Description 

(9.cont.) 

36-51 * 

116 

Meters of Range or Height, or 
millimeters/second of Range 

Rate. 

■ .* 

52-54 

13 

Millimeters of Range or Height, 
or micrometers/sec of Range Rate 


55-80 


For Range or Range Rate 

10. 

55-60 

16 

Sigma for Range in millimeters , 
or Range Rate in micrometers/ 
second.. 

11. 

61-66 

16 

Satellite ID for relay 
satellite. 

12. 

67 

.11 

Transponder Channel 

- 

- 


1 - A channel 

2 - B channel 

3 - C channel 


55-80 


For Altimeter 

10. 55 

-56 

12 

Meters of Height Sigma. 

11. 57 


11 

10 ^ meters of Height Sigma. 


12 . 


58-60 


13 


Centimeters/second of Height 
Rate Sigma. 



Field 

Columns 

Format 

Description 

13. 

61 

m 

11 

Subsatellite geography 
indicator 



• . 

0 « No information avail 

able. 

1 ** water (ocean or sea) 

2 * land (includes lakes 


14. 


62-80 Height Rate Observation 

62-77 116 Millimeters/second of Height 

Rate measurement. 


Micrometers/second of Height 
Rate measurement. 


78-80 


13 


w 



C.6 GEOS-C DATA FORMAT 


April 8, 1974 


INTRODUCTION 

The GEOS-C satellite will be tracked by many different 
systems and networks of tracking stations. It is intended 
that this format will accomodate data taken from all track- 
ing systems and networks as well as altimeter data taken 
by the satellite. As a consequence of its design character- 
istics this format not only meets the specifications of the 
GEOS-C project but also the- precise orbit determination 
requirements of all currently existing satellites. 

Since the format is a card image representation, it 
is ideally suited for transferring data between the different 
types of computer facilities supporting the GEOS-C project 
research. The only character images used by the format are: 

integer numbers 0-9 

letters of the alphabet A-Z 

blanks, and 

minus signs 

All decimal points are implied, therefore eliminating the 
need for that character representation. 
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Api jli «| i.y/'< \ 

Descriptions of columns 1-32 for all measurement types • 

Columns Subset Description 

1-7 Satellite -ID 

8-9 Measurement Type 

10-14 Right Ascension and Declination 
• 10 - Optical 

11 « Laser Optical 
15-19 Hour Angle and Declination 
15 = DSN S-Band Radar 
20-29 Range or Range Difference 

20 = Laser 

21 “ CrBand Radar 

“ 22 « C-Band VLBI Range Difference Radars 

23 15 DSN S-Band Radar 

24 = USB Radar .. 

25 = GRARR S-Band Radar 

26 «= GRARR ATS-R C-Band Radar 

27 = GRARR VHF Radar 

30-39 Range Rate or Range Rate Difference 

30 * Combined DSN/USB Radars 

31 = C-Band Radar 

32.= C-Band VLBI Range Rate Difference 
Radars 

33 = DSN S-Band Radar 

34 = USB Radar 

35 = GRARR S-Band Radar 

36 « GRARR ATS-R C-Band Radar 

37 = GRARR VHF Radar 

38 « TRANET Doppler 

39 83 Geoceiver Doppler 
40-44 Altimeter Height 

45-49 Available for additional measurement 
types such as Altimeter Height Rate. 

C- 26. 2 


April 8, 1974 


Columns 


10-11 


x 


Subset Description 

50-54 l and m Direction Cosines 

51 e Minitrack Equatorial Mode 

52 » Minitrack Polar Mode 

55-59 Available for additional angles 
or cosines 

60-63 X-Y Angles (East-West) 

60 « USB Radar 

64-69 X-Y Angles (North-South) 

64 «• USB Radar 

65 = GRARR S-Band Radar 

66 = GRARR ATS-R C-Band Radar 

67 = GRARR VHP Radar 

70-79 Azimuth and Elevation Angles 

70 = Laser 

71 ** C-Band Radar 

73 = DSN S-Band Radar 
80-89 Available for additional Range type 
measurements. 

90-99 Available for additional Range 
Rate type measurements. 


10 


11 


Time System Indicator 

0 = Ground Received Time 

1 = Satellite Transponder/Transmitter Time 

2 = Ground Transmitted Time 

3 = Satellite Receiver Time 

0 = UT-0 

1 *= UT-1 

2 - UT-2 

3 « UT-C 

4 « A . 1 

5 « A. 3 (A. T. B . I .II,) 

6 c A-S (Smithsonian) 
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April 8, l y/4 


Columns 

Subset 

Description 

12-16 


Station ID 


17-32 

17-18 

GMT o£ Observation 
year 
YY 

•w 


19-21 

Day of Year 



22-26 

DDD 

Time of Day (Seconds from 

midnight 


27-32 

GMT) 

Fractional Part of Seconds 

(in 



microseconds) 
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•* V 

Columns 33-80 for Right Ascension/Dcclination and Hour Angle/ 
Declination Data (Types 10-19) 

Columns Subset Description 


33 


34-35 


36-54 


55-57 




• Beacon activity indicator 
O' = Beacon Active 
1 *» Beacon Inactive or no Beacon 

Equator and Equinox designation 

34 0 = Mean equator and equinox 
1 * True equator and equinox 

t 

35 Date of equator and equinox 

1 = Standard equator and equinox 
(Jan. 0.0, 1950) 

_2 = Jah. 0.0 of year of observation 

3 - Equator and equinox of date 
(instant of observation) 

4 ** DODS reference date 

(0 hr. Sept. 18, 19S7) , 

* 

Observation Data 
36-38 R.-A. or H.A. (hours) 

39-40 R.A. of H.A. (minutes of time) 

41-4S R.A. or H.A. {XX. XXX Seconds of time) 

46 Sign of declination 

47-48 Declination (degrees of arc) 

49-50 . Declination (minutes of arc) 

51-54 Declination (XX. XX arc seconds) 

Preprocessing Indicators 

55 0 - Data has been corrected for annual 

aberration effects 

1 e Data not corrected for annual aberration 

56 0 s Data has been corrected for parallactic 

refraction 

1 c Data not corrected for parallactic 
refraction 
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April 8 t 1974 


Columns 

58-61 

62-65 

66 

67-80 


Subset Description 

57 0 *= Data has been corrected for diurnal 

aberration 

1 * Data not corrected for diurnal 
aberration 

Standard deviation in R.A. or H.A, 

XX. XX Arc Seconds multiplied by 
cosine of the declination 

i 

Standard deviation in declination 
XX. XX Arc Seconds 

^Preprocessing Report 

0 = Report not indicated 

1 - 9, A - Z Values to be. assigned. 

* 

Not Used, 
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April 8, 1974 


Columns 33-80 for range and range difference* measurements 
(Data Types 20-29) 

Columns Subset Descript ion 


33-35 


33 


34 


35 


Preprocessing Indicators 

0 = Data has been corrected for ionospheric 

refraction effects. 

1 = Data not corrected for ionospheric 

refraction. 

0 « Data has been corrected for 

tropospheric refraction effects. 

1 = Data not corrected for tropospheric 

refraction. 

2 *= Data has been corrected for tropospheric 

refraction using the correction formulas 
for international laser data (see cols. 
76 - 80 ). 

3 = Data not corrected for tropospheric 

refraction. Columns 76-80 contain 
coefficient for use with international 
laser formulas. 

0 = Data has been corrected for transponder 

delay effects . 

1 = Data not corrected for transponder delay 


36-54 Observation Data 

36-45 Range or Range Difference (Km) 

46-54 Range or Range Difference 

XXX.XXXXXX Meters 


55-56 Preprocessing Indicators 

55 Preprocessing Report 

0 *= Report not indicated 

1 - 9, A - 2 values to be assigned 
*Range Difference - Time Delay VLBI 
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Columns 


57-61 

62-68 

69-73 

74 


75 

76-80 


April 8, 1974 

Subset Description 1 

S6 Transponder Channel for GRARR 

1 = Channel A or 1st sidetone 

2 e Channel B or 2nd sidetone 

3 ■ Channel C or 3rd sidetone 

or 

Transponder type for other systems 

1 - coherent 

2 - non-coherent 

Ref erence .station number for range diff 
erencing or transmitting station number 
for station- to-satellite- to satellite- 
to-station data. 

Relay satellite-ID for satellite- to- 
satellite data. 

Measurement standard deviation 
XX. XXX Meters 

Range ambiguity indicator 
..For DSN S-Band 

1 - MARK 1A (824809582.0 m) 

2 - TAU (151285510.38518 m) 

For USB 

1 - 824809582.0 m 
For GRARR S-Band 

1 - 18737031.3 m 

2 - 4684257.8 m 

3 - 936851.6 m 

Not used 

Tropospheric refraction correction 
XX. XXX Meters 

or, 

Coefficient of Tropospheric Refraction 
international lasers (see col. 34) 

XX. XXX Meters 
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Columns 33-80 for range rate, doppler, and range rate difference* 
measurements (Data types 30-39) 


Columns Subset 


Description 


33-35 

. 33 


34 


35 


Preprocessing Indicators 

0 03 Data has been corrected for ionospheric 
refraction effects. 

- 1 *= Data not corrected for ionospheric 
refraction; 

0 = Data has been corrected for tropospheric 

refraction effects. 

1 = Data not corrected for tropospheric 

refraction. 

‘Receiver Mount Type 

1 = X-Y (East -West) 

2 » X-Y (North-South) 

3 e Azimuth-Elevation 

4. = Hour angle -declinat ion 


56-42 Counting interval for average range rate 

data (Types 30,33,34). 

XXXXX.XX Seconds 

or. 

Tropospheric refraction correction for 
other range rate data types. 

XXXXX.XX cm/Second 


43-54 Observation Data 

43-49 Range rate or range rate difference 

XXXX.XXX Km/Second 

50-55 Range rate or range rate difference 

XXX. XXX Millimeters/Second 

*Range Rate Difference - Fringe Rate VLBI 
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Columns 

56 

57-61 

62-68 

69-73 

74 

75-77 

78-80 


Subset Description 


April o, iy/4 


Preprocess inn Report 

0 » Report not indicated 

1 - 9, A - Z Values to be assigned. 


Reference station number for range 
differencing, or transmitting station 
number for station-to-satellite-to- 
satellite -to -station data or average 
range rate data. 


Relay satellite -ID for satellite -to- 
satellite data. 


Measurement Standard deviation 
XXX. XX milli -meters/second 

Transmitter or reference station Mount 
type 

1 - X-Y (East-West) 

2 = X-Y (North-South) 

3 - Azimuth-Elevation 

4 = Hour angle -declination 

Receiver antenna axis displacement 
XX. X Meters 

Transmitter or reference station 
Antenna axis displacement 
XX. X Meters 
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Columns 

Columns 

33-41 

42-44 

45-47 

i 

48-49 

50-54 

55-58 


Ap i 1 J> ± j * 


33-80 for Altimeter Data (Types 40-44) 


Subset 


Description 


Altimeter Height Measurement (0) 
X, XXX, XXX. XX Meters 

Sigma Measurement 
X.XX Meters 

Tropospheric Refraction 
= 0 - none applied 
f 0 - value applied 
X.XX Meters 

Ionospheric Refraction 
*= 0 - none applied 
f 0 - value applied 
.XX Meters 


50 

51-54 


Geoid Height 
Sign of Height 
Value of Height 
XX. XX Meters 


55 

56-58 


Tide Height 
Sign of Tide 
Value of Tide 
X.XX Meters 
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April o. 


iy /4 


Columns 

59-65 

66-72 

73-77 

78-79 

80 


Subset Pescript ion 

Satellite Latitude* 

59 Sign of Latitude 

60-65 Value of Latitude 

xxexxxx 

Satellite Hast Longitude 
• XXX?XXXX 

Sea Surface Height (0-h) 

Altimeter Height Measurement 
less Satellite Height (h) 
above Spheroid* 

73 Sign of Height 

74-77 Value of Height 

XX. XX Meters 

i . 

Blank 

Preprocessing Report 

0 *■ Report not indicated 

1-9, A-Z = Values to be assigned 
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ApiJLi o, 


iy / 4 


Columns 

33-80 for 

Minitrack Data (Types 50-54) 

Columns 

Subset 

Description 

33-34 

33 

Refraction Indicators 

0 « Data has been corrected for 

ionospheric effects 

1 = Data not corrected for ionosphere 


34 

*0 *= Data has been corrected for 



tropospheric effects 
1 « Data not corrected for troposphere 

35 

> 

Not Used. 

36-54 


Observation Data 


36 

Not Used. 


37-45 

Cosine % 

- 


XXXX.XXXXX Mils 


46-54 

Cosine m 


XXXX.XXXXX Mils 
Not Used. 

Standard Deviation in Cosine £ 

XX. XX Mils 

Standard deviation in Cosine m 
XX. XX Mils 

Preprocessing Report 

0 K Report not indicated. 

1 - 9, A - Z Values to be assigned. 


55-57 

>8-61 

62-65 

66 
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Columns 33 
(Types 60- 

Columns 

33-34 


35 

36-54 


55-57 

58-61 

62-65 


April 8 , 1974 

80 for X-Y angles and Azimuth and Elevation angles 
79) 

Subset Description 

Refraction Indicators 

•33 0 B Data has been corrected for 

ionospheric effects. 

1 « Data not corrected for ionosphere 
34 0 = Data has been corrected for 

tropospheric effects. 

1 *= Data not corrected for troposphere 

Not Used. 

Observation Data 

Azimuth or X angle (degrees) . Sign of 

✓ 

X angle appears in column 36. 

Azimuth or X angle (arc minutes) 

Azimuth or X angle 

XX. XXX Arc Seconds 

Sign of Y angle 
Elevation or Y angle (degrees) 

Elevation or Y angle (minutes) 

Elevation or Y angle 

XX. XX Arc Seconds 

Not Used. 

Standard deviation in X angle or azimuth 
XX. XX Arc Minutes 

Standard deviation in elevation or Y angle 
XX. XX Arc Minutes 


36-38 

39-40 

41-45 

46 

47-48 

49-50 

51-54 


C-26.14 



April 8, 19/4 


Columns 

66 

67-71 

72-76 

77-80 


Subset Description 

Preprocessing Report 

0 = Report not indicated 

1 - 9, A - Z Values to be assigned. 

Tropospheric refraction correction to X 
angle 

XXX. XX Arc Minutes ’ 

Tropospheric refraction correction to 
Y angle or elevation 
XXX. XX Arc Minutes 

Not used. 
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C. 7 RAP DATA FORMAT 


(To Be Determined) 
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C. 8 GROUNDTRACK TAPE FORMAT 

The Groundtrack tape is written by GEODYN specifi- 
cally for input the GROUNDTRACK graphics program. The 
Groundtrack tape format is as follows: 


Columns 

Format 

Description 

1-6 

16 

Date in YYMMDD. 

9-12 

14 

Hours and minutes in HHMM. 

15-21 

F7.4 

Seconds in SS.SSSS. 

24-29 

A6 

Station name 

32-46 

F15.8 

Satellite Geodetic Latitude 
in degrees . 

49-63 

FI 5. 8 

Satellite East Longitude in 
degrees . 

66-80 

F15.5 

Satellite height in meters. 


Total - 80 bytes. 
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C.9 


0RB1 TAPE FORMAT 


An ORB1 tape is a binary satellite ephcmeris tape 
having a logical record of 350 double precision words. 
There are three types of records on the tape - a title 
record, data records, and a sentinel record. 



1 


ORBITAL TAPE FORMAT -1 (TEB) 


WORD 

NUMBER 


BINARY TITLE RECORD FORMAT 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

27 

28 
29 

80 

81 

82 

91 


92 


Form of: data identification = 76796291 
Satellite number 


Date 

Day count of year 
Seconds of Day 
Date 

Day count of year 
Seconds of day 


U.T . 

Start time of satellite 

Data 

U.T. 

End time of satellite 


J Data 

At f interval between satellite data items in seconds 
Date Coordinate system 

Day count of year Reference data time and position 
Apparent sidereal time in radians 

Spares 

C, , atmospheric drag parameter 

d 2 
Area of satellite in cm. 

Mass of satellite in grams 

Lunar perturbations indicator 

(= 1 indicates lunar pert, were used in determining orbit 

*= 0 indicates solar pert, were not used) 

Solar perturbations indicator 

(= 1 indicates solar pert, were used in determining orbit 

* 0 indicates solar pert, were not used) 


101 


epoch time (C.U.T.) 


102 

a, semi-major axis at t 

(C.U.L.) 

103 

e, eccentricity at t 


104 . 

v, true anomaly at t 

(radians) 

105 

X 


satellite position 

106 

y ! 

■ 

vector r at t (C. 

107 

z < 
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f 


WORD 

NUMBER 

108 

109 

110 
111 
112 
114 
.116 

117 

118 
120 
121 
122 

123 

124 

125 

126 

191 

192 

193 
.194 

195 

196 
200 


ORBITAL TAPE FORMAT -1 (TEB) (Cont.)< 


BINARY TITLE RECORD FORMAT 

x ; Satellite velocity 

y Vector r at t Q (C.U.L./C.U.T.) 

z 

r, magnitude of r at t^ (C.U.L.) 
r, magnitude of r_ at t Q (C.U.L./C.U.T.) 

M f mean anomaly at t (radians 

argument of perigee at t Q (radians) 
i, inclination at t (radians) 

fi, right ascension of ascending node at t Q (radians) 
n, mean motion at t (radians/C .U .T . ) 

E, eccentric anomaly at t (radians) 

i, rate of change of argument of perigee at t Q 

(radians/C. U.T.) 

rate of change of R.A. of 
ascending node at t (radians/ C .U .T. ) 

P, period at t Q (C.U.T.) 

Height of perigee at t (C.U.L.) 

Height of apogee at t (C.U.L.) 

Year (last 2 digits) 

Month 
Day 
Hour 
Minute 

Seconds x 1000 

Indicator of theory used to compute orbit 
(* 1 indicates PE orbit generator 
e. 2 indicates Gill Integration orbit generator 
* 3 indicates Brouwer orbit generator) 


Date and 
Time of Epoch 


e 


030 



WORD 

NUMBER 


DATA RECORD 
* 

DESCRIPTION 


1 

2 

3 

4 
.5 
’ 6 

7 

8 
9 

10 

11 


Date (YYMMDD) Time of first 

Day Count of Year (DDD) satellite position 
Seconds of Day (SSSSS) velocity vector - t 
At, interval between position-velocity vectors 
Spare 

x \ Satellite position vector 

y > in km. 


z / 

n 


Satellite velocity vector in 
km/sec . 


z t 


12-305 


49 other satellite position velocity vectors 
at times t+At, T+2At, t+3At,... t+49At 


306-350 Spares 


Notes: If the last data record contains less than 50 * 

position-velocity vectors, the remainder of the 
record is filled with nines (99999999 .DO) 


SENTINEL RECORD 

The first word of the sentinel record contains 
nines (99999999 .ODO) . Words 2-3S0 are spares. 
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C.10 PCE DATA TAPE FORMAT 


A PCE data tape may be written by GEODYN in orbit 
generator mode or read by GEODYN in data reduction mode. 
The format is as follows: 


Columns Format 


Description 


2-3 12 


4-9 

16 

10-13 

14 

14-21 

F8.4 

22-45 

D24. 16 

46-55 

DIO. 4 


Measurement type 

16 - X \ 

17 - Y 

18 - Z 

19 - X 

20 - Y 

21 - Z 


Inertial 

Geocentric 

Cartesian 

Coordinates 


22 - a. \ 

23 - e 

24 - i 

25 - « 

26 - w 

27 - M 



Osculating 

Geocentric 

Keplerian 

Elements 


Date of data (YYMMDD) 


Hours and minutes 


Seconds 


Measurement 

Measurement standard deviation 


Units for measurements and sigmas are: 


X,Y , Z , a 

• • • 

X,Y,Z 

e 

i >U, W »M 


meters (a - m) 
meters/second - cm/sec) 
dimensionless (<r - ppm) 
degrees ( o - arc sec) 
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C* 11 RV TAPE FORMAT 


An RV tape is a satellite ephemeris tape which 
is output by GEODYN specifically for use by the DELTA 
analysis and graphics program. An RV tape is a binary 
tape written without format control. The RV tape word 
format is as follows: 


Word 

IXE£ 

Description 

1 

DP 

Days from Jan 0.0 of the reference 
year for the arc. 

2 

I 

Date in YYMMDD . 

UTC 

3 

I 

Hours and minutes in HHMM. Time 



System, 

4 

R 

Seconds . 

5-10 

DP 

Satellite inertial position-velocity 
vector in meters and meters/second. 

11 

DP 

Satellite latitude in degrees. 

12 

DF 

Satellite longitude in degrees. 

13 

DP 

Satellite height in meters. 


Total = 92 bytes. 

The last record has a value of 999, 0D0 in the 
first word. 
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C. 12 GEODYN SCRATCH FILES 


GEODYN uses seven, scratch files. Some of these 
files are direct access (D.A.), but most are sequential 
(Seq.). The uses of these files are described below: 


File Type 

10 Seq. 

11 Seq. 

12 D. A. 

13 Seq. 

14 Seq. 

15 D. A. 


Seq . 


Description 

Used in the processing of doppler data in 
GEOS format. 

Used in the processing of all GEOS format 
data. 

Used for storage of all data in data 
reduction mode. 

Used for the storage of solar and geo- 
magnetic flux data. 

Used for the storage of a priori information. 

Used for the storage of updated parameter 
values. 

Used for storage of that segment of the 
normal matrix pertaining to individual arcs. 


16 



March 1, 197 3 

C. 13 GEODYN SIMULATED DATA TAPE 

GEODYN has the capability of writing an output tape 
with the calculated observations in place of real observa- 
tions. The tape is a binary unformatted tape designed for 
the purpose of testing the GEODYN estimate process. The 
following format is a nominal format which corresponds to 
internal variable names only. This format should not be 
interpreted rigidly for all measurement types. 


Word 

T ip. e . 

Description 

1 

I 

Date in YYMMDD . 

2 

I 

Time in HHMM . 

3 

D.P. 

/ 

Seconds. 

4 

D.P. 

First calculated measurement. 

S 

D.P. 

Second calcualted measurement. 

6 

D.P. 

Standard deviation of first measurement. 

7 

D.P. 

Standard deviation of second measurement. 

8 

I 

Station number. 

9 

1*2 

Measurement type. 

10 

1*2 

Number of measurements (1 or 2) . 

11 

I 

Satellite ID. 

12 

D.P. 

Tropospheric refraction index. 

13-14 

D.P. 

Transit times for average range rate data. 

15 

1*2 

Transponder channel. 


The measurements and standard deviations are in 
•‘natural" units . 

length - meters 
time - seconds 
angles - radians 

The measurement types are the same as those for the 
SIGMA card. 
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C.14 


SC4020 PLOT TAPE FORMAT 


The SC4020 Plotter Driver Tape contains a contiguous 
stream of 36 bit SC4020 commands ♦ This stream of commands 
is blocked 682 words of length 36 bits each to efficiently 
Utilize the F-53-1A off-line buffer. The particular command 
by-command format of any given Driver Tape is solely a 
-function of the particular information being plotted. 



•APPENDIX D 

BUILT-IN STATION 

STATION* CdCDETIC LATITUDE 


KA ME NU-vfetbM 

UtO 

MN 

SEC GND S 

MIL 4 

l 

2e 

30 

29. 151 A 

G6M3 

2 

2c 

37 

57* 6003 

3CA3 

3 

02 

2 1 

4, 2329 

ANS 3 

A 

i 

0 

59*6273 

CV13 

i 

* 27 

AS 

51* 7824 

ACN3 

o 

-J -57 

— 19* 3126 

MAD 3 

7 

40 

27 

19* 7268 

' C*OT 

* 

-24 “t A 

-23. 93CO 

G*M3 

9 

10 

13 

37*oese 

HS< 4 

10 

-2t -3C 

-0* 24 45 

H*ft 3 

It 

c 2 

7 

34* 0557 

GOSS 


05 

20 

29* 4543 

GV4J 

IS 

27 

*7 

46*731 3 

TEA 3 

w 

27 

2 <3 

12*771 3 

3TC3 

13 

36 

59 

54*4452 

»ioo 

all 

35 

23 

22* 1639 

ECKJ 

.512 

35 

17 

59* it 54 

VEND 

3 14 

3d 

1 A 

51*5255 

MArO 

j14 

3 £ 

25 

33. 1545 

*C»0 

3** l 

-3 1 

-22 

-55.26^0 

TXD3 

«<t£ 

-33 

-?* 

—3*2436 

joj? 

:: 1 

• *i • 

-5 1 

-23.7 = 0 l 

GOOD 

30 1 

<*C 

25 

4 3. 5s 53 

CED3 

OCc 

AO 

? 7 

11. 1359 

WTD 

371 

25 

23 

45* 9716 

<OU«'J*J 

1 <J u 0 

3 

15 

3* 922 1 

JPCINT 

lw Wl 

33 

25 

49* 8C 33 

PTHY"!^ 

UCJ 

26 

32 

‘ 53. 1513 

DLlTD* 

1003 

0 

- J7 

-22. 3462 

LlM\Pi 

10C3 

-11 

— — 3 

- J7. 0634 

SNTAG6 

ICC J 

— 33 

— 3 

-5 im 307 t 

nE*-l I 

1012 

*7 

44 

29* 2567 

COLEG5 

1C 1 J 


52 

17*6614 

GFU»T<3 

1 J 1** 

■* C ‘ 

1 

21.5499 

WNKFL-J 

1-15 

3 1 

25 

U6.4160 

J03U9w 

10 1 3 

-23 

- 53 

- 1.450 e 

*CJ AV£ 

1017 

33 

1 3 

a.7* OQS 4 

CCM£4 A 

10 13 

-J 1 

—23 

-25. £91 6 

ULA5<* 

1- 1 V 

= 4 

id 

3^ • 7775 


1 ‘J £ 1 

3 3 

2 ” 

a.'' 9 p 0 3 6 

j ft*iv»< 

1 V 

2C 

jo 

53* 1512 

MAC jAi 


-l 5 

0 

-32.291 9 

XCC‘-’"*< 

lo2— 

-3 1 

-23 

-25* 5? 1 c 

1CJITJ 

iv«"J 

C 

-37 

-22* 14,2 

ILX»A y 

• wU* 

-l l 

--D 

- yr • 0 5 *<• 4 

IS AT Ao 

t J •• 

-33 


— = a »'jC7t 

is -JAV 

; •* 

2 3 

1 r 

*7* “054 

iuj-su- 

l- - i 

-23 

-53 

— 1* of OP 


f 


OSITIONS 

HAST LONGITUDE SMME»UIO Hctr.HT 


DEG 

4N 

SECONDS 

( MgTEPS J 

279 

13 

23.2243 

-4S.000 

261 

45 

oA. 0263 

-46.000 

255 

20 

3C. 5710 

—66. 000 

2 56 

l % 

50*6524 

-27.000 

3 44 

21 

52.6612 

206.000 

345 

4 0 

20*9220 

555.000 

3E5 

4? 

53.0274 

321 .030 

1 13 

13 

31.21 63 

5.030 

1 44 

44 

14.9202 

76.000 

1 43 

53 

35.'* 110 

1 149. 00J 

2 CO 

20 

5* 6933 

1123.C JO 

243 

7 

34.3299 

920.000 

2 49 

1 6 

43.5135 

-29.000 

2 52 

37 

16. M 68 

-41 .COO 

263 

9 

25. 4569 

-6. 300 

2 43 

9 

1.5593 

933.003 

2 43 

l 1 

34.7211 

936.000 

243 

12 

17. 3501 

1040.030 

243 

6 

37.1409 

979.000 

l 36 

53 

|4. 2974 

143.030 

1 49 

53 

62. 39 10 

673.000 

27 

At 

t. 6213 

1410.030 

355 

4 3 

3.0565 

925.300 

356 

37 

55.3475 

775. 000 

2 79 

26 

1.0750 

-41 .000 

3C7 

1 1 

41.0^07 

-18.000 

262 

34 

45.^303 

—53. 693 

273 

3 

4.2737 

-42. too 

15 J 

25 

16. 4307 

3553. 913 

2 52 

5C 

a9. 1431 

50.733 

2 65 

I * 

53.7002 

713.950 

3C7 

1 D 

46.2670 

46.000 

2 12 

J 

36. 1525 

l 56. 3 67 

2 £2 

5 9 

19.6 1 sO 

203. 1 62 

359 

1 i 

5.0626 

90.410 

27 

42 

26.21 53 

1540.977 

243 

5 

59.0170 

576.254 

l It 

52 

15. 1 J63 

130.403 

2 t2 2 5 

30*0077 

263. 125 

262 

34' 

40.7303 

—6 2. t>3 3 

. 2 76 d 

4.2737 

-42. 000 

’ *7 

1 7 

£9 • 3065 

1359.777 

1 26 

5 2 

15. 1 963 

130.403 

2 51 

25 

15.4 167 

3555. 913 

2 l’2 

60 

59.1-31 

53.7C3 

259 

19 

£2.7502 

71 J . 6 ; j 

243 

5 

59. 01 70 

176 . 2 5«* 

27 

42 

26.2153 

1540. *77 


reproducibility of the 
ORIGINAL PAGE IS POOR 
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itATiow CcUCSTIC LATITUDE 


NA’4^ NU.v-Js.rt 

JtC 

mn % 

SECONDS 


1032 

«*7 

44 

29. 2e67 

1C0LEG 

U-3 

C >4 , 

S2 

i7.e?i a 

IGFO^K 

l J i** 

-* £ 

1 

21. f *5? 

t *n<Fl 

1335 

j i 

26 

46.41 fcO 

1ULA3K 

1C io 

C A 

t a 

36.7775 

1UJSMN 

•10j7 

JL 

12 

7.29< 3 

io^ipu 

U jd 

-32 

— 27 

-J2. 656? 

tfcJOMA 

10** 2 

5 = 

12 

7. 3057 

ITALIAN 

10AJ 

-IS 

0 

-32.3515 

OdJAj 

tUl 

-3 = 

-27 

— 32. 656 5 

.MAC G 4.7 S 

1 1 

— IS 

-1 

— 1 A • 47 n 1 

aosiANS 

11 £C 

Js 

1 1 

45. 6 55 C 

SMTAGA s 

1127 

—3 2 

w » 

-4. 2CS A 

AUASKP3 

1 i c 5 

3 ** 

5? 

19*051 l 

CA^VCN S 

* I =2 

-24 

— 5— 

-11.4592 

MACSAWV 

1=22 

- 1 > 

-1 

-17. 2^74 

CCSA AMV 

1 £ 

J = 

1 1 

41«4$4C 

SNTASAV 

imf 

— J j 

•4 

-6.773s 

alakpv 

lg£j 

z A 


17. 3FS5 

-AOVJSV 

ia = 2 

-24 

— ~Z*~ 

-16.4155 

ANT IGA 

2 i (Jq 

1 7 

3 

55. 1626 

GFNVL5 

5 323 

J 2 

23 

*,9.4714 

GPVILL 

£j3h 

32 

25 

32. 4554 

•jS'fa; 

2-, 0 3 

2 9 


22.6137 

SSDF7D 

3-21 

42 

27 

17.521 f 

SEMM£3 

2 -,0 2 

50 

— 6 

SO* C553 

SWANt 3 


17 

24 

1». 1053 

Gkw T 

£+03 

2 1 

25 

4S. 66 l £ 

Cl/ 5 4C3 

2-* Cc 

1 i 

5 

25. 4435 

TFNDAJ 

2*07 

1C 

*4 

35. 7464 

OF A 

.w51 

A? 

i 5 

3b. 7360 

T»l iC< 

2- 5 2 

3c 

7 

25. 54«9 

PCTmGJ 

2 + oJ 

5 l 

24 

57. 5143 

ATHSG7 

2 ♦ & J 

37 

5 3 

26. 18 2 7 

TC*MO : * 

J— c ♦ 

•w 0 

2 9 

14. 436« 

CHOF'JJ 


35 

-*0 

** 4723 

kindly 

2471 

3 2 

22 

£B, *3471 

£Ca ADS 

> tj. 

3-4 

b7 

aO • j 1 1 4 

hi.n r V- 

2 3 •« J 

52 

0 

f . « 512 5 

jt;P i 

2- -> i 

27 

4 

4 

15.601 £ 

Al.ST-7 ON 

Jj a 7 

J * 

> m 

15. 100 1 

HC^EST 


2a 

!0 

25. S?5c 

<z*-y *v* 

;<u2 

■4 1 

7 

69. 1 364 

5T301- 

*0 •» 3 

2 o 

25 

5*. i oyc 

i T^Z A< 

---1 

- = 

3 

= 3.7 M3 

Z T ' 

-.5.3 


- :•- 

-3-. lift 


-.-i >3 

5 e 

i j 

J & . 7 f. 5 2 

4N T 

. J a 1 

1 7 


37. 5Ct J 


D-2 


EAST LONGITUDE SPHENOID HEIGHT 


3 EG 

MM 

SECONDS 

1METEW3 1 

3C7 

16 

46. 2670 

46.000 

2 12 

9 

36.5525 

156.367 

2 €2 

5 9 

19.5160 

203. 1 32 

359 

13 

3.0325 

90*410 

2 12 

23 

30.6377 

233.125 

277 

.7 

41.27 32 

549.7J3 

1 AS 

57 

14. 9123 

949*573 

277 

7 

40.9741 

350.061 

47 

17 
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Altimeter Data 

Arc Description Cards 

Arc Set 

Area of Satellite 
Atmospheric Density 

Average Range Rate Input 

/ 

Biases 

Electronic 
Measurement Biases 
Standard Deviation 
Station Timing Biases 


Page 

C-22 ,23,24,25,26 
1 . 2 - 1, 2 

1. 0- 1; 1. 2-1 ;Section 1.2 
1.2-46,47 

1 . 1 - 1 
C-17-26 

1.1- ll;1.2-26,27 

1.2- 12,13,14 

1.2- 13 

1.2- 12,13,14 


Binary Residual Tape 
B-Matrix 
B-Matrix Tape 
Built-in Stations 
Central Body 
Common Set 

Constant Timing Corrections 
Constrained Stations 
Convergence Criterion 
Coordinate Systems. 

Core Requirements 
Correlation Coefficients 
DATA Card 
Data Reduction 


1.2-36 ; 2 . 0-8 ; C-6 

1.1- 5;1. 1-5.1 

1.2- 15,16;2.0-9;C-2,3,4,5 
Appendix D 

1.1- 5. 4 

1. 0- 1:1. l-l;Section 1.1 

1.2- 38,39 

1.1- 7;1 .1-25 

1.1- ll;l.l-37 

1 . 2 - 4 : 1 . 2 - 8 ; 1 . 2 - 10 ; 1 . 2-37 

2 . 0 - 2, 3 

1 . 1 - 6, 7 

See Termination Cards 

i f ii; 1.2 -34; 1.2 -36, 37; 1.2 -53, 54; 

1.2- 65 , 66 ; 2 . 0 -.7 ; 3 . 0 - 2 , 3 ; 4 . 0-6 , 7 , 8 ; 

4. 0-9, 10, 11; 4. 0-12, 13; 4. 0-14, IS, 16 


1-1 




poor 



' AUgUSt J!., IV'/ ‘S 


I 

INDEX (Cont.) 


c 



Page 

Data Tape Output 

Appendix 

Deletion of Data Types 

• 

DODS Format 

• 

1.2-67 

GEOS Format 

1.2-70 i 

Denormalization Equation 

1. 1-17, 18 ; 1. 2-42, 43 

Diagnostic Error Messages 

See Error Messages 

DODS Data Base 

1.2-6;2.0-6;S.0-l,2,3; 
Appendix A 

DODS Data 

1.2-S;1.2-50,52;1.2-66;1.2-67; 
2 . 0-8 ; C-7 , 8 , 9 , 10 

Drag 

See Force Model 

Earth 

See Force Model 

Earth Constants 

• 1.1-8,9,10 

Editing 

l.l-ll;See Measurement Editing 

Electronic Biases 

See Biases -Electronic 

Elevation Angle 

l.l-12;1.2-30 

Ephemeris Output 

1. 2-33, 34 ;1.2-3S, 36,37 ; 2. 0-6,7 

Ephemeris Tape 

5.0-34,35;C-ll,12 

Epoch Time 

1. 2-1 ,2;1, 2-S ,6,7 

Error Messages 

Section 5.0 

Diagnostic Errors- 
General 

S. 0-4,5, 6, 7, 8, 9 

DODS Data Base 

5. 0-1, 2, 3 

Ephemeris Tape 
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Example Jobs 

B-Matrix 

Data Reduction 
* 

DELTA 

DODS Data Base 
GEORGE 
GROUNDTRACK 
0RB1 Conversion 
Orbit Generator 
PCE Data Tape 
Flux Data 
Force Model 
Drag 
Earth 

Jupiter 

Mars 

Moon 

** Saturn 

Solar Radiation 
Pressure 

Su. 

Tides 

Venus 


* 

Pace 

Section 4.0 ' 

4.0- 14,13,16 

4. 0- 23, 24, 25; 4. 0-6, 7, 8; 4. 0-9, 10,11; 

4. 0- 12, 13;4. 0-14, 15, 16 

4.0- 9,10,11 

4.0- 12,13 

4. 0- 6, 7, 8, 9 

4. 0- 2, 3, 4, 5 

4.0- 9,10,11 

4.0- 9,10,ll;4.0-17,18 

4.0- 17,18 

1.1- 13,14;2.0-11 

1 . 1 - 1 

1.2-16;1.2-21, 22;1.2-46,47 

1.1- 8, 9, 10; 1.1- 14. 1,14. -2, 14. 3, 14. 4, 

14.5.14.6.1.1- lS,16,17,18;l.l-34; 

1.2-16;1.2-23,24,25 ; 1.2 -40, 41, 42, 43 

1.1- 4,4.1:1.2-17,18 

1.1- 4,4.1;1.2-17,1S 

1.1- 4.4.1;1.2-17,18 

1.1- 4,4.1:1.2-17,18 

1.2- 16; 1.2-46,47;1. 2 -SS, 56 

1.1- 4,4.1;1.2-17,18 

1.1- 31,32 

1.1- 4,4.1;1.2-17,18 
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Formats 

Format Specification 
• arid Code • 

I/O File* Formats 
Option Cards 

* "FORTRAN Logical Units 
GEODYN 

Geopotential Coefficients 

Geopotential Model 
GEOS Data 

Gravity Perturbations 

GROUNDTRACK Tape 

Hardware and Software 
Restrictions 

Inner Iterations 

- Input Cards 

Integration 

Iteration Convergence Test 
Iteration RMS 
Iterations 


Page 


1 . 3-1 

Appendix C;See also Table of 
Contents 

1 . 1 - 3 ; 1 . 2-11 
*5ee JCL 

i pii 

1.1- l;See Force Model -Earth; 

1 . 2 - 16 

See Force Model -Earth 

1 . 2 - 5 ; 1 . 2-50, 52 ; 1 . 2-65,66; 

1 . 2 - 68 , 69, 70; 2. 0-8 ;C- 13, 14, IS, 16 


See Force Model 


1 . 2-36 ; 2. 0-8 ; C-27 


2.0-11 

See Iterations 

1 . 0 - 1 ; 2 . 0-6 


1.1- 26,26.1,27;l.l-33,34; 

1 . 2 - 57,5S,59;1.2,61,61.1 


1.2-7; 


1 


.2-43. 


l.l-ll»See Convergence Criterion 
See RMS ' 


ii; 1.1-5; 1.1 -11; 1.1-19, 20; 1.1-37; 

1. 2-5, 7 ; 1 . 2-16 jl. 2-34 
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Job Control Language (JCL) 
Core Requirements 
Examples 

FORTRAN Logical Units 
Overlay 
Job Submittal 
JPL Ephemeris 
Jupiter 

Keplerian Element Ephemeris 

Logical Units 

Lunar Perturbations 

Magnetic Flux 

Mandatory Cards 

Mars 

Mass of Satellite 
Mass Ratios 
Master Blocks 
Master -Stations 
Measurements 

Deletion 

Editing 

Partials 

Standard Deviations 

Two Station Measure- 
ments 


Page 

Section 2.1 

2. 0- 2,3 

2. 0- 4,5 

2. 0- 6, 7, 8,9 

2 . 0 - 10 
Section 3.0 
2 , 0-6 

See Force Model 
1. 2-36 ; 2 . 0- 7 

See JCL-FORTRAN Logical Units 

1.1- l;See Force Model-Moon 

1.1- 13,14 

1.0-1 

See Force Model 

1.2- 46,47 

1. 1- 5;1 . 2-18 

1.1- 29,30 

1.1- 7 ;l.i 25 ; 1. 2 - 16 




1.2-37 


See Standard Deviations -Measurement 


1.2-51 
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Minitrack Data 
Moon : . 

Optical Data Corrections 
Option Cards 
Optional Cards 
Orbital Elements 


C-16 

See Force Model 

1.2-39,40 

1.1- 1;1.1-3;1.2-1;1.2-11 

1 . 0-1 

1.2- 1,2;1.2-8,9,10;1,2-16 
See also Position Vector and 
Velocity Vector 


Orbital Tape Format 


C-29,30,31 


Orbit Generator 


i,ii; 1.2 -33, 34; 1.2 -35, 36, 37; 1.2 -54; 

1.2-65;3.0-2,3;4.0-9,10,11;4.0-17,1S 


0RB1 Tape 
Order of Arcs 
Output Specifications 
Outer Iterations 
Overlay Structure 
PCE Data 

Planetary Gravity Effects 
Polar Motion Tables 
Position Vector 
Preprocessing of Data 


1.2-31, 32;2.0-9;C-28 

1.1- 23 

See Printout Requests 
See Iterations 
2 . 0-10 

1.2- 5; 1.2 -50, 52; 1.2 -54; 1.2 -66; 

2.0- 8 ;C-32 

1.1- l;See Force Model 

2 . 0-11 

1. 2- 2 ;1. 2-8 ,9 

1.2- 38,39,40 
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ORIGINAL PAGE' IS POOR 
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Printout Requests 

Reference Time 
Refraction Index 
^Residual Printout 
Restart 

Right Ascension of Greenwich 
RMS , 

Run Description Cards 
RV Tape 

Satellite-to -Satellite Data 
Saturn 

Scratch Files - GEODYN 
Program 

SC4020 Plot Tape 

Selection/Deletion Cards 

Selection of Data Types - 
GEOS Format 

Set-up 

Sigma Defaults 

Simulated Data Tape - 
GEODYN Program 

Software Restrictions 

Solar Flux 

Solar Perturbations 

Solar Radiation Pressure 


Page 

1.1- 19, 20; 1; 2-31, 32; 1.2 -33, 34; 

1.2- 35, 36, 37; 1.2-44 ,4 5; 2. 0-6, 7 

1.2- l,2;1.2-4 

1.2- 39 

1. 1- 19, 20;1. 2-44,45 

1 . 1 - 21,22 

2.0-11 

1. 1- 11;1. 2-29 

1 . 1 - 1, 2 

1.2- 33, 34;2.0-9;C-33 
C- 2 2, 23, 24, 25, 26 
See Force Model 

C-34 

t-36 

1. 2- 1,1.2-65; 1.2 -67, 68 ,69,7 0,71 

1.2- 68,69,70 
Section 3.0 

1.2- 50,51 

1.2- 5;1. 2-53,54 ;1. 2-66 ; 2.0-8 ;C-55 
2.0-11 

1.1-13,14 

See Force Model -Sun 
See Force Model 
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Biases 
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C and S Coefficients 

Drag Coefficients 

Measurement 

Solar Radiation Pres- 
sure Coefficients 

Station Positions 

Surface Density 

Station Adjustment Sigmas 

Station Coordinate Adjust- 
ments 

Station Coordinate Cards 

Station Coordinates 
Station Estimation 
Station Positions 

Station Position Estimation 

Step Size 

Sun 

Surface Densities 
Termination Cards 
Common Set 
Arc Sc- 

Op lion Cards 

* Selection/Dele- 

tion Cards 


1*1-16, 17;1. 2-41, 42 

1 . 2-22 

1.2- 28 ,29; 1* 2-48,49 ,50, 51, 52 

1.2- 56 

1.1- 24,25 

1.1- 28,30 

1.1- 35 

See Correlation Coefficients 

1 • 1-1 ; 1 . 1 -3 5 , 36 , 36 . 1 ; See Station 
Positions 

See Station Positions 

See Station Position Estimation 

1.1- l;l.l-6,7 ;1. 1-24,2 5; 

1 .1- 35 ,36 ,36 . 1 ; Appendix D 

1.1- 24,25 

See Integration 
See Force Model 

1. 1- 1 ;1. 1-1; 1.1-27, 28, 29, 30 

1.1- 1;1. 1-37 

1.2- l;1.2-65,66 

1.2- 71 



Terininum Elements 

Tidal Perturbations 

Time Difference Tables 

Timing Corrections 

Title Cards 

Tracking Complement 

Transit Time Correction 

Tropospheric Refraction 
Correction 

Two-Body Earth Computations 
Two-Station Measurements 
Unconstrained Stations 
Variance-Covariance Matrix 
Velocity Vector 
Venus 

VLBI Format 


A U 8 U 


See Force Model-Tides 
2 . 0-11 
1.2-38,39 

1.1- 2;1. 2-1,3 

’See Station Positions 

1.2- 58 . 

1.2- 38 

1.2- 25 

1.2- 51 '■ 

1. 1- 7, 1.1-25;1. 2-16 

1.2- 62,63,64 ' 

1 . 2 - 2 ; 1 . 2-10 
See Force Model 
C-17-26 
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